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CHAPTER I 


INTRODUCTION 

Since the launch of the first commercial communications satellite, 
the INTELSAT I, in April 1965, satellite communications have become more 
and more popular internationally [1,2,3]. Not only does the number of 
satellites continue to increase, but the demand for future satellite 
services is growing rapidly [3,4]. This situation creates a problem of 
how to provide enough satellite communications capacity to satisfy all 
the potential users. The study of this report is devoted to 
contributions toward solving this problem by developing methods that can 
efficiently utilize the spectrum and orbit resources. 

The most popular satellite orbit for the civil communications 
services is the geostationary orbit [5,6]. The idea of using this orbit 
was proposed by Arthur C. Clarke [1]. However, the geostationary orbit 
can accommodate only a limited number of satellites for a given 
frequency channel because satellites that have the same frequency 
channel must be properly separated from one another in space for 
acceptable interference protection [6], When the geostationary orbit is 
considered "crowded" with satellites, it is crowded in terms of 
electromagnetic compatibility. This requirement greatly limits the 
satellite capacity of the geostationary orbit. 

A method of increasing the communications capacity is to increase 
the spectral band available for satellite communications use. The 
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spectral region of interest lies between the maximum usable frequency 
(MUF) for reflection by the ionosphere and the first oxygen absorption 
line (about 60 GHz) [7]. By international agreement, this spectral 
region is divided into many bands for various services, e.g., the 
broadcasting-satellite service (BSS), the fixed-satellite service (FSS), 
the military service, navigation, weather detection, etc [2,8]. As a 
result, only a limited number of spectral bands are available for the 
civil communications services. 

As an example of the large demand for satellite communications, 
when the 11 to 12 GHz spectral band for civil satellite service was 
opened for use on an international basis, every administration requested 
some frequency channels [4]. Even those administrations that do not 
need, or cannot afford, a satellite at the present time requested 
channel assignments for future use. Therefore, lack of enough spectrum 
and orbit resources becomes a serious problem with regard to planning 
for future communications satellite traffic, i.e., orbital and frequency 
allocations. 

The goal of the orbital and frequency planning task is to use the 
limited amount of resources to provide enough communications capacity to 
satisfy every potential user. The main concern of the planning task is 
to ensure that mutual interference between different satellite systems 
is acceptable; a poorly planned scenario would be likely to result in 
unacceptable carrier-to-interference ratios (C/I), and hence 
unacceptable signal quality, for at least some users. To avoid this, 
the satellite orbital locations and frequency channel allocations should 
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be carefully planned. This planning task becomes very difficult when 
the number of satellites and channels to be assigned becomes large. 

Moreover, since the demand for present and future communications 
satellites is large, it is important to use the available spectrum and 
orbit resources efficiently. In other words, the objective is to 
achieve the maximum information-transfer capacity for the resources 
allocated. 

In 1977, the spectral band 11.7 to 12.5 GHz was assigned to Region 
1 and the band 11.7 to 12.2 GHz to regions 2 and 3 for the planning of 
the BSS. The administrations in Regions 1 and 3 completed the planning 
in the 1977 World Administrative Radio Conference (WARC-77) [9]. The 
plan was based generally on spacing the satellites uniformly at 6 
degrees, when they are assigned the same frequency channel [10]. The 
administrations in Region 2 delayed the planning process in WARC-77, 
with the intention of making the most efficient use of the 
geostationary orbit and the spectral band [11,12]; then proposed a plan 
at the Regional Administrative Radio Conference in 1983 (RARC-83) 
[13,14]. Planning for the fixed-satellite communications service (FSS) 
has been deferred to conferences in 1988 or beyond. Therefore, a method 
is needed to solve the assignment problem in as near an optimal way with 
regard to this objective as possible. 

Much work has been done on this problem. Several studies discuss 
some important factors that should be considered in the assignment 
problem. A study by D. J. Withers identifies three areas that should be 
exercised in order to achieve effective resource utilization [15]. The 


3 


first area is engineering for an interference-limited environment by 
properly using the antenna characteristics, signal modulation, and 
multiple-access techniques to minimize the interference power. The 
second area is effective inter-system coordination. It includes: 1) 

proper pairing of the up- and down-link bands, 2) a standard frequency 
translation between up- and down-link bands, 3) an agreed scale of 
permissible single-entry interference noise allocation, 4) maximization 
of the satellite service arc by easing the satellite elevation angle 
constraint. The third area is the reduction of Inhomogeneity in 
orbit-spectrum sharing. Also in the Final Acts of the WARC-77, the 
importance of placing satellites as close as possible for efficient 
orbit utilization is expressed [16]: the satellite spacing should be 

small, while still keeping the mutual interference acceptable. 

Some studies deal with the orbital assignment problem alone. A 
Japanese study tackles the problem of orbit utilization through a 
non-linear programming optimization procedure with the objective of 
minimizing the total orbital arc used for a scenario [17], The basic 
approach is to relate the satellite geocentric angular separation to the 
interference power; the problem is formulated as a non-linear 
programming problem and the sequential unconstrained minimization 
technique is applied to solve it numerically. Another Japanese study 
modified the above program so that it can find the optimal orbital 
assignment for a new satellite when It is Inserted into an existing 
scenario [18]. The result is optimal in the sense of finding the best 
location for the new satellite, while making some modification to the 
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existing scenario, such that the satellite separations or, equivalently, 
the interference powers still meet the requirements. The best location 
for inserting the new satellite is chosen so that the total orbital arc 
occupied by the final result is minimized. To carry out rigid planning 
of many satellites, an evolutional model is used by repeatedly using the 
modified program to assign satellite locations one by one, and the 
ordering of insertion must first be chosen. It is not obvious, and 
appears unlikely, that the quality of the result is independent of this 
ordering. 

Some studies deal with the frequency assignment problem alone. The 
methods of map-coloring and dot-linkage have been proposed to achieve 
the most conservational use of the spectrum resource [19,20,21,22]. A 
Japanese study tackles this problem by rearranging the frequency 
assignment of a given scenario that has an initial optimal orbital 
assignment [23]; e.g., orbital assignment is obtained from [17] with the 
frequency assignments assumed the same for all satellites. Therefore, 
in this Japanese study the frequency assignment, identified as a 
permutation problem, is handled independently of the orbital assignment; 
it is formulated as an integer program and the optimal permutation is 
obtained via the branch-and-bound method. When the new frequency 
assignments are made, the required satellite orbital separations may be 
reduced. The objective is the minimization of the total orbital arc 
occupied by the final scenario. Note that the combination of the two 
Japanese studies [17,23] becomes a complete package that solves the 
orbital and frequency assignments in two steps. 
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Some studies deal with both the orbital and frequency assignments. 
Some of these are aimed more at obtaining a better understanding of the 
problem than at actually solving it [24,25]. An application-oriented 
research group in Canada developed a software package for planning 
synthesis in connection with RARC-83 [26,27,28]. It is a multiple-stage 
process, and one of the objectives spelled out in the report is the 
minimization of the total orbital arc occupied by the scenario. First 
the minimum required orbital separations that meet the single-entry 
protection ratio are calculated for all pairs of satellites. These data 
are calculated in the following four cases: co-polarization and 

co-channel, co-polarization and adjacent channel, cross-polarization and 
co-channel, cross-polarization and adjacent channel. The initial 
scenario is an equal-spaced orbital assignment with co-channel, 
co-polarization frequency assignment. Then, the computer program can 
let the planner make changes both manually and automatically in the 
channel and/or polarization assignments, and manually in the orbital 
assignments. The program always produces a scenario, even if it turns 
out not to meet the required protection ratios. The result is a local 
optimum and not necessarily a global optimum. 

Two methods have been proposed by a research group at the Ohio 
State University. The first method uses an extended gradient search 
technique to improve an existing scenario [29]. A sum of negative 
exponentials of the aggregate effective C/I ratios is used as the 
objective function to be minimized; therefore the procedure seeks to 
maximize the smallest of all such C/I ratios. The gradient direction of 
the objective function at the point representing the existing scenario 
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is calculated; then the objective function values at some discrete 
points along the negative of this direction are found. The point 
yielding the most favorable (minimum) objective function value is chosen 
as a new scenario to start another search process. The iterative 
process stops when a better scenario can not be found. The second 
method is called the cyclic coordinate search method [30], In this 
method each orbital and frequency variable is varied in turn. Each time 
a set of points is examined, the point yielding the most favorable 
objective function value is selected as the new coordinate value for 
that variable. A cycle is completed when all the orbital and frequency 
variables have been varied once. The cyclic process is repeated with 
suitably adjusted step sizes until it reaches a solution where the 
improvement of C/I results halts. The detailed description of the two 
methods is given in Chapter III. 

Transmission of a signal from the Earth terminal and its reception 
at a satellite constitutes an up-link; and transmission of a signal from 
a satellite and its reception at the Earth terminal constitutes a 
down-link. By international agreement, the up-link signals and the 
down-link signals are not in the same spectral band, so that the signals 
will not interfere with each other [31]. In this study only the 
regulation of the down-link traffic is considered. The up-link problem 
can be implemented similarly in another spectral band. Note that the 
above studies [17,26,27,28,29,30] also deal only with the down-link 
communications traffic regulation. 

In Chapter II, some of the important factors and parameters 
involved in the carrier and interference power calculations are 
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discussed. Note that in this study an effective carrier-to-interference 
power ratio is used to evaluate the feasibility of a scenario; in 
Section II. F it will be shown that this is equivalent to the more usual 
margin calculation. 

In Chapter III, the objective function to be minimized in the 
extended-gradient and cyclic-coordinate search techniques to improve a 
given scenario is analyzed [29]. It is shown that this function has 
large values when two (or more) orbital or frequency assignments are 
collocated; this indicates that both the signal quality requirement and 
the permutation of the orbital /frequency assignments need to be looked 
into in order to find the globally optimal scenario. Furthermore, some 
arguments are given which indicate that the objective function is likely 
to be a function with only one local minimum for a fixed permutation of 
orbital and frequency assignments; even though a definitive proof has 
eluded us, a set of numerical examples are presented which support these 
arguments. This suggests that a sufficient condition for obtaining the 
globally optimal solution by the extended-gradient search method is that 
it should terminate as an ordinary gradient search with the optimal 
permutation in orbital and frequency assignments. 

In Chapter IV, a different approach, the AS concept, is presented. 
It is shown that the single-entry C/I protection requirement is 
equivalent to a required minimum satellite separation; hence the highly 
non-linear C/I requirement can be viewed as constraints on the satellite 
locations. With this approach the orbital assignment can be formulated 
as a mixed-integer program and solved by the branch-and-bound method; or 
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it can be formulated as a linear program with non-linear side 
constraints and solved by a version of the simplex method for linear 
programming with restricted basis entry. 

Chapter V intends to show that an important body of information for 
choosing FSS service areas is the communications -demand density. It is 
proposed that the service areas of an FSS system should be specified 
according to the communications-demand density in conjunction with the 
concept of small -beam design; the frequency re-use scheme can be 
implemented through small beams and well -separated service areas. A 
case study demonstrates that the communications supply for the United 
States could be significantly increased if the service areas are 
specified according to these principles. The AS concept presented in 
Chapter IV is used in this case study. 
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CHAPTER II 


DESCRIPTION OF PARAMETERS AND FACTORS IN C/I CALCULATIONS 
A. INTRODUCTION 

The feasibility of a scenario is evaluated according to the signal 
quality [32], which is usually expressed in terms of the signal -to-noise 
ratio (S/N) [33], For instance, the unweighted signal quality 
requirement for 625-line, color-television signal is suggested to be a 
S/N ratio of 33 dB for 99% of the worst month [34]. The signal power 
(S) is measured in the baseband channel after modulation improvement and 
baseband processing [35,36]. For the purpose of planning the 
broadcasting-satellite service (BSS), the requirement is that the 
pre-detection carrier-to-noise ratio (C/N) at the receiver input should 
equal or exceed 14 dB for 99% of the worst month or 10 dB for 99.9% of 
the worst month [37,38]. 

The noise power includes the receiver thermal noise (N^) and all 
the interference powers from other communications systems [39]. In 
order that the interference powers from other communications systems 
will not further degrade the C/N level, the total (or aggregate) 
interference power level should be weak enough so that its contribution 
to the total noise power is negligible. 
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There are several criteria how weak the total interference power 
should be [40,41]; two which are commonly used are the interference-to- 
thermal noise ratio (I/N^) and the carrier-to-interference power ratio 
(C/I) [40]. In this study the C/I criterion is used to evaluate the 
feasibility of a scenario, with a minor modification. Account must be 
taken of the fact that a given interference level at the same frequency 
as the desired carrier is more damaging than the same level of 
interference at a far removed frequency. This can be done by 
multiplying each interference power at a non-carrier frequency by a 
relative protection factor (less than unity) before adding the 
interference powers to obtain an effective ratio (C/I e ). This is 
discussed in more detail in Section II. F, where it is also shown that 
this formulation is fully equivalent to the more usual, but less 
convenient, representation in terms of margins. 

In the following sections some of the important factors and 
parameters involved in the C/I e calculation are discussed. 

6. SATELLITE 

A satellite is a relay terminal, its basic function is to receive 
signals from some Earth stations and to re-transmit them to other Earth 
stations. It has receiving and transmitting antennas, and a frequency 
translator [42,43] or, in some cases, more sophisticated signal - 
processing circuits which include frequency shifting. 

A satellite may use any orbit to travel [2], but the geostationary 
orbit is the most popular orbit for civil communications. In this orbit 
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the satellite position is almost stationary with respect to any place on 
Earth; thus complicated satellite tracking systems at the Earth 
terminals may be avoided. This orbit is directly above the equator at a 
height of 37,165 kilometers above the Earth surface [2], By 
international agreement this will be the main orbit for civil 
communications services [5,6]. Therefore, in this study only this orbit 
is considered. 

Because the satellite is away from the Earth, it is possible to 
control its position and attitude only to a certain precision. A 
satellite drifts away from its designated orbital location and its 
position needs to be adjusted from time to time [44,45], With present 
technology, a satellite orbital location may be maintained within 0.1 
degree in the north-south and east-west directions, resulting in 0.14 
degree of maximum excursion [44]. As for attitude control, the 
transmitting antenna pointing error may be kept within 0.1 degree, and 
the tolerance in the rotation about the beam axis is typically two 
degrees [44,45]. For the calculations of this report, the pointing 
error may be taken into account through the minimum elliptical beam 
calculation to be discussed in Section II.D, or through the antenna 
discrimination function to be discussed in Section II. E. The satellite 
rotational error may be taken into account through the minimum 
elliptical beam calculation to be discussed in Section II.D. 
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C. SERVICE AREA 


A service area is a designated area on the Earth surface to which a 
corresponding satellite directs its signals. Actually an antenna 
transmits its signal in all directions according to its pattern 
function, even though its target area is of limited size. By 
international agreement, a service area should be illuminated by its 
satellite main beam within the -3 dB contour of that beam [46,47]. This 
means that the received power density at any point in the service area 
should be within 3 dB of the power density at the antenna aim point. 

For simplicity, a service area is represented by a set of test 
points at its boundary. To evaluate the C/I e results of a scenario, one 
calculates the C/I e values at all the test points of all the service 
areas, and compares them with the C/I e requirement level. Because the 
interference is likely to be the worst on the service area boundaries, a 
scenario with satisfactory C/I e values at all the test points should 
guarantee that the C/I e values will be good at all the places inside all 
the service areas. The test points are also used to generate the 
minimum elliptical beam parameters to be discussed in Section II. 0. 

In this study a satellite and its corresponding service area(s) are 
viewed as a communications system. A service area may have several 
satellites, but all of these satellites are treated independently of one 
another. 
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0. MINIMUM ELLIPTICAL BEAM 


The antenna main beam (-3 dB contour) of a satellite to its 
designated service area should be shaped, using available antenna 
technology, to fit the shape of the service area as suitably as possible 
[48]. However, a reference antenna pattern for shaped beams is still 
unavailable, i.e., there exists no internationally agreed method for 
predicting, for regulatory purposes, the discrimination to be expected 
from shaped beams outside their respective service areas. In order to 
analyze a scenario in terms of C/I e results, the interference powers are 
therefore calculated by assuming that a service area is illuminated by a 
minimum elliptical beam from the satellite position that covers all the 
test points of that service area [47]. 

This minimum ellipse is specified by five parameters: the 

longitude and latitude of the beam aim point, the orientation angle, and 
the major- and minor-axis beam widths. In Figure 2.1, the aim point, 
which is on the Earth surface, is denoted as A. The antenna beam plane 
is a plane perpendicular to the satellite beam axis, SA, and passing 
through the aim point A. The vectors AM and AN denote the major and 
minor axes. The beam widths in the directions of these two axes, viewed 
from the satellite, are <|fo ax and i^i n , respectively. The orientation 
angle, which is not shown, is the angle measured anti -clockwise in the 
antenna beam plane from a line parallel to the equatorial plane to the 
major axis of the ellipse. The more detailed descriptions of these 
parameters can be found in a National Tel ecommuni cation and Information 
Administration (NTIA) document [49]. 
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Figure 2.1. Minimum ellipse configuration. 
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To generate the ellipse data, one needs to specify the nominal 
satellite location and the locations of a set of test points; other 
parameters that can be specified are the minimum beam width, minimum 
elevation angle, satellite pointing error and rotational error. The 
resulting elliptical cone covers all the test points that are used to 
generate this ellipse; in a worst-case calculation that includes the 
specified satellite pointing error and rotational error. A computer 
program developed by Akima was used to generate the five ellipse 
parameters for the calculations of this report [49]. 

With these parameters, the half -power beam width (HPBW), ^ 0 » with 
respect to any test point, T, on the Earth surface can be calculated by 
means of a somewhat involved procedure. First the vector AP of Figure 
2.1 is calculated, where P is the intersection of the line ST and the 
antenna beam plane. Then, the angle a between the vector AP and the 
major axis AM is calculated. The length of AR is calculated from 

AR = AM-[cos 2 (a) + A^sin 2 ^)]" 1 ' 2 , (2.1) 

where A^, called the axial ratio of the elliptical beam, is the ratio 
between AM and AN. Then it is assumed that the ratios between the 
angles t max » t to and their corresponding arcs AM, AN, AR are the 

same; this is a very good approximation when the angles are not large. 
With this approximation, the angle Is calculated from 

2 22 - 1/2 

^to = ^max*^ cos + \ »sin (a)] , (2.2) 
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where A can be seen as the ratio between >1/ and <i> . . 
r r max r min 

The off-axis angle, i|^, toward the point T is calculated by the 
cosine law 

cos(^ t ) = (ST 2 + SA 2 - AT 2 )/(2* ST* SA) , (2.3) 

where ST, SA, AT denote distances in Figure 2.1. 

The off-axis angle and the corresponding HPBW will be used 
to calculate the directivity of the antenna toward the specific point, 

T. The antenna pattern envelopes used to calculate the directivity will 
be described in Section II. E. 

The detailed description of all the parameters and the formulations 
of all the calculations can be found in the Spectrum Orbit Utilization 
Program (SOUP) manual [50]. A streamlined SOUP code is listed in 
Appendix A. It is much less complicated than the original SOUP code. 

The basic calculation is exactly the same, but with fewer options. For 
instance, the streamlined version does not consider propagation loss. 

E. ANTENNA REFERENCE PATTERNS AND PROTECTION RATIO 

Two sets of antenna reference patterns are used in the study. They 
are representations, adopted by the International Radio Consultative 
Committee (CCIR), of the envelopes of real antenna patterns. Their use 
therefore should result in a near worst-case interference power 
calculation. The first set includes the satellite transmitting antenna 
reference patterns and the ground receiving antenna reference patterns 
suggested in the International Radio Consultative Committee Conference 
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Preparatory Meeting in 1982 (CCIR-CPM-82) for the Regional 
Administration Radio Conference for Region 2 in 1983 (RARC-83) 
[10,51,52]. These are the BSS patterns. The gain, G, of the 
transmitting antenna is calculated from [53,54] 

G = e-[0r/* max )/(180/223)] 2 - A p , (2.4) 

where the symbols mean 

e : the beam efficiency, taken as 0.6 in this study, 

'Wnax : the beam width of the major axis, in radians, 

A r : the axial ratio of the elliptical beam. 

The ground receiving antenna gain in the BSS calculations is taken as 
40.2 dB, corresponding to a circular-beam antenna of 1-meter diameter, 

12 GHz carrier frequency and 0.6 beam efficiency. The reference 
patterns are shown in Figures 2.2 and 2.3, and are used in the 
calculation of the objective function values in Chapter III. 

In Figures 2.2 and 2.3, both the transmitting and receiving 
antennas have two reference patterns with mutually orthogonal 
polarizations. The transmitting antenna transmits signals of the 
designated polarization according to the co-polarization pattern, at the 
same time it also transmits orthogonally polarized signals according to 
the cross-polarization pattern. The receiving antenna receives the 
wanted signals according to the co-polarization pattern by aligning its 
polarization parallel to that of the wanted signal. Any signal that is 
co-pol ari zed with the wanted signal is received according to the 
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A: co-pol ar component (dB) 

-12(<l7*o) 2 

-30 

-17. 5-25- 1og(4»/i); 0 ) 

after intersection with curve 

B: cross-polar component (dB) 

-30 

-40-40* 1 ogC(^/^ 9 )-D 

after intersection with curve 

C: minus the on-axis gain (dB), 

: off-axis angle 

^ 0 : half-power beam width 


for 0 < ip/\p 0 < 1.58 

for 1.58 < i|»/i|/ 0 < 3.16 

for 3.16 < ty/ipQ 

C: as curve C 


for 0 < < 1.56 

for 1.56 < \p/ip 0 
C: as curve C 

G=42.5 dB in this illustration) 


Figure 2.2. BSS satellite transmitting reference patterns. 


19 





0.1 0.2 0.5 1 2 5 10 20 50 100 


o 


co-pol ar component without 

side-lobe suppression i 

(dB) 

0 

for 

’I'/'I'O < 

0.25 

-12(M 2 ) 2 
-11.3-25* log(ip/ip) 

for 

0.25 < < 

0.94 

for 

0.94 < < 

18.88 

-43.2dB (-3dBi ) 

for 

18.88 < t|// ip 0 


cross-polar component (dB) 

-25 

for 

0 < < 

0.25 

-30-40* logCl-log(^/ip )] 

for 

0.25 < Tp/4> 0 < 

0.44 

-20 

for 

0.44 < < 

1.28 

-17 . 3-25 *1 og (\p/ip 0 ) 

for 

1.28 < il>/ik 0 ^ 

3.22 

-30 until intersection with co-polar component curve; 

then as for co-polar component 



^ : off-axis angle, 
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Half-power beam width = 1.8 degrees 
Nominal on axis gain = 40.2 dB 

Flat portion for ip/^ 0 less than 0.25 takes account the antenna 
pointing error 


Figure 2.3. BSS ground receiving reference patterns. 
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co-polarization pattern, and any signal that is cross-polarized with the 
wanted signal is received according to the cross-polarization pattern. 
These unwanted signals constitute the interference power at the receiver 
input. 

In the second set of reference antenna patterns, intended for the 
fixed-satellite service (FSS) calculations, the satellite transmitting 
antenna reference pattern given in Figure 2.4 is a modified version of 
the fast fall -off reference pattern from CCIR-CPM-82/RARC-83 [551. This 
modified reference pattern had been suggested as useful for the FSS 
systems [56,571. The transmitting antenna gain is again calculated from 
Equation (2.4). The ground receiving antenna reference pattern, shown 
in Figure 2.5, is a modified version of the reference pattern from the 
International Radio Consultative Committee (CCIR-82) Report 391-4 [581. 
The modification is suggested in a CCIR-82 Recommendation and in the 
CCIR-CPM-82/RARC-83 Report [59,601, and should become the standard in 
the year 1987. The receiving antenna gain is 43.2 dB, corresponding to 
a circular-beam antenna of 4.5-meter diameter, 4 GHz carrier frequency 
and 0.6 beam efficiency. Note that there is no cross-polarization 
pattern in Figure 2.4 because such a pattern has not been adopted by the 
international committee. 

The relative protection ratio, PR(dB)-PR 0 (dB) , used in Chapter III 
is taken from a CCIR-CPM-82 Report [611 and is shown in Figure 2.6. The 
value of PR 0 is the co-channel protection ratio; the value of the actual 
protection ratio PR(dR)-PR 0 (dB) depends on the carrier frequency offset 
between the wanted and unwanted signals as well as on the modulation 
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Figure 2.4. FSS satellite transmitting reference pattern. 
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Figure 2.5. FSS ground receiving reference pattern. 
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Relative protection ratio 
(PR-PRo) (dB) 



method of the' signals [62], The Carson bandwidth will be taken as 25.2 
MHz, which is appropriate for the BSS TV/FM case. 

F. RECEIVED-POWER CALCULATION 

Refer to Figure 2.7 for the geometry, but note that the radius of 
the Earth relative to that of the geostationary orbit is exaggerated 
greatly for clarity. In this figure, the satellites assigned to service 
areas A, B are designated S^, Sr, respectively. The aim points of S/\, 

Sr are the points a, b. The point d is one of the test points in A. 

The test points are normally chosen on the boundary of their service 
areas because interference is likely to be the worst there. The minimum 
ellipses of A, B from S^, Sr are labeled E(A), E(B). For the test point 

d, the off-axis angle of the S ^ signal and the corresponding HPBW in the 

direction toward d are ^tc» Hco* the subscript c is meant as a mnemonic 
for carrier. The distance from S A to d is x. These values are used to 
calculate the carrier power received from S ^ at d. Also, for test point 
d, the off-axis angle of the S R signal and the corresponding HPBW are 

designated ttio» respectively. Since the receiving antenna at the 

test point d is pointed at S/\, its off-axis angle toward Sr is tri • The 
distance from Sr to d is y. These values are used to calculate the 
interference power received from Sr at d. 

Referring to Figure 2.7, the carrier power, C, at d is calculated 
by means of the Friis transmission equation [63,64]: 
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Figure 2.7. Configuration of received-power calculation. 
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(2.5) 


W D A^tc)* G d * c2 
' ' ^-(4») 2 -x 2 
where the symbols mean 

PA : the transmitting power of Sa, 

Ga : the transmitting antenna gain of Sa, 

Da : the co-polarization transmitting antenna discrimination 
from Sa in the direction toward d, 

G(j : the receiving antenna gain at d, 
c : the velocity of light, 
fA : the carrier frequency. 


In Chapter III, the effective isotropic radiated powers are assumed 
constant for all the satellites. In Chapters IV and V, it is assumed 
that the carrier power densities at the aim points are equal for all the 
service areas [65]; therefore, all the satellite transmitting powers are 
adjusted to meet this requirement. 

The interference power received at d must be calculated with care 
because there may be a polarization mismatch between the wanted and 
unwanted signals. Both of the orthogonally polarized signals 
transmitted from Sg must be decomposed into two components, one that is 
parallel to the wanted signal of the receiver and the other that Is 
orthogonal. With the proper choice of the antenna reference patterns, 
the received interference power from each component is calculated from 


I 


P B* G B ,D B^ti ),G d* D d(Vi )* c2 
(4tt ) 2 -y 2 


( 2 . 6 ) 


27 



where the symbols mean 

PB : the transmitting power of Sb, 

Gb : the transmitting antenna gain of Sb, 

Ob : the transmitting discrimination from Sb to d, 

G(j : the receiving antenna gain at d, 

Od : the receiving antenna discrimination from d in the 

direction toward Sb when the antenna is pointed at S/y, 
fB : the interference frequency. 


The total interference power is the summation of all the components. 
For instance when the signals of S/\ and Sb are co-polarized, the 
received interference power at d from Sb is calculated from [64,66] 


*cp ^c.rc + *tx,rx + ^tc,rx + *tx,rc^*°x * 


(2.7a) 


when they are cross-polarized, the received interference power is 
calculated from 


I 

xp 


I. + I. + 
tc,rx tx,rc 


^tc,rc + *tx,rx^* D x 


(2.7b) 


Here the subscripts mean 

cp : wanted and unwanted signals are co-polarized, 
xp : wanted and unwanted signals are cross-polarized, 
tc : signal transmitted according to the co-polarized reference 
pattern, 

tx : signal transmitted according to the cross-polarized 
reference pattern, 

rc : signal received according to the co-polarized reference 
pattern, * 
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rx : signal received according to the cross-polarized reference 
pattern. 


The term D x is the rain depolarization factor n66»67]. Since the 
propagation effect is not considered in this report, the term n x is 
taken to be zero. Also, Equation (2.7a) is further simplified as 


I 


cp 


^c.rc 


(2.7c) 


because the cross-polarization component, which arises from the cross- 
polarized patterns of both the transmitting and receiving antennas, is 
negligible compared to the co-polarization component. Finally it should 
be noted that Equation (2.7b) is an approximation. In an actual case, 
these terms should add as phasors, not in a power sense; but to perform 
that calculation the relative phases of the tc, tx patterns and that of 
the rc, rx patterns would have to be known. So, without the term D x , a 
worst-case formula would be 

*xp ~ ^ *tc,rx + *tx,rc^ * (2.7d) 

When the carrier frequencies of the wanted and unwanted signals are 
different, a frequency filtering factor must be included in the 
effective interference power calculation. The proper expression for 
such a filtering factor can be easily obtained as follows. When 
evaluating the effect of the unwanted signal, the usual procedure is to 
calculate a term called the protection margin from 
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( 2 . 8 ) 


M.(dB) = (C/I.)(dB) - PR.(dB), 
where the symbols mean 

i : index of frequency channel of the wanted signal, 

Ii : total interference power in channel i, 

PR-j : protection ratio against interference power in channel 
i. 

The interference power is acceptable when the value of M-j(dB) is 
positive. (When there are several interference signals of the same 
carrier frequency, their total interference power is obtained by 
decomposing every polarized signal into parallel and orthogonal 
components with respect to the wanted signal, then calculating each 
interference power and summing them, as discussed above.) When there 
are several unwanted signals of different carrier frequencies, the 
equivalent protection margin, M, that evaluates the over-all effect of 
the interference power is calculated from [68,69] 

-Mi (dB)/10 

M(dB) = -10* log [ I 10 n ], (2.9) 

i 

where each term M^ is the protection margin in one frequency channel, 
and the summation is over all frequency channels. Equation (2.9) can be 
re-written as 
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M(dR) = -l0.1og(zl0- {(C/I i ) W- PR iW }/10 ) 


-10 .log ( E 10 - { ( C / I i)( dR )-i:PRi( dB )-PRo(dR)]>/10 #10 PRo(dB)/10 


-10.1og( E > (^)-CPRi (dB)-PR o ( d B)3}/10 _ pR 

i o v 

, , Ui(dB)+[PRi(dB)-PR o (dB)}/10 

C(dB) - 10*log( E 10 1V 1V ov ' ) - PR (dB) 


( 2 . 10 ) 


If the equivalent total, or aggregate effective interference power, I e , 
is defined from the expression of the equivalent protection margin, M, 
as 

M(dB) = (C/I e )(dB) - PR 0 (dB) 


» C(dB) - Ie (dB) - PR 0 (dB) , (2.11) 

then the equivalent aggregate interference power can be expressed as 

I (dB) . 10.10 9 ( * 1 o tI i (dB)+CPRl (dB)-PR o (dB)3>/10 ) ^ 

6 i 

Therefore, the quantity [PR.. (dB)-PR Q (dB)], which is the relative 
protection ratio in Figure 2.6, can be Interpreted as a filtering factor 
operating on the interference power when the carrier frequencies of the 
wanted and unwanted signals are different, and denoted as 

F ^wanted ’^unwanted ^ dB ^’ ThiS a PP roach wil1 be taken throughout 
Chapter III. For a scenario of many satellites and where each satellite 

has many frequency channels, the aggregate effective interference power. 
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I® n j, received in channel n at test point j in service area k is 
calculated from a power summation of terms of the type given in Equation 
( 2 . 12 ), 


I® n j(dB) = 10* log{ z Z 10 


[Iim,knj(dB)+F(f n ,f m )(dB)]/10 


ieK/k meN. 


Z 10 


[Ikm.knjW + Fffn.fnHdBN/W 


me N^/n 


(2.13a) 


or 


T knj 


ieK/k meN. F(f n’ f m ) + me u /n ^m.knj F( W 


(2.13b) 


where the symbols mean 


K/k : the index set of all the satellites, excluding 
satellite S^, 

Ni : the index set of all the channels assigned to 

satellite i , 

lim.knj : the single-entry interference power from channel m of 
S-j , received at channel n of test point j in service 
area k, 

F (fn»fm) : the relative protection ratio between the carrier 

frequency f n and interference frequency f m , as shown 
in Figure 2.6. 


The aggregate (C/I e ) value in a channel at test point d is obtained 
by calculating the values of C using Equation (2.5) and I e using 
Equation (2.13b) and dividing. 
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G. SIGNAL QUALITY REQUIREMENT 


When the signal quality requirement is stated as "S/N be no less 
than 30 dR for 99 % of the worst month" [381, the main concern is to 
obtain satisfactory performance except during rare, large attenuation of 
signals due to heavy rains [701. The prediction of statistical rain 
rate distributions as a function of geography will not be adressed in 
this study. In this study the carrier and interference powers are 
calculated without considering atmospheric absorption and rain 
attenuation. As for the assumption that the carrier power densities are 
the same at all satellite antenna aim points, the satellite transmitting 
powers can be adjusted to allow for these propagation effects so that 
the power density requirements are still satisfied; but this has not 
been done in the calculations which will be presented. 

Most satellites function as a repeater [711: they receive a 

signal, change the signal carrier frequency, and transmit it back to 
Earth. In such a design, any interference power generated on the 
up-link remains in the signal when re-transmitted in the down-link [721. 
In this study, only the interference power generated in the down-link is 
considered; any interference power from the up-link is not included in 
the calculation. For BSS, it has been proposed in international 
telecommunication meetings that the C/I e requirement value in the 
up-link be 10 dB higher (better) than that of the down-link [73,74,751. 
Therefore for the BSS the interference power in the up-link should make 
an insignificant contribution to the overall interference power, and can 
be neglected. For FSS, the C/I e requirement values for the up- and 
down-links may not differ much [76,771. Therefore for the FSS case the 
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C/I e requirement value for each half link must be specified with care so 
that the total interference power does not degrade the signal quality 
excessively. 

In evaluating the feasibility of a scenario, first one calculates 
the C/I e values at all the test points; then compares them with the 
pre-determined C/I e requirement value. A scenario is said to be good 
when all the resulting C/I e values exceed the requirement value [78,791. 
In Chapter III, a C/I e requirement value of 30 dR will be used to 
evaluate the quality of the extended gradient search results; the 
precise required value depends on the modulations, but 30 dB is typical 
[78]. In Chapter IV, the requirement value is arbitrarily set at 25 dB 
in the numerical calculation; the precise value is relatively 
unimportant because the main goal of the numerical calculations in 
Chapter IV is to demonstrate the feasibility of the method. In Chapter 
V, the requirement value is set at 20 dB in the numerical calculation, 
where the U.S. is used as a case study, for the following reason. For 
the U.S. domestic satellite systems in the 6/4 GHz band, four-degree 
spacing between two satellites was used initially to regulate the 
orbital assignments. However, in order to better utilize the orbit 
resource, the Federal Communication Commission (FCC) has decided to 
adopt two-degree spacing in the satellite planning task in the 6/4 GHz 
band, to be effective beginning in 1987 [801. Since the method proposed 
in Chapter V is not meant only for the U.S., the C/I e requirement is set 
at 20 dB in the numerical calculations so that it could appear more 
reasonable on the international basis, and would be compatible with the 
assumed antenna parameters, as discussed in Section V.E.l. 
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CHAPTER III 


THE OBJECTIVE FUNCTION OF THE EXTENDED-GRADIENT AND 
CYCLIC-COORDINATE SEARCH METHODS 


A. INTRODUCTION 

The main purpose of the study in this chapter is to explore the 
mathematical nature of the orbital /frequency assignment problem by 
investigating an objective function used in the extended-gradient and 
cyclic-coordinate search methods to solve the satellite planning 
problem. It is shown that in order to find the global optimum solution 
one must deal with both the permutation of the orbital /frequency 
assignments and the signal quality requirement, e.g., the C/I ratio, in 
the optimization process. It is also shown that for a given set of 
scenarios with fixed orbital permutation (frequencies fixed) the 
function defined by the smallest single-entry C/I value has at most one 
local maximum. This strongly suggests that the objective function has 
only one local minimum for a given orbital permutation, as supported 
by some numerical evidence; this indicates that the extended gradient 
search process is highly likely to find the global optimal solution if 
it terminates as an ordinary gradient search with optimal 
orbital /frequency permutation. 
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B. GENERAL DISCUSSION OF THE GRADIENT AND CYCLIC-COORDINATE SEARCH 
METHODS AND THEIR OBJECTIVE FUNCTIONS 

1. Introduction 

The purpose of including this section is to introduce the basic 
principles of the gradient and cyclic-coordinate search methods, and to 
demonstrate how they perform for a particular objective function. Using 
the gradient search method for the assignment problem was proposed by 
Professor Clarence H. Martin of the Department of Industrial and Systems 
Engineering [29]; using the cyclic coordinate search method was proposed 
by Professor Clark A. Mount-Campbell of the Department of Industrial and 
Systems Engineering [30]; the objective function used in both methods 
was formulated by Professor Curt A. Levis of the Department of 
Electrical Engineering and Professor Clarence H. Martin of the 
Department of Industrial and Systems Engineering [29]; the numerical 
calculations discussed in Section III.B were performed by Professor 
Charles H. Reilly and Mr. David J. Gonsalvez of the Department of 
Industrial and Systems Engineering [30]. 

2. The Gradient and Cyclic-Coordinate Search Methods 

The gradient and cyclic-coordinate search methods are two 
techniques commonly used by systems engineers to find an optimal 
condition for system performance [81,82]. In applying either method, an 
objective function, which is a function of a set of decision variables, 
is constructed in order to rank candidate solutions. The decision 
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variables represent the controllable operating conditions of the system; 
in the present application, they are the satellite orbital location and 
frequency channel assignments. The objective function should be a 
measure of the performance of the system under the conditions specified 
by the values of the decision variables. For instance, one solution, or 
specification of values for the decision variables, should represent 
more attractive operating conditions than a second solution if the 
former solution yields a greater (smaller) value when evaluated in an 
objective function which is to be maximized (minimized). The optimal 
solution (operating conditions) should be that solution which provides 
the greatest (smallest) objective-function value, when the objective 
function is to be maximized (minimized). 

With such an objective function, the gradient search is performed 
in the following way [29, Rl]. Assume that the function is concave 
(convex) and is to be maximized (minimized). First a point that 
represents an initial operating condition is located, and the gradient 
components of the objective function at that point are calculated. Then 
a proper step size is chosen in the gradient (negative-gradient) 
direction so as to reach another point representing another operating 
condition. Because the objective-function value predicts the system 
performance, a sufficiently small step in the gradient (negative- 
gradient) direction should always lead to a new condition that is better 
than the initial condition, at least according to the chosen objective 
function. Then, the same procedure is repeated with the improved 
operating condition as the new initial point and the search step size 
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properly adjusted according to the magnitude of the gradient: the 

smaller the gradient, the smaller the step size. This procedure should 
be repeated until a point is reached at which the gradient search does 
not yield improvement in system performance. This final solution will 
correspond to the optimal operating condition. 

When the objective function is not concave (convex), global 
optimality of the solution is not guaranteed; instead the solution may 
converge to a local optimum, with the choice of the initial starting 
point influencing strongly which local optimum is selected. To reduce 
this influence and enhance the likelihood of finding a global or near- 
global optimum, a modified gradient search procedure, called an extended 
gradient serach, was used in the present application [86], Consider the 
objective function as one to be minimized. First, the gradient 
direction at the initial point is calculated. Then from this point a 
search line is extended in the negative-gradient direction to the 
boundary of the feasible region; the objective function is calculated at 
a set of ten equally-spaced points along that search line. The point 
with the most favorable objective-function value is chosen as the new 
starting point to do another calculation. If the starting solution is 
the most favorable solution, another ten equally-spaced points between 
the starting solution and the first tested point are examined, the point 
yielding the most favorable objective-function value is chosen as the 
new starting solution. The procedure stops when no solution can be 
found that yields a more favorable objective-function value. 
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Since the search line is extended, the method may allow the search 
process to go from a region with one local minimum to another region; 
this reduces the likelihood of the search being "trapped" at an 
undesirable minimum, i.e., one much greater than the global minimum. 

The cyclic coordinate search method, with a given objective 
function, proceeds as follows [30,82], First, an initial operating 
condition is assumed. Then the decision variables are varied in turn, 
one at a time. Each time a set of points, where the corresponding 
variable is varied along its feasible coordinate range, is examined; the 
point yielding the most favorable objective-function value is 
identified, and the search continues from this new solution with another 
decision variable and the same search process. A cycle is completed 
when every decision variable has been allowed to vary; once a cycle is 
completed, another cycle can begin. When there is no more improvement 
in the objective-function value, the process is repeated within a 
smaller region and with smaller step-sizes. The whole search process 
terminates when no more improvement is obtained with a step size 
commensurate with the accuracy to be obtained. 

3. An Objective Function for Satellite Orbital /Frequency 
Assi gnments 

For the satellite planning problem, a suitable objective function 
is formulated as [29] 
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Z = E E E Z|< n j 
keK neN|< jeJ|< 

= E E E exp (a-[C. .(dB)-I® .(dB)]} , (3.1) 

keK neN k jeJ k knj KnJ 

where the symbols mean 

K : the index set of all satellites, 

N|< : the index set of all frequency channels assigned 

to satellite k, 

J|< : the index set of test points of service area k, 

a : arbitrary parameter used to avoid overflow in computer 
calculation, 

Cknj : the carrier power at channel n of test point j in 

service area k, in dB, 
e 

Iknj • the effective interference power at channel n of test 
point j in service area k, in dB. 

It is assumed that each satellite is associated with one service area; 

thus the index k that represents a satellite also represents its 

corresponding service area. Note that it is the satellites that are 

counted in the index set K; one service area may be served by several 

satellites and each satellite is treated as an individual unit. The 

e 

value of Ck n j is calculated from Equation (2.5). The value of I^nj 1s 
calculated from Equation (2.13a), reproduced here 


40 


I 


e 

knj 


(dB) 


ID* 1 og, o { E E 10 [I< Mnj(dB)* F(f n, f m)(dBn /10 t 
ieK/k meN.. 


i o nkm,knj(d B ) +F (^n»^m)( dR )]/l°j 

meN^/n 


where I-j m> knj> calculated from Equation (2.6), is the interference from 
channel m of satellite i into channel n of test point j in a service 
area served by satellite k. In Equation (3.1) the exponentially 
weighted summation is over all the frequency channels, at all the test 
points, in all the service areas. In Equation (3.2), the double- 
summation adds the interference from all other satellites, and the 
single-summation adds the interference from other channels of the same 
satellite. 

The satellite locations, denoted by o-j for satellite 1, and the 
carrier frequencies, denoted by f -j ^ for channel h of satellite 1, are 
the decision variables. The limits of the orbital variables are usually 
determined by elevation angle constraints [83]; sometimes an eclipse- 
protection requirement may impose additional restrictions [84], The 
limits of the frequency variables are determined by the limits of the 
available spectral band [4]. 

It is clear that the orbital variables are continuous variables; 
however, this is not necessarily true for the frequency variables. In 
past international conferences, generally the available spectral bands 
were each divided into channels of equal bandwidth, with each channel 
specified by its center frequency [85], This is likely to be true for 
future conferences also; then the frequency variables will be discrete. 


41 



In the gradient search method, the derivatives with respect to all the 
decision variables need to be calculated; since it is impossible to 
differentiate with respect to a discrete variable, it is assumed here 
that in the frequency assignment the center frequencies of all the 
channels are allowed to vary continuously while each channel still has 
the same, fixed bandwidth. With this assumption, differentiation with 
respect to the frequency variables is allowed, and the frequency 
protection ratio of Figure 2.6 can be used. It is hoped, but cannot be 
guaranteed, that the optimization of the continuous-varying channel 
problem will lead to at least a near-optimum of the discrete channel 
problem. 

Note in Figure 2.6 that there is a plateau in the frequency 
protection ratio. When two satellites are assigned frequencies with the 
frequency offset in this range, the gradient search method would find 
variation of these frequency assignments not useful because it would not 
change the objective-function value. This is definitely not the result 
the system planner wants, because separating the frequency assignments 
sufficiently could produce better C/I e results. To avoid this problem, 
the plateau is deformed to form an isosceles triangle with small slopes 
(±0.05 dB/unit 8); this modification allows the frequency assignments to 
be separated if the frequency offset is located in the plateau region. 

From Equation (3.1) it is clear that the value of Zknj is small for 
large (C/I e )knj( clB ) values. Since a good scenario should have large 
C/I e at all test points, the global minimum of the objective function is 
likely to be a good solution for the assignment problem. 
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The gradient components of Z are 


VZ 


a a a 

( aoj z * ao 2 z ’*’* af n 


(3.3) 


where the partial derivatives have been taken for all the decision 
variables; also 


T- ~ kcK neN^ 3^1^ 


"knj 


•3( C / I e )| tn j (dB) 
ao. 


af 


3 I E I i . - 9(C/I e»knj< dB » 

Z = kcK ncN k JcJ k ""•> 


lh 


lh 


(3.4) 


(3.5) 


Referring to Figure 2.5 and Equations (2.5), (2.6), moving clearly 
changes ^ c , <i> n - , and x, moving Sr changes also parameters , \p r ^ , and 
y. Every satellite acts both as a desired and interfering source, so it 
affects the parameters associated with S ^ in some terms, and those with 
Sr in others, and both in a few terms. As one moves S/\ and Sr, the 
effective isotropic radiated powers, i.e., the products P/\*G/\, Pr*Gr 
of satellites and Sr respectively, are kept constant; but the carrier 
powers calculated from Equation (2.5) change somewhat at the test 
points, causing some contribution to 9Z/3o-|. Also, the minimum ellipses 
must be recomputed for new satellite locations; the change in ellipse 
changes ^tco » Hio an d therefore D ^ (^tc)» B B (tti) since \)/\ depends on 
ttc/'f'tco and similarly for Dr. 

Here one sees that in Equation (3.4) or (3.5), Z| <n j is the 
weighting factor for the terms 3/30] [(C/I e )| <n j(dB)] and 
3/3f-|h[(C/I e )knj( dB )]* since the term z knj is a negative exponential 
function of (C/I e )knj( dB ^ n wil1 be lar 9 est for the values k, n, and j 
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for which (C/I e )|< n j (dB) is smallest; so the weight serves to emphasize 
the contribution from those test points and channels which need 
improvement most urgently. With such an objective function, the 
gradient search method tends to relocate most strongly the 
orbital /frequency assignments that are responsible for the worst C/I e 
terms, and the result is an increase of these C/I e values. This, of 
course, was the rationale for choosing the exponential function. A 
numerical example, below, will illustrate this point. 

4. Numerical Exercise Using the Extended Gradient Search Method 

A numerical exercise will now be given to show that the properties 
associated with the objective function are indeed as discussed, and that 
the performance of the extended gradient search method with respect to 
this objective function is as predicted. 

In this exercise seven administrations are under consideration; 
they are Argentina, Bolivia, Brazil, Chile, Paraguay, Peru and Uruguay 
(denoted as ARG, BOL, BRZ, CHL, PRG, PRU, and URG respectively) with the 
geographic relation shown in Figure 3.1 and the test points listed in 
Table 3.1. It is assumed that every administration has requested one 
orbital location and three contiguous frequency channels, with cross 
polarizations for adjacent channels, for its satellite; the orbital 
locations and the carrier frequencies of the leading (lowest) channels 
are the decision variables, and these leading channels all have the same 
polarization. The feasible orbital arc was taken from 90 to 110 degrees 
west, and the spectral band from 12,200 to 12,300 MHz, for all 
satellites; the bandwidth of each channel is assumed to be 12 MHz, with 
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Figure 3.1. 


Geographic relation of the seven South American 
administrations. 
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Table 3.1 


■ 

Test Points of Seven Administrations 


ARGENTINA PARAGUAY 


LON . 

LAT. 

LON. 

LAT. 

-66.2 

-21.8 

-57.6 

-25.3 

-62.8 

-22.0 

-58.6 

-27.3 

-53.8 

-27.2 

-56.2 

-27.2 

-56.7 

-36.9 

-54.7 

-25.5 

-63.8 

-54.7 

-54.2 

-24.1 

-68.3 

-54.8 

-58.1 

-20.2 

-73.2 

-50.9 

-59.1 

-19.3 

-71 .4 

-39.0 

-62.2 

-20.5 

-70.5 

-31.4 

-62.7 

-22.2 

-68.6 

-24 . 8 

-58.7 

-27.2 


BOLIVIA 


PERU 


-65 .0 

-12.2 

-70.4 

-18.3 

-65.5 

-9.8 

-69.0 

-12.3 

-69 .0 

-11.2 

-70.5 

-9.4 

-60.0 

-16. 1 

-74.0 

-7.6 

-57.5 

-18.0 

-70.0 

-2.7 

-67.5 

-22 . 7 

-75.2 

-0.0 



-80.3 

-3.4 



-81.3 

-4.4 

BRAZIL 


-81.2 

-6.1 



-76 . 1 

-13.4 

-60.5 

4.5 



-52.0 

3.0 



-46.0 

-1 .5 

URUGUAY 


-35.0 

-7.5 



-42 .0 

-22.5 

-56.2 

-34.9 

-53.0 

-32.5 

-54.9 

-35.0 

-56.3 

-29.5 

-53.5 

-34.0 

-70.8 

-10.5 

-53.2 

-32.7 

-73 .0 

-7.0 

-55.6 

-30.8 

-69.0 

1.0 

-56.9 

-30. 1 



-57.6 

-30.2 



-58.2 

-31.9 

CHILE 


-58.4 

-33.9 



-57.9 

-34.5 

-69.5 

-17.5 



-67.1 

-23.0 



-70.0 

-34.2 



-71 .7 

-43.2 



-68.4 

-52.3 



-72.8 

-51.3 



-75.7 

-46.8 



-74.0 

-28.9 



-70.4 

-18.3 
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a 2.58 MHz guardband between two channels. The carrier and interference 
powers are calculated from Equations (2.5), (2.6) and (2.13a); the 
antenna reference patterns are from Figures 2.2 and 2.3; the frequency 
protection ratio is from Figure 2.6; the ellipse data are calculated 
from [49], 

In this calculation the initial scenario is that all satellites are 
collocated at 110 degrees west, and the initial frequency assignments 
extend from 12,235 to 12,265 MHz in 5 MHz intervals. The intermediate 
solutions at all the search steps and the final solution after ten 
search processes are shown in Figure 3.2 to demonstrate how the search 
process proceeds. In Figure 3.2(c), only the worst aggregate C/I e 
values for these administrations are shown. In these figures the 
administrations are denoted by numbers according to alphabetic order: 1 

for Argentina, 2 for Bolivia, 3 for Brazil, 4 for Chile, 5 for Paraguay, 
6 for Peru and 7 for Uruguay. The improvement of the C/I e results is 
clearly seen as the iteration process proceeds. The orbital assignments 
are almost steady after iteration 7, and the frequency assignments are 
almost steady after iteration 5; this indicates that the search process 
may have reached the vicinity of a local minimum. 

The improvement mechanism of the extended-gradient search method 
can be observed as follows. Note that the C/I e results of 
administrations 3, 4, 6 and 7 are the worst after iteration 1; then at 
iteration 2, the orbital and frequency assignments of these four 
satellites make a very significant change, while that of the other three 
satellites are almost unchanged. This is exactly the purpose of the 
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extended gradient search method 


FREQUENCY (MHZ) 

.00 12220.00 12240.00 12260.00 12280.00 12300.00 



(b) Frequency assignment variation 


Figure 3.2. Continued. 
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exponentiation in the objective function as discussed in Section 
III.B.3. 

The purpose of modifying the ordinary gradient search method can be 
seen as follows. It is known that the objective function takes on large 
values when some of the orbital and/or frequency assignments are 
identical (for reason to be explained in Section III .C) ; hence 
exchanging the position of two satellites means jumping over such a 
region. In Figure 3.2(a) the orbital order established in iteration 1 
is disturbed by the extended search process: the order of satellites 2 

and 5 is changed at iteration 7. Thus the extended gradient search 
method can move from a region with one local minimum to another region. 

5. Numerical Exercise Using the Cyclic Coordinate Search Method 

The performance of the cyclic coordinate search method is 
demonstrated here; the same objective function is used to solve the same 
assignment problem. The initial scenario is changed as follows: all 

satellites are at 110 degrees west and 12,250 MHz. The results are 
shown in Figure 3.3. The improvement in the C/I e results is obvious. 

C. EMPIRICAL EXAMINATION OF THE OBJECTIVE -FUNCTION TOPOGRAPHY 

1. The Importance of the Objective-Function Surface Topography to 
the Gradient and Cyclic-Coordinate Search Methods 

The topography of the objective function greatly influences the 
outcome of the search methods. An objective function has only one 
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FREQUENCY (MHZ) 

12200.00 12220.00 122110.00 12260.00 12280.00 12300.00 



0 12 3 4 


ITER 

(b) Frequency assigm 



local maximum (minimum) if it is strictly concave (convex) or 
pseudo -concave (pseudo-convex). (See Appendix B for the definitions and 
properties of concave, strictly-concave and pseudo-concave functions.) 
The gradient and cyclic-coordinate search methods should lead one to a 
point very close to this maximum (minimum), if not directly locate it. 

If the objective function is not strictly concave, strictly convex, 
pseudo-concave, or pseudo-convex, it may have several local maxima or 
minima. The greatest maximum (least minimum) is called the global 
maximum (minimum), or simply maximum (minimum). The best operating 
conditions correspond to the global maximum (minimum), when the 
objective function is to be maximized (minimized). For such a function, 
both the gradient and the cyclic-coordinate search procedures will 
almost certainly give improvement to an initial operating conditions. 
However, both procedures may eventually be trapped at a local optimum 
instead of reaching the global optimum. 

Therefore, knowledge of the mathematical properties of the 
objective function is important in determining whether the gradient and 
cyclic-coordinate search methods are to succeed. For instance, when the 
objective function has only one local maximum (minimum), or if the local 
maxima (minima) all correspond to nearly equal objective-function 
values, then it will be relatively easy to obtain a near-optimal 
solution by either method. If, on the other hand, there exist many 
local maxima (minima) at which the objective-function value is much 
smaller (greater) than its global maximum (minimum) value, then there is 
a greater likelihood of arriving at a solution which is much poorer than 
the true optimum. 
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Therefore, the topography of the objective function used for the 
satellite orbital and frequency assignment synthesis, i.e., Equation 
(3.1), will be examined next. 

2. Relation Between Variations of Orbital /Frequency Separations 
and Single-Entry C/I c Value 

A very important concept should be pointed out first: for two 

satellites and their corresponding service areas, the single-entry C/I e 
values at all the test points increase as the satellite orbital/ 
frequency separations increase. 

In the C/I e calculation several terms are relatively invariant when 
satellite locations change. The received carrier powers are calculated 
from Equation (2.5), and the interference powers are calculated from 
Equations (2.6) and (2.13a). The satellite effective isotropic radiated 
powers (EIRP), defined as the product of the transmitting power and the 
transmitting antenna gain, are assumed the same for both satellites, 
regardless of orbital locations. Because the service areas are covered 
by the main beams within the -3 dB contour, the carrier transmitting 
discrimination factor at the test point is always larger than, and close 
to -3 dB. Because the geostationary orbit radius is 6.6 times the Earth 
radius, the propagation distances from the satellites to all the test 
points vary little when satellite locations are changed. 
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One factor that dominates the variations of the single-entry C/I e 
values is the receiving discrimination term ) in the interference 

power, and that variation depends chiefly on the satellite spacing. Any 
relative change of the two satellite locations will change the ground 
antenna off-axis angles toward the interfering satellite. Through the 
reference pattern in Figure 2.3, this angular change Induces a 
substantial change in the interference power. The further the two 
satellites are separated, the lower the discrimination factors, and thus 
the lower the interference. Another factor, the interference transmit 
discrimination term Dgf'lti), also varies when the satellite locations 
are changed; however, the change of Dg^) is much less than the change 
of D c |(t r j) because (seen from the geostationary satellite toward the 
Earth) changes much less than tri (seen from the ground upward to the 
sky) when the location of the interfering satellite is changed. The 
single-entry C/I e values increase when the corresponding satellite 
spacing increases until the discrimination factor n ( j('P r i) is in the far 
side-lobe of the reference pattern; then the C/I e values become almost 
constant as the spacing keeps on increasing. This will be called 
"quasi -monotonic" variation: C/I e values decrease (increase) 

continuously or remain constant when satellite spacing decreases 
(increases) continuously. 

For a given service-area pair, the single-entry C/I e value depends 
mainly on the magnitude of the satellite spacing, and slightly on the 
mean satellite orbital location. (It will be shown in Chapter IV that 
it varies only slightly for a large range of this mean satellite 
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location.) Therefore, for the purpose of the discussion in this 
section, an approximation is made by assuming that, for a given 
service-area pair, the single-entry C/I e values at all the test points 
depend only on the magnitude of the satellite spacing and are 
independent of the mean satellite location. 

Another factor that dominates the variation of the single-entry 
C/I e values is the frequency offset in the frequency assignments. As 
explained in Section II. F, the frequency filtering factor of Figure 2.6 
must be included in the effective interference calculation when the 
carrier frequencies of the wanted and unwanted signals are different. 

The larger the frequency offset, the lower the filtering factor, and the 
lower the interference. Therefore, the variations of the single-entry 
C/I e values also depend strongly on the separation of the frequency 
assignments: for given orbital assignments, the further the frequency 

assignments are separated, the higher the single-entry C/I e values. 

3. Topographic features of the objective function 

a) Objective-function topography of three-satellite example, 
orbital variables only 

For simplicity a three-satellite system will be used to demonstrate 
some key points. The three satellites are Sj, S?, S 3 , and these symbols 
will also be used to designate their orbital locations. A single 
channel will be considered so that no frequency variables are involved; 
in this case I and I e are identical. 


58 


The configuration space for this example is described here. First 
one draws the three orbital variable axes and indicates the feasible 
region, as in Figure 3 . 4 . The line segment HB measures the feasible 
range for Sj, HO for S2, HF for S3. In this figure the three satellites 
are assumed to have the same feasible orbital range, i.e., HB = HD = HF; 
therefore the feasible region is a cube. Any point in this cube 
represents an assignment scenario; its coordinates Si, S2, S3 are the 
satellite orbital locations of this assignment. For example, the point 
H represents the assignment in which all three satellites are collocated 
at one end of the orbital arc; the point A is for three satellites 
collocated at the other end; the point B is for Si located at one end 
while S2 and S3 are collocated at the other. Note that associated with 
every point in the configuration space there is an objective-function 
value which can be calculated by Equation ( 3 . 1 ). 

Several important features need to be mentioned. First, any point 
on the line AH indicates a three-satellite collocation. Then, all the 
points In the shaded plane ABHE correspond to S2, S3 collocation; those 
in the plane ACHF to Si, S2 collocation; those in the plane ADHG to Si, 
S3 collocation. Points within any one of the six sub-regions separated 
by the three collocation planes have the same satellite permutation. 

For example within the sub-region bounded by the planes AHB, AHC, ABC 
and HBC, all the points have the satellite permutation Si>S2>S3 ; this is 
indicated by the notation 1 - 2-3 in Figure 3 . 4 . Since the objective- 
function value is large when the C/I values are small, the objective 
function should have the highest values along the line AH because the 
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C/I values are the smallest for collocated satellites. Near the 
pairwise collocation planes ABHE, ACHF, ADHG the objective function 
should be moderately high, while between these planes the objective 
function falls off to form valleys because all the satellites are spread 
out. Thus we can visualize the objective-function topography as 
dominated by a system of ridges corresponding to satellite collocations; 
these ridges are connected to each other at the line AH where three 
satellites are collocated. Of course, the ridges will be high if the 
service areas are close together and hence the transmitting 
discrimination factor Dp in Equation (2.6) Is large; they will be small 
when the service areas are well separated. 

An arbitrary plane in this cube may be chosen to show the 
objective-function values corresponding to points on this plane. The 
plane chosen here is the plane CJKLFMNP shown in Figure 3.5(a). The 
objective function might have the shape shown in Figure 3.5(b). The 
base plane is divided into six sub-regions by its intersection with the 
shaded planes in Figure 3.4. Each intersection line represents the 
collocation of two satellites, and each sub-region represents one 
permutation of the three satellites. The objective-function topography 
sketch shows these sub-regions separated by the ridges representing 
two-satellite collocation. 

The objective-function values of a real set of scenarios will now 
be calculated here to confirm the above arguments. The three 
administrations are: Sj for Peru, S 2 for Bolivia, S 3 for Paraguay. The 

orbital locations chosen for the calculation constitute the shaded plane 


61 









RSTU shown in Figure 3.6; it is inside the region of permutation 
Si>S 2 >S 3 , and is a plane for which the value of S 2 is constant at 60°W. 
The value of Sj ranges from 60 to 64 degrees west, while the value of S 3 
ranges from 57 to 60 degrees west. The objective-function value is 
calculated from Equation (3.1), the values of C and I are calculated 
from Equations (2.5), ( 2 . 6 ) and (2.13a), the antenna reference patterns 
are taken from Figures 2.2 and 2.3, the ellipse data are calculated from 
[49]; the test points are those given in Table 3.1 for these 
administrations. The results are listed in Table 3.2, and the 
topography is plotted in Figure 3.7. As predicted, the maximum of the 
objective function in this plane occurs at R, the three-satellite 
collocation, the objective function rises to apparent ridges above the 
lines RS and RU which lie in the planes ABH and ADH of two-satellite 
collocations, and away from the lines of satellite collocation the 
objective function falls off and forms a valley. 

The typical topography of the objective-function surface for points 
lying on a plane inside any region of fixed permutation is shown in 
Figure 3.8. The point R' will be on the line AH, the lines R'S' and 
R'U' will be in the planes of two-satellite collocation, the points T' 
will be chosen in the R'S'U' plane. This figure illustrates the tyggcal 
shape and locations of the tip, ridges and valley. 
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Table 3.2 

List of Objective-function Values 


Z 


2 2exp[»-(C/I) H (dB)], a 

K J J 


5. 


S->(BOL) = 60 degrees. 


S : (PRU) 

s 3 (prg) 

z 


ln(z) 

60. 

60. 

0.2625 

E+4 

7.873 


59.5 

0.8918 

E+3 

6.793 


59. 

0.7644 

E+3 

6.639 


58. 

0.7606 

E+3 

6.634 


57. 

0.7601 

E+3 

6.633 

60.5 

60. 

0.1529 

E+4 

7.333 


59.5 

0.1084 

E+3 

4.685 


59. 

0.4488 

E+2 

3 .804 


58. 

0.4460 

E+2 

3 .798 


57. 

0.4455 

E+2 

3.797 

61. 

60. 

0.1010 

E+4 

6.918 


59.5 

0.5477 

E+2 

4.003 


59. 

0.27 81 

E-l 

-3.582 


58. 

0.1111 

E-l 

-4.500 


57. 

0.1104 

E-l 

-4.506 

62. 

60. 

0.9673 

E+3 

6.874 


59.5 

0.5322 

E+2 

3.974 


59. 

0.1490 

E-l 

-4.206 


58. 

0.1498 

E-4 

-11.109 


57. 

0.6218 

E-5 

-11.989 

63 . 

60. 

0.9613 

E+3 

6.868 


59.5 

0.5289 

E+2 

3 .968 


59. 

0.1454 

E-l 

-4 .23 1 


58. 

0.82 50 

E-5 

-11.705 


57. 

0.1852 

E-6 

-15 .502 

64. 

60. 

0.9596 

E+3 

6.867 


59.5 

0.5277 

E+2 

3.966 


59. 

0.1440 

E-l 

-4.241 


58. 

0.7938 

E-5 

-11.744 


57. 

0.1058 

E-6 

-16.062 
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Figure 3.8. Typical topography of objective function. 
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A point made in Section III.B .3 can also be illustrated by 
referring to Figure 3 . 7 . The gradient direction at the point x is 
almost perpendicular to the line RS. This means that when this gradient 
direction is followed, the orbital location of Sj is almost unchanged 
while S3 is placed further away from S2. Note that points close to the 
line RS correspond to nearly collocating the satellites S2 and S3; when 
these satellites are close to each other, the test points of their 
service areas have low C/I values. This means that the corresponding 
weighting factors in Equation ( 3 . 4 ) are high, and the gradient search 
process will push the two responsible satellites S2 and S3 apart. 

b) Objective-function topography of n-satelllte case 

The general shape of the objective-function topography for an 
n-satellite case will now be addressed. Note that the orbital and 
frequency variables are two different classes of decision variables so 
that their effects on the objective-function topography should be 
discussed separately. 

In the first step only the orbital variables are discussed; it is 
assumed that the frequency assignments are the same for all the 
satellites. For a case of n satellites to be assigned orbital 
locations, the n orbital variables constitute an n-dimensional 
configuration space, and it is divided into n! regions of different 
orbital permutations (i.e., satellite orderings) by n(n-l)/2 (n-1)- 
dimensional hyperplanes of two-satellite collocation. Each region of a 
fixed orbital permutation is surrounded by (n-1) hyperplanes of two- 
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satellite collocation and two (n-l)-dimensional boundary hyperplanes, 
each corresponding to one satellite at the boundary of its feasible 
orbital arc. The objective function has large values when satellites 
whose service areas are not well separated are collocated. For two 
neighbouring regions of different permutations, the corresponding 
permutations differ only in exchanging the relative positions of two 
satellites. Therefore the hyperplane between these two regions 
corresponds to the collocation of these two satellites, and the 
objective function may have large values in this hyperplane. All of 
these hyperplanes will connect at the line where all satellites are 
collocated; as a result the objective function has many local minima 
within the complete feasible region, and they can be characterized by 
their specific orbital permutations. 

In general a satellite is associated with an orbital location and a 
set of frequency allocations. Basically, the effect of the frequency 
variables on the topography depends on the frequency protection ratio 
pattern shown In Figure 2.6. There is one Important similarity between 
this pattern and the co-polarization antenna patterns: the relative 
protection ratio value is non-increasing as the frequency offset 
increases. So the spreading of the frequency assignments has the same 
effect on the objective-function value as the spreading of the orbital 
assignments: the further the frequency assignments are separated, the 

less the Interference, and the smaller the objective-function value. 
Hence the objective function may have large values when the frequency 
assignments are collocated (or nearly collocated since the plateau in 
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Figure 2.6 has a certain bandwidth), and local minima may occur where 
frequency assignments are spread out. 

The total number of local minima of the objective function in a 
case of n orbital and n frequency variables can be deduced as follows. 
Theoretically, the n orbital variables create n! regions with at least 
one local minimum in each region, then in each region the n frequency 
variables further create n! sub-regions with at least one local minimum 
in each sub-region. Therefore, there might be at least (n!) 2 local 
minima in this case. However, in reality the number of local minima is 
likely to be smaller, e.g., the collocation of two satellites may not 
result in bad C/I values to produce a hyperridge of large objective- 
function values when the service areas are well separated. 

Also, another important conclusion can be deduced from the above 
discussion: the local minima of the objective function may be 

characterized by the permutations of their orbital /frequency 
assignments. 
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c) Possibility of one local min 1mm for a fixed permutation of 
orbital /frequency assignments 

(1) Introduction 

Figures 3.6 and 3.7 also show that the objective function as a 
convex (or pseudo-convex) function with only one local minimum within a 
region of fixed orbital permutation. In this case this result should be 
obvious because, for a three-satellite example with only orbital 
variables, the objective function can be minimized by spreading the 
outside satellites as far as allowed by the feasible arc, i.e., until 
they reach the boundaries. It is interesting to speculate whether, for 
a general case of n satellites, there is only one local minimum in a 
region of fixed orbital /frequency permutation. If this conjecture could 
be shown to be true, or approximately true in the sense that all the 
minima for a fixed permutation have approximately equal objective- 
function values, then an ordinary (not extended) gradient search 
procedure would be sure to find an optimal, or at least near-optimal 
solution for a given orbital /frequency permutation. However, a 
definitive proof of the conjecture has eluded us. The remainder of this 
chapter presents evidence that it is likely to be true for the orbital 
variables. Therefore, in this section it is assumed that every 
satellite has been assigned the same frequency channel, and the 
discussion is confined to a region of fixed orbital permutation. 
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(2) Locations of Global Maximum and Local Minima of 
v Objective Function 

First, it is shown that the global maximum of the objective 
function corresponds to a scenario of all-satellite collocation by 
showing that the objective-function value decreases quasi -monotoni cal ly 
along any linear trajectory starting from al 1 -satel 1 ite collocation. 

Definition 1 

For two scenarios x, y in the n-dimensional configuration space, 
the linear trajectory between them is specified by the set of scenarios 
z such that 

z = ax + (l-a)y , (3.6) 

where a is a parameter with value 0 < a < 1. 

Lemma 1 

Assume that all the satellites have continuous feasible orbital 
arcs. Then given any two scenarios x, y, there exists a linear 
trajectory between them which is completely inside the feasible region; 
i.e., the feasible region is a convex set. 

Proof 

A scenario z in the n-dimensional configuration space is expressed 
by a l*n row matrix z={z\, Z2, ... z n ), where the components are the 
satellite locations of satellites 1, 2, ... n respectively. 

Along a linear trajectory of starting scenario x and final scenario 
y, a scenario z satisfies 
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z = ax + (l-a)y , 
where 0<a<l, or 

( ^ 1 » Z 2 , • •• Zp ) = a(xj» X 2 , ••• Xp) + 

(l-a)(yi, y2» ••• Yn) » 
or 

z-j = axi + (l-a)yi for l<i<n. 

Since all the satellites have continuous feasible orbital arcs, one has 
z-j eFi if xi , y-j eF^ for l<i<n, where Fi is the feasible arc for satellite 
i. Hence any scenario z along this linear trajectory is well defined, 
and the linear trajectory is completely inside the feasible region. So 
the feasible region is a convex set. 

Lemma 2 

Within a region of fixed orbital permutation in the n-dimensional 
configuration space all the pairwise satellite separations vary linearly 
along a linear trajectory. 

Proof 

From Lemma 1, for scenario z the pairwise satellite separation 
between satellites i, j is | z-j -z j | . From Equation (3.7c) one has 

(z.-Zj) = a(x i -x j ) + (l-aMy.-yj) , (3.8a) 

where 0<a<l. Since all the scenarios are of the same orbital 
permutation, (x-j-x-j) and (yi-yj) must have the same sign. Then, since 


(3.7a) 

(3.7b) 

(3.7c) 
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both a and (1-a) are positive, (z^-zj) will also have this same sign. 
Therefore, Equation (3.8a) can be rewritten as 

|zi-zj| = a | xi -xj | + (l-a)|yi-yj| , (3.8b) 

i.e., all the pairwise satellite separations vary linearly when a 
scenario is varied along a linear trajectory. 

Theorem 1 

The global maximum of the objective function is on the line of 
all -satellite collocation. 

Proof 

A linear trajectory that starts with all -satellite collocation is 
illustrated in Figure 3.9, note that the ordering of the satellites 
remains the same. Since the initial pairwise spacings are all zero, 
from Lemma 2 all the satellite spacings increase linearly in the 
trajectory. As discussed in Section III.C.2, all the single-entry C/I 
values increase quasi -monotoni cal ly when the corresponding pairwise 
satellite spacings increase linearly. As a result, all the aggregate 
C/I values, which involve the summation of all their contributing 
single-entry interference powers, must also increase quasi -monotoni caly 
in the process. Thus the objective-function value, which is the 
summation of the exponential of the negative of the aggregate C/I 
values, must decrease quasi -monotoni cal ly. Hence, the global maximum of 
the objective function is on the line of all -satellite collocation. 
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Figure 3.9. Linear expansion of orbital assignments. 
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Theorem 2 

Local minima of the objective function occur at boundaries of the 
feasible region such that all of the pairwise most widely separated 
satellites are at opposite limits of their respective feasible arcs. 

Proof 

Let the satellites be numbered according to their orbital 
positions, so that xi<x 2 <...x n . Let the lower limits of their feasible 
arcs be denoted e^ (1=1,2 n) and the upper limits uif (i=l,2,...,n). 

Assume xi*ei, i.e., satellite 1 is not at the lower limit of its 
feasible arc. Then 

z — a(ei,X 2 ,X 3 x n ) + (l-a)x , 0<a<l (3.9) 

defines a linear trajectory for which all the |zj-Zj| (j=2,3 n) must 

increase continuously away from the other satellites. All other 
spacings remain unchanged on the trajectory. As discussed in Section 
III.C.2, this increases the single-entry C/I values involving satellite 
1 quasi -monotoni cal ly while all other single-entry C/I values remain 
unchanged; therefore the objective-function value will decrease 
quasi -monotonical ly along the trajectory. Therefore it can not be a 
local minimum for the scenario x with xj+ej. This proves the theorem. 


76 



(3) Topography of the worst single-entry C/I value 

In the following section the topography of an auxiliary function A, 
defined as the the worst (smallest) single-entry C/I value, will be 
investigated. Although it is not the objective function Z, it is hoped 
that these two functions may have sufficient resemblance that some of 
the properties about Z can be extrapolated from A. First the function 
is defined as 

A = Min [(C/I). (3. 10) 

ieK, keK, jeJ k 1,KJ 

where the symbols mean 

(C/I)i,kj : the single-entry C/I value at test point j of 

service area k and the single-entry interference 
power is from satellite i, 

Jk : index set of test points of service area k, 

K : index set of satellites in a scenario. 

The value of A for a given scenario Is obtained as follows. First 
all the single-entry Interference from any one satellite to all the test 
points of other satellite systems are calculated using Equation (2.6), 
and these values are denoted as Ii } kj when i is the interfering 
satellite and j is the test point in service area k. The carrier power 
at any test point is calculated using Equation (2.5), and is denoted C|<j 
when j is the test point in service area k. The value of (C/IJ-j^j is 
calculated by taking the ratio of C^j and Ii ,k j * The value of A is 
obtained from choosing the smallest of all the (C/O^kj terms for that 
scenario (note that a scenario is a point in the configuration space). 
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Note that, if this function is used as the objective function for the 
extended-gradient or cyclic-coordinate search method, then the most 
favorable assignment is that for which the objective function is the 
global maximum. 

Considering only the orbital variables, we wish to show that this 
function has at most one local maximum (hence a global maximum when it 
exists) for each ordering (permutation) of the satellites, and that such 
a maximum must lie on a boundary of the feasible region. Some lemmas 
and definitions will be useful in the proof of this theorem. 

Lemma 3 

Within a region of fixed permutation in the feasible region, all 
the (C/I)jj|< terms are quasi -monotonic functions along any linear 
trajectory (orbital variables only). 

Proof 

As discussed in Section III.C.2, for a service-area pair, the 
single-entry C/I values at all the test points vary quasi -monotoni cal ly 
as the satellite spacing varies. According to Lemma 2 all the pairwise 
satellite spacings either increase or decrease linearly, or remain 
constant along a linear trajectory within a region of fixed permutation. 
Therefore all the single-entry (C/I )i ,kj functions along that trajectory 
vary quasi -monotoni cal ly. 

Definition 2 

For two curves that have a countable number of common points, a 
vertex of the two curves is a common point of the two curves. 
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Let a function Y(x) be defined from a set of quasi -monotonic 


functions Y-j(x), iel in a finite range as 

Y(x) = Min [Yi(x)] . 
iel 


(3.11) 


Lemma 4 

If the local maximum of Y(x) exists, it must occur at a vertex or 
at the boundary of the allowed range. 

Proof 

Assume that Y(x 0 ) is a local maximum (minimum) of Y(x), and it is 
not a vertex of Y-j(x), iel, nor is it at the boundary. Then Y(x) is 
equal to one of the Y^(x) in the vicinity of x 0 ; hence Y(x) is monotonic 
in the vicinity of x 0 . Being a local maximum (minimum) means that 
Y(x)<Y(x 0 ) (Y(x)>Y(x 0 )) for x in the vicinity of x 0 , which contradicts 
the fact that Y(x) is monotonic in the vicinity of x 0 . Hence the local 
maxima (minima) of Y(x) occur at the vertices or at the boundary. 

Lemma 5 

Y(x) can not have more than one local maximum. 


Proof 

Suppose that, as shown in Figure 3.10, the function Y(x) has two 
separated local maxima at xj and X3, then there must exist a local 
minimum at some point X£ between xj and X3. By lemma 4, this local 
minimum must be at the vertex of two Yi(x), say Yj(x) and Y 2 ( x ) . The 
point X 2 being a local minimum of Y(x) requires, in the vicinity of X 2 , 
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X 2 X 3 


X 


1 


Figure 3.10. Hypothetical case where two maxima occur in a linear 
trajectory. (Proved impossible by contradiction) 
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Y(x) > Y(x 2 ) 


for x > and x < x^ 


(3.12) 


Let Yi(x) be the function which coincides with Y(x) for x<X£ (i.e., the 
curve 1i) and Y 2 (x) the function which coincides with Y(x) for x>X 2 
(i.e., the curve I 2 ) in the vicinity of X 2 » then one has 

Yj(x) > Y(X 2 ) for x < X 2 . (3.13) 

Yj(x) being quasi -monotonic means that, from Equation (3.13), 

Yi(x) < Y(X 2 ) for x > X 2 (3.14) 

in the vicinity of X 2 . The same reasoning that led to these two 
equations, when applied to Y 2 (x ) , gives in the vicinity of X 2 , 

Y 2 (x) > Y(X 2 ) for x > X 2 , (3. IB) 

Y 2 (x) < Y(X 2 ) for x < X 2 . (3.16) 

Remembering that Y(x) coincides with Y 2 (x) (by definition of Y 2 (x) 
above) for x>X 2 , Equation (3.15) can be rewritten 

Y (x) > Y(x 2 ) for x > X 2 , (3.17) 

and together with Equation (3.14) this gives 

Yj(x) < Y(x) for x > X 2 (3.18) 

in the vicinity of X 2 . This contradicts the definition of Y(x) in 
Equation (3.11). Thus Y(x) can have at most one local maximum. 
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Theorem 3 

Along any linear trajectory within a region of fixed permutation in 
the feasible region, the auxiliary function A has at most one local 
maximum (orbital variables only). 

Proof 

The discussion here is confined to a linear trajectory within a 
region of fixed permutation. From Lemma 3, all the single-entry C/I 
terms are quasi -monotonic functions on such a trajectory. The function 
A is defined as the minimum of all these C/I terms; thus the function A 
may be identified with Y(x) in Lemma 5 if the trajectory parameter a is 
identified with x in that lemma. Therefore, from Lemma 5, the function 
A has at most one local maximum (hence a global maximum when it exists) 
along a linear trajectory. 

Theorem 4 

The auxiliary function A has at most one local maximum within a 
region of fixed permutation in the configuration space (orbital 
variables only). 

Proof 

The discussion here is still confined to a region of fixed 
permutation. It will be shown that Theorem 3 is contradicted if the 
function A has two local maxima within such a region. Suppose that the 
function A has two local maxima at x^ and X 2 » then the linear trajectory 
that passes through both xj and X 2 will have two local maxima at and 
X 2 ; however, this is a clear contradiction of Theorem 3. Hence within 
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this region the function A has at most one local maximum; hence any 
maximum must be global in the region. Because the C/I variations are 
only quasi -monotonic, there may be sub-regions in which the value of A 
is constant; if this value is not exceeded elsewhere in this region, A 
will not have a distinct maximum in this region. 

Theorem 5 

The local maxima of the auxiliary function A are located at the 
boundary of the feasible region in the configuration space (orbital 
variables only). 

Proof 

By the same argument as in Theorem 2, a scenario located in the 
interior of a feasible region cannot be a local maximum because it is 
always possible to find a scenario located at the boundary that has 
better or equal single-entry C/I values. Thus any maximum of A has to 
be at the boundary. 

(4) Numerical test 

As discussed in Section III.B.3, the objective function Z is 
formulated to emphasize the worst aggregate C/I e terms by exponentiating 
the negative of the aggregate C/I e value; thus it is reasonable to say 
that there is a direct relationship between the function Z and the 
negative of the worst aggregate C/I e term, or several such terms if they 
are of approximately equal value. In view of Theorem 4, this suggests 
that the function Z is likely to have only one local minimum within a 
region of fixed permutation. Since a rigorous proof has eluded us, a 


83 



numerical example of a four-administration model was calculated to 
investigate x this feature: the four administrations are Argentina, 

Bolivia, Paraguay, and Peru; the feasible orbital range is from 62 to 68 
degrees west. The objective-function value was calculated from Equation 
(3.1), the values of C and I were calculated from Equations (2.5), 

(2.6), and (2.13a); the antenna reference patterns are from Figures 2.2 
and 2.3; the ellipse data were calculated from [49], Since it is known 
from Theorem 2 that the local minima must be located at the boundaries 
of the feasible region in the configuration space, the objective 
function values at all the 24 boundary planes of this numerical example 
were calculated; the equal -height contour plots are shown in Appendix C, 
where Argentina, Bolivia, Paraguay and Peru are denoted as ARG, BOL, 

PRG, and PRU respectively. In each figure, two satellites are located 
at opposite ends of the feasible arc, and the (varying) orbital 
locations of the other two are indicated by the coordinates. The upper 
triangle represents one permutation and the lower triangle another 
permutation. Clearly there is only one local minimum for each 
permutation. The results suggests that for a general case (n 
satellites) it is likely that there is only one local minimum within a 
region of fixed orbital permutation. 

D. DISCUSSION AND CONCLUSION 

Even though the objective function in Equation (3.1) is not the 
only one that might be formulated to solve the orbital /frequency 
assignment problem, it still should be representative of objective 
functions designed to maximize low C/I e ratios. Therefore many of the 
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mathematical features of the assignment problem can be inferred from 
this function. 

Referring to Figures 3.4 and 3.5 and the discussions in Section 
III .C , the objective function of Equation (3.1) often has large values 
at points of either orbital or frequency assignment collocation, and the 
local minima can be characterized by the permutations of their orbital/ 
frequency assignments. This feature indicates that there might be on 
the order of (n!) 2 local optimum solutions for a n-satellite case of n 
orbital variables and n frequency variables; in reality the number of 
local optima is likely to be smaller because the collocation of two 
satellites may not result in bad C/I e ratios when the service areas or 
their frequencies are well separated. Still, one thing is clear: 
permutation of the orbital /frequency assignments is an important part of 
the problem. Therefore a technique for finding the globally optimal (or 
a near-optimal) scenario must be able to choose the proper permutation 
of the orbital /frequency assignments, as an important step toward 
satisfying the C/I e requirement. 

Also, the discussion in Section III.C.3.C, which shows that the 
objective function is likely to have only one local minimum within a 
given orbital permutation, suggests that in order for an extended 
gradient search method to obtain the globally optimal solution, it 
should terminate as an ordinary gradient search with the optimal 
permutation in orbital /frequency assignments. 
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CHAPTER IV 


OPTIMAL ORBITAL ASSIGNMENTS BY MEANS OF THE AS CONCEPT 
A. INTRODUCTION 


In this chapter a direct correspondence between the single-entry 
C/I e protection requirements and the necessary satellite spacings is 
exhibited. This relationship is used to formulate linear constraints to 
enforce single-entry C/I e protection requirements between all pairs of 
satellites. Two new formulations are then developed under the 
assumption that frequency assignments are the same for all the service 
areas. One is a mixed integer linear program, solved by a branch-and- 
bound procedure. The other is a linear program with both linear and 
nonlinear side constraints, the simplex method with restricted basis 
entry can be used to find an approximate solution when this formulation 
is used. As a consequence, the cumbersome nonlinear C/I e expression 
used for the synthesis formulation in Chapter III is avoided. 

The 1977 World Administrative Radio Conference (WARC-77) suggested 
that for maximum orbit utilization, space stations should be placed as 
close to each other as is consistent with keeping the mutual 
interference to acceptable levels [16]. This concept was explored by a 
Canadian study group by relating the single-entry C/I protection 
requirement to the satellite spacing [26,27,28]; they formulated a 
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orbital /frequency assignment program, for use at the 1983 Regional 
Administrative Radio Conference (RARC-83), based on the satellite 
spacing requirement. The concept of using this relationship in 
conjunction with linear optimization is original, to our best knowledge. 
The author is indebted to Professor Charles H. Reilly of the Department 
of Industrial and Systems Engineering for the mixed-integer and linear 
programming formulations and to Mr. David J. Gonsalvez for the actual 
programs. 

The organization of the chapter is as follows. First the 
relationship between the single-entry C/I e protection requirement and 
pairwise satellite spacing will be formulated, and the relative 
importance of system parameters will be discussed* Next, exact and 
approximate methods of calculating the required spacing will be 
presented. This will be followed by a heuristic discussion of the 
relationship between single-entry and total acceptable protection 
requirements, in order to establish single-entry requirements which are 
highly likely to lead to satisfaction of the total acceptable C/I e 
protection requirement. Next the principle of the methods will be 
elucidated with a very simple hypothetical four-service area example to 
show how satellite ordering (permutation) enters into the linear 
optimization process. Finally the results of both the mixed integer 
program (MIP) and restricted-basis entry linear program (RBLP) 
formulations will be presented for a scenario of six South American 
administrations. 


c-x. 
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B. RELATION BETWEEN SINGLE-ENTRY C/I e PROTECTION REQUIREMENT AND 
REQUIRED SATELLITE SPACING 


The relationship between a single-entry C/I e requirement and the 
required satellite spacing will now be derived. The general 
configuration of two service areas A, B and the locations of their 
satellites Sr is shown in Figure 4.1 which is identical to Figure 
2.7. Refer to Section II. F for detailed description of all the 
parameters. It is assumed that frequency assignments are the same for 
both satellites. The received carrier power, C, in channel n (carrier 
frequency f n ) at test point d is given in Equation (2.5) as 


C = 


P A* G A* D A^tc ),G d* c< 

f 2 * (4 it) 2 -X 2 
n 


(4.1) 


the effective single-entry interference power, I n , from Sg to channel n 
of test point d can be inferred from Equation (2.6) as 


i e = z 
" m 


VW^-VV^ri^ 

f m 2 * (4ir) 2 *y 2 

VVV t ti K v D Ai ) - c2 F <W 

. £ 


(4ir ) 2 *y 2 


m f 


2 ’ 


(4.2) 


m 


where the summation is over all frequency channels assigned to S B . The 
term "single-entry" is defined as the aggregate of emissions from any 
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one satellite entering any receiver in the wanted service within the 

6 

channel to fee protected [87]. The single-entry I n values are different 
for different channels, and the worst one is in the channel at the 
center of the assigned band; it is denoted as I e , and is used to 
evaluate the Interference effect in this chapter. 

The exact single-entry C/I e ratio (in center channel n) at test 
point d is therefore 


C P aV x? V^tc) 

l e ~ P B* G B /y2 ’ 


f Z *[z(f ,f )/f 2 ] . 

n L m ' n* ir m J 


(4.3) 


Even though the only independent variables in Equation (4.3) are 
the satellite orbital locations, there are many hidden parameters and 
relationships. The minimum ellipse (its size, orientation, and aim 
point) for is a function of satellite location. It is assumed in 
this chapter that the carrier power flux densities at the aim points are 
equal for all the satellite systems. For any ellipse size there is a 
corresponding G/\ value, and the value of P^ must be adjusted to give the 
required power density at the aim point. As for DA('Ptc)> the value of 
itc depends on the satellite location, and there is also an implicit 
parameter <ltco which is a function of the ellipse (refer to Figure 2.1). 
The same considerations apply to satellite Sg, and Gp, Pg, Db. Strictly 
speaking, the values of G<j, and thus D<j, are different for different 
channels because they are functions of carrier frequencies; in this 
chapter the values of Gj for all channels are assumed the same, which is 
reasonable for a narrow band frequency assignment. 
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The exact calculation of Equation (4.3) is complicated, it is 
therefore worthwhile to seek an approximation to Equation (4.3). 

Because the geostationary orbit radius is 6.6 times the Earth radius, 
the propagation distances x and y are approximately equal. Then, since 
the power densities at the respective aim points are assumed equal [65], 
it follows that 

P A* G A P B‘ G B 

)(2 a y2 • (4.4) 

When the service-area shape can be fitted reasonably with an ellipse, 
the test point d should be located on or near the -3 dB contour [46,47], 
giving 


D A<*tc> ’ V2 ' 


Equation (4.3) can therefore be approximated as 


Cl 1 

~e " 2.D B (+ t1 )-D d (* p1 ) * (VW 


(4.5) 


(4.6) 


The factor l/{ f n 2 * [EF(f n ,f m )/f m 2 ]} is a constant since it is 

m 

assumed that the frequency assignments are the same for all satellites 
(note that f n is at the center channel). An internationally agreed 
F(f n ,f m ) reference function exists only for the BSS, but not for the FSS 
which is the subject of the study in this chapter. However, it is known 
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that numerous U.S. FSS satellites have 24 transponders in the 6/4 GHz 
band occupying the whole 500 MHz bandwidth, each transponder channel is 
36 MHz wide and the guard band between two channels is 4 MHz (cross- 
polarization discrimination allows frequency re-use in a satellite). In 
this chapter only the co-channel interference will be considered. Thus 

the value of l/{ f n 2 * [EF(f n ,f m )/f m 2 J is taken to be one. Then Equation 

m 

(4.6) can be approximated as 

C 1 

l e “ 2 * D B^ti)* D d^rl) ’ (4 * 7) 

The variation of the C/I e values with respect to the satellite 
orbital locations will now be discussed. Because the service areas are 
stationary and the satellite orbit radius is 6.6 times the earth radius, 
the value of tti » and hence the value of Db ('P ti ) » changes little when 
the location of Sg, the interfering satellite, is changed by a small arc 
length. When the service-area pair and system parameters are given, the 
term 4>ri » and thus n,j(^ r -j), becomes the only factor that can 
significantly affect the C/I e values; this is done by changing the 
satellite spacing to vary \|> r i , and thus the Dd(^ r i ) value. 


92 



To have an acceptable single-entry C/I g ratio, the main 

contributions are seen to come from the terms Dg ) and ^('J'p.j)* The 

term D d ( vp* - ) comes from the separation between the two service areas; 

the further they are separated, the less the value D R (t t ^). The very 

fact that it comes from the geographic separation of the service areas 

makes it a valuable resource for the system planner. The term 0^(4^) 

comes from the satellite spacing; the larger the spacing, the larger the 

value ., and the smaller the value n.(\p .). The most important 
n an 

feature about the satellite spacing is that this quantity can be 

controlled by the system planner. 

For a given single-entry C/I fi protection requirement, the system 

planner should first look for the term D B (<|^) for its contribution to 

the margin between C and I g ; when this is not enough, he then has to 

look for the term D,(^ .) to make up the difference. For satellite S. 

d'rv a 

InratoH at 1 And lnratart tn the Met nf the threshnlH satellite 

. B - "A’ 

spacing that lets the resulting worst single-entry C/I g value equal the 
protection requirement is the local required spacing for the satellite 
pair, and is designated as B(l+)’ tlere s ^ n + means east 

1. For Sg to the west of it would be AS A pyg^i Wl 'th any less 
spacing at least one of the test points would have an unacceptable 
single-entry C/I g ratio, while with any more spacing all the single- 
entry C/I g ratios would be better than the requirement. Thus, the As(l) 
function can be viewed as the reflection of the single-entry protection 
requirement, eg., | S A -S g |> as ab would guarantee satisfactory single-entry 
C/I e ratios at all the test points in service areas A and B. 
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Since a satellite spacing requirement given by as ( 1 ) is equivalent 
to the equivalent single-entry protection requirement, it can be taken 
as a constraint on the relative locations of the corresponding 
satellites. So, instead of thinking about a scenario having 
satisfactory single-entry C/I e values, one may think about the scenario 
satisfying all the as(1) constraints. One important feature about the 
As (1 ) value is that, in producing the margin between C and I e , it fully 
utilizes the service-area separation to minimize the necessary satellite 
spacing. This approach corresponds directly to the idea proposed in 
WARC-77: for maximum orbit utilization satellites should be placed as 

close to each other as is consistent with keeping the mutual 
interference to acceptable levels [16]; hence the set of As(l) values is 
exactly what the system planner should use in order to achieve this 
goal. Using the As(l) value in the orbit planning task has a great 
advantage in terms of numerical calculations. While the C/I e expression 
involves many geometric equations and is highly nonlinear with respect 
to the orbital variables, the as ( 1 ) set can be calculated once and for 
all for each service-area pair and then used as constraints. This 
greatly changes the aspect of the orbital assignment task, and the 
methods to solve it; this will be apparent in Section IV. F. 

It is suggested in WARC-77 that single-entry protection 
requirements can be used as a guide for determining sharing criteria 
[87]; still the total interference from all sources must be calculated 
to evaluate the scenario definitively. In this chapter, satellite orbit 
planning methods are developed based on single-entry protection 
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requirements; it will be shown in Section IV. D that this is likely to 
lead to adequate aggregate protection when the single-entry C/I e 
protection requirement is larger than the total acceptable protection 
requirement by 5 dB. First, one needs to turn one's attention to 
algorithms for calculating the As values just defined. 

C. CALCULATION OF AS VALUES 

1. Exact Method 

An exact determination of the required satellite spacings requires 
solution of Equation (4.3) for i|^ r -j for various locations of S^. An 
explicit solution has not been found; still, a numerical value can be 
obtained by evaluating the right side of Equation (4.3) (or equivalently 
Equations (4.1) and (4.2)) for increasing separations until the required 
C/Ia value results. Specifically, the rigorous calculation of the 
threshold satellite spacing for two service areas A, B may be done by 
the following algorithm. (Refer to Figure 4.2 for the geometric 
relations.) 

(1) Set the location of satellite S/\ at orbital location 1. 

(2) Move the location of satellite S B incrementally from 1 toward 
the east. After each move, use the streamlined Spectrum/Orbit 
Utilization Program (SOUP) code in Appendix A to calculate the 
C/I e values at all the test points of the two service areas. 
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Figure 4.2. Configuration of as(1) value calculation. 



(3) There exists a satellite spacing beyond which all the C/I 

e 

values exceed the single-entry protection requirement, and 
below which at least one C/I g value is worse than the 
protection requirement. This is the satellite spacing for the 
two service areas at the prescribed locations and is denoted as 
as A(1) r( 1 +) ’ where S1 9 n + means east of 1. 

(4) Repeat procedures (l)-(3) for S D west of S.; the resulting 

spacing is denoted as as a ^ ^ where the sign - means west 

of 1. 


(5) Repeat ( 1 ) -( 4) with a set of new locations for satellite S 
until its feasible arc has been covered. 

The above calculation is rather time consuming. For each test 
position, the minimum ellipse data of the service area has to be 
generated. 


A 


2. Approximate Method 

To ease the computational burden, an approximate method was 
adopted. This method is based on the assumption that the power 
densities at the ground receiver locations do not change when the 
location of either satellite is moved away from 1 by a small arc length. 
The procedure is: 

(1) Collocate the two satellites at 1. Use the streamlined SOUP 
code in Appendix A to calculate all the ground receiver C/I e 
values. 
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(2) Pick out the test point with the worst C/I e value, calculate 
the. margin between this C/I e value and the single-entry 
protection requirement. At this point the service-area 
separation factor has been accounted for, so this margin should 
be made up by separating the two satellites from collocation to 
produce the receiving discrimination loss Dj (^ r i ) toward the 
interfering satellite. 

(3) Use the receiving antenne reference pattern in Figure 2.5 to 
calculate the necessary off-axis angle that provides this 
margin; this is the topocentric angular separation (viewed from 
the test point with the worst C/I e value) the two satellites 
should have. 

(4) Spread the satellites symmetrically apart from 1 step by step 
until the separation gives this necessary topocentric angle at 
this test point. The final geocentric separation is the 
approximate as(1) value when S ^ and Sr are in the vicinity of 
orbital location 1. 

The computer code for this calculation is listed in Appendix D. 

An example of the As (1 ) calculation is shown in Table 4.1 with the 
two administrations being Bolivia (BOL) and Paraguay (PRG). The 
satellite transmitting and ground receiving reference patterns are given 
in Figures 2.4 and 2.5. Initially the satellites are located at 90 
degrees west. The interference calculation is carried out, and the test 
point with the worst C/I e value, -0.66 dB, is at 62.2 degrees west 
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Table 4.1 


Example of as(1) Calculation Procedure 


COUNTRY 

SATELLITE 

(LON. ) 


FREQUENCY ■ 

( MHz ) 


EOL 

-90. 

.00 


4000.00 



PRG 

-90. 

.00 


4000 . 00 



TEST COUNTRY : EOL 


SATELLITE : 

-90.00 



TEST 

POINT 


I NT. SAT. 

C/I (dB) 

MARGIN 

(dB ) 

LON. 

LAT. 







-65.00 

-12.20 


PRG 


28.86 

-1.14 


-65.50 

-9.80 


PRG 


27.51 

-2.49 


-69.00 

-11.20 


PRG 


27.78 

-2.22 


-60.00 

-16.10 


PRG 


B .02 

-21.98 


-57.50 

-18.00 


PRG 


2.58 

-27 .42 


-67.50 

-22.70 


PRG 


14.88 

-15. 12 


TEST COUNTRY : PRG 


SATELLI 

TE : 

-90.00 



TEST 

POINT 


INT. SAT. 

C/I (dB) 

MARGIN 

(dB) 

U) 

r- 

LO 

1 

-25.30 


EOL 


7.75 

-22.25 


-58.60 

-27.30 


BOL 


10.58 

-19.42 


-56.20 

-27.20 


BOL 


14.54 

-15.46 


-54.70 

-25.50 


BOL 


13.22 

-16.78 


-54.20 

-24.10 


BOL 


1 1 .20 

-18. G0 


-58. 10 

-20.20 


BOL 


1.15 

-28.35 


-59. 10 

-19.30 


BOL 


0.08 

-29 . 92 


-62.20 

-20.50 


BOL 


-0.66 

-30.66 


-62.70 

-22 . 20 


EOL 


-0.25 

- 30 . 2 5 


-SS . 70 

-27.20 


EOL 


10. 15 

-1S.C5 


WORST MARGIN 

IS -30.1 

66 

dB AT 

PRG 

; ( -62.20, 

-20.50) 


REQUIRED SATELLITE SPACING : 4 

.00 

AT -90.00 

FOR C/I 

30.0 


99 



longitude, 20.5 degrees south latitude in Paraguay. The single-entry 
C/I e protection requirement is chosen to be 30 dB, hence it requires 
30.66 dB attenuation from the ground receiving discrimination to provide 
the necessary margin between C and I e . Then the receiving reference 
pattern is used to calculate the required off-axis angle, i.e., the 
topocentric angle, of the two satellites as seen from the test point 
with the worst C/I e value. The satellite geocentric spacing is obtained 
from step (4) of the approximate method, the result is shown in Table 
4.1 which shows that, in the vicinity of 90 degrees west, the necessary 
satellite spacing for these two service areas is four degrees. This 
spacing should result in at least 30 dB single-entry C/I e values at all 
test points, with the worst one at 30 dB. This result can not be 
guaranteed precisely, since the test point which was the worst for the 
original satellite locations may not be the worst for the new locations, 
and since new ellipses were not generated for the new satellite 
locations. 

Using the same antenna parameters and 30 dB protection requirement, 
some of the as(1) values for several service-area pairs at different 
mean satellite locations are listed in Table 4.2. The calculation is 
made in 10- or 20-degree increments. 

An example of the curve As(l) vs. mean satellite location 1 is 
shown in Figure 4.3, the two administrations are Paraguay and Uruguay. 
The As (1 ) value is fairly constant within a large range of satellite 
locations, and starts to increase when the satellite elevation angle is 
small . 
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Table 4.2 

As(l) Values of Six South American Administrations 


country pair 

longitude 

AS 

country pair 

longitude 

AS 

ARG BOL 

-70 

4.00 

CHL PRG 

-20 

3.05 


-80 

4.02 


-30 

2.45 


-90 

4 . 05 


-40 

1 .69 


-100 

4.12 


-50 

1.35 


-110 

4.17 


-60 

1.14 





-70 

l .03 

ARG CHL 

-70 

4.18 


-80 

1.14 


-80 

4.05 


-90 

1.25 


-90 

4.00 


-100 

1 .46 


-100 

4.02 


-110 

2.00 


-110 

4.19 


-120 

3.32 

ARG PRG 

-70 

4.24 

CHL PRU 

-70 

3.84 


-80 

4.28 


-80 

3.83 


-90 

4.32 


-90 

3.85 


-100 

4.28 


-100 

3.83 


-110 

4.32 


-110 

3.94 

ARG PRU 

-70 

0.94 

CHL URG 

-20 

2.52 


-80 

1.04 


-30 

1.44 


-90 

1.15 


-40 

0.93 


-100 

1 .25 


-50 

0.80 


-110 

1.41 


-60 

0.43 





-70 

0.42 

ARG URG 

-70 

4.18 


-80 

0.41 


-80 

4.14 


-90 

1.03 


-90 

4 06 


-100 

1.28 


-100 

4.06 


-110 

1.59 


” 1 1 0 

3 . 94 


- 1 20 

2.10 

BOL CHL 

-70 

4.13 

r RG r R Li 

-20 

0. 45 


-80 

4.20 


-40 

0.47 


-90 

4.28 


-60 

0.48 


-100 

4.39 


-70 

0. 49 


-110 

4.57 


-80 

0.49 





-90 

0-50 

BOL PRG 

-70 

4.00 


-100 

0.85 


-80 

3.99 


-110 

1 . 10 


-90 

4.00 


-120 

1 .76 


-100 

4.03 





-110 

4.04 

PRG URG 

-20 

2.35 





-40 

2.20 

BOL PRU 

-70 

3.87 


-60 

2.13 


-80 

3.95 


-70 

2.16 


-90 

3.99 


-80 

2.19 


-100 

4.10 


-90 

2.20 


-110 

4.26 


-100 

2.34 





-110 

2.46 

BOL URG 

-20 

0.43 


-120 

2.64 


-40 

0.42 





-60 

0.40 

PRU URG 

-20 

0.45 


-70 

0.39 


-40 

0.43 


-80 

0.38 


-60 

0.41 


-90 

0.38 


-70 

0.40 


-100 

0.84 


-80 

0.37 


-110 

0.94 


-90 

0.37 


-120 

1 .09 


-100 

0.33 


-110 0.32 

-120 0.36 
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Although the as ( 1 ) values are calculated up to two decimal 
fractions, in practical situations only the first decimal fraction is 
meaningful because the typical satellite station keeping inaccuracy is 
n.l degree [44]. The reason for keeping two decimal fractions is only 
to show their variations more precisely. 


3. Validity of Approximate Method 


The approximate method gives the correct result if the power 
densities at all the test points do not change when the satellites are 
moved by a small arc length. Strictly speaking, the power-density 
invariance assumption is not correct because the carrier power densities 
are designed to be constant only at the beam centers [65], However, in 
practise, it is a highly acceptable assumption. To demonstrate this 
point more clearly. Equation (4.3) is written as 

C/I e = K(1 A’ 1 B )/D d ( Vi ) * (4.8) 


where 

P A* V* 2 

K(1 A’V = P B * Gg/y 2 


D A^tc) 1 

¥*ti> ’ f n 2 *^( f n’ f m)/ f m 2 l 


(4.9) 


and 1 and lg denote the orbital locations of S/\ and Sf}. The 
approximate method consists of calculating K for collocation, i.e., 
K(l,l), then assuming that this value remains constant as satellite 
locations are shifted a few degrees from 1. The discussion in Section 
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IV. B (specifically, the paragraphs covering Equations (4.4) to (4.7) and 
the one immediately following Equation (4.7)) does indicate that this is 
a very good approximation. 

Some examples support these arguments and show that the C/I e values 
at the ground receivers do not change by more than 1 dB when the 
satellites are moved by two degrees; an example will now be given to 
demonstrate this point. In this example the Bolivia and Paraguay 
satellites are located at 92 and 88 degrees west, respectively, a 
four-degree separation as suggested from the result in Table 4.1. The 
C/I e results from the streamlined SOUP calculation are given in Table 
4.3. Note that the worst C/I e value is very close to the 30 dB 
single-entry protection requirement, which indicates that the as(1) 
value of the approximate method is very close to that of the exact 
method. 

4. Relation Between Service-Area Adjacency and as(1) Value 

For any adjacent service-area pairs the As(l) values for a given 
protection requirement are approximately independent of the sizes and 
shapes of the service areas. This can be seen best from Equation (4.7); 
for adjacent service areas A and B, the test point d in service area A 
with the worst C/I e is on or near the common border with B, and its 
^B(^ti) be approximately 1/2 (-3dB). This determines the required 
^(♦ri)* and ♦ri* and hence as (1 ) . As another way of looking at it, the 
carrier and interference power densities along the common border are 
nearly equal because they are designed to be approximately 3 dB below 
the respective beam-center power densities, which are designed to be 


104 



Table 4.3 

C/I e Results to Show Validity of Approximate Method 


COUNTRY 

SATELLITE (LON.) 

FREQUENCY 

( YHz ) 

COL 

- 

92 .00 

4000.00 


PR.G 

- 

88.00 

A 000. 00 


TEST COUNTRY : BOL 

SATELLITE 

: -32.00 


TEST 

LON. 

POINT 

LAT. 

I NT. SAT. 

C/1 (dB) 

KARGIN 

-65.00 

-12.20 

PRG 

59.67 

29.67 

-65.50 

-3.80 

PRG 

56.32 

28.32 

-63.00 

-11.20 

PRG 

58.71 

28.71 

-60.00 

-16.10 

PRG 

38 . 82 

8 . 82 

-57.50 

-18.00 

PRG 

33 . 23 

3.23 

-67.50 

-22.70 

PRG 

46.03 

16.03 

TEST CO'JNTR 

Y ; PRG 

SATELLITE 

: -68.00 


TEST 

LON. 

POINT 

LAT. 

INT. SAT. 

0 / T ( rl R } 

MAP C T N 

-57.60 

-25.30 

BOL 

38.33 

S . 33 

-58.60 

-27.30 

BOL 

41.27 

11.27 

-56.20 

-27.20 

BOL 

44.88 

14.88 

-54.70 

-25.50 

BOL 

43.33 

13.33 

-54.20 

-24 . 10 

BOL 

41.26 

11.26 

-58.10 

-20.20 v 

BOL 

31.79 

1 .79 

-53.10 

-19.30 

BOL 

30.75 

0.75 

-22.20 

-20.50 

BOL 

30.06 

0.06 

-62.70 

-22.20 

BOL 

30.39 

0.39 

-58.70 

-27.20 

BOL 

40.86 

10.86 
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equal. Therefore, the required margin between C and I e Is provided only 
by the satellite spacing. Since this is true for all adjacent service- 
area pairs, the required satellite spacings are nearly the same. For 
the receiving reference pattern in Figure 2.5 and a 30 dB single-entry 
C/I e protection requirement, these As(l) values are approximately four 
degrees, as can be seen from several cases in Table 4.2. 

The As (1 ) values for two non-adjacent service areas are likely to 
be smaller, and depend on their shapes, sizes and the separation between 
them. This is evident from Equation (4.7) since 0g(<pt-j) is likely to be 
smaller numerically (also in dB) in this case compared to the adjacent 
case. Heuristically, because of the service-area separation, the 
interference power densities in these two service areas are likely to 
have a deeper transmitting antenna discrimination loss relative to the 
carrier power densities. Therefore, the system needs less receiving 
antenna discrimination loss, Dd(<|>n)> to achieve the required C/I e 
ratio. This means that less satellite spacing, or a smaller as(1) 
value, is needed. This is shown in several cases in Table 4.2. 

D. RELATION BETWEEN SINGLE-ENTRY AND TOTAL ACCEPTABLE C/I p PROTECTION 
REQUIREMENTS 

In satellite communications, a scenario is evaluated by computing 
an equivalent margin (see Equations (2.8) and (2.9)), which takes into 
account all the interference, at all the test points. In Chapter II 
this was shown to be equivalent to comparing the aggregate C/I e values 
with the total acceptable protection requirement. It is therefore 
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essential that an assignment algorithm is based ultimately on the total 
acceptable protection requirement. However, it was also suggested in 
WARC-77 that the single-entry protection requirement can be used as a 
guide for determining sharing criteria; of course the total interference 
from all sources still must be calculated to evaluate the scenario fully 
[871. The extended-gradient and cyclic-coordinate search methods 
discussed in Chapter III are based on the aggregate C/I e values. In 
contrast, the methods in this chapter, using the As concept, are based 
on single-entry C/I e protection requirements. Clearly such an approach 
will be acceptable only if there exists some relationship between the 
single-entry and aggregate C/I e values. 

Although so far no such relationship has been established 
rigorously, it is generally felt that the single-entry C/I e protection 
requirement does not need to exceed the total acceptable requirement by 
more than a few decibels. In Order to ensure tSiat a scenario made on 


the basis of single-entry C/I e protection requirement would result in 
acceptable aggregate C/I e values, WARC-77 suggested that the 
single-entry C/I e protection requirement be higher by 5 dB than the 
total acceptable protection requirement [78], Note that numerically 5 
dR is equal to 3; this suggestion is therefore based on the assumption 
that at the test point which has the worst aggregate C/I e value the 
aggregate interference will not exceed three times the strongest 
single-entry interference. This assumption is supported by the 
characteristics of the satellite transmitting and ground receiving 
antenna reference patterns. Referring to Figures 8.4, 2.5 and Equation 
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(2.6), because these two reference patterns are highly directional, 
among all the received single-entry interference powers only the few 
that come from the main or near side-lobe of the transmitting reference 
patterns and are received in the main or near side-lobe of the receiving 
reference pattern are relatively strong; the others are generally weak 
enough to be negligible. However, it is difficult to prove rigorously 
that 5 dB extra protection requirement is absolutely enough to cover the 
difference between the single-entry and aggregate C/I e values; hence the 
aggregate C/I e results must be calculated to evaluate the feasibility of 
a scenario. In the numerical examples in this chapter, 5 dB extra 
protection was used for all service areas. 

In any case, the methods to be described below do not depend on the 
validity of 5 dB, or any universal number. If the final analysis by the 
streamlined SOUP program shows more than 5 dB extra protection is needed 
for some service areas, then a more appropriate value may be chosen to 
compute the As ( 1 ) values for an improved synthesis. 


E. PERMUTATIONAL ASPECT OF THE ORBITAL-ASSIGNMENT PROBLEM 

With a data base of As (1 ) values for all the service area pairs, 
conceptually the planning problem could be solved by choosing a proper 
ordering or permutation of the satellites and completing the scenario 
by locating every satellite in its feasible orbital range while making 
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sure that the as (1 ) constraints are satisfied. As an example, consider 
a four-service area case with the geographic relation as shown in Figure 
4.4 with the objective of finding a scenario that uses the least orbit 
resource. The satellites are denoted as $a, Sr, Sq and Sr. Assume that 
the values of asrq(1) and asrd(1), the required satellite spacings for 
the non -adjacent service areas, are constant and are two degrees; 
similarly assume that the as(1) for all the adjacent service areas are 
constant and are four degrees. A scenario that meets this objective 
would be the ordering Sr-Sr-Sq-Sa with minimum required satellite 
spacing because it requires only eight degrees of orbital arc. 
Arrangements that do not include both the Sr-Sq and S B -S D satellite 
adjacencies would have to require at least a ten-degree arc. In this 
example, with the objective of conserving orbit resource, the basic 
concept is to have adjacent satellites serve non-adjacent service areas; 
here it utilizes the service-area separation to reduce the need for 


satellite spaci ng , 


arid achieves the maximum orbit utilization suggested 


in reference [16]. 
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Figure 4.4. Geographic location of four service areas. 


110 


In such a brute-force procedure, it is likely that a feasible 
solution can not be found for a particular ordering. It could be that 
the available orbital arc for the task is used up before the allocation 
of all the satellites has been completed. Then, another ordering has to 
be tried. 

Note that there are m! possible permutations for m satellites. 

This number becomes astronomical when m is large. If the goal is to 
find a scenario which is optimal by some criterion, in principle all 
the m! permutations have to be tested. Even if each test is simple and 
fast, the overall workload is still enormous. 

F. ORBITAL ASSIGNMENT OPTIMIZATION FORMULATIONS 

1. as(1) Constraint and Objective Function 

The orbital assignment problem can now be formulated as the 
optimization of a yet unspecified objective function, subject to the 
As ( 1 ) constraints on pairwise satellite spacings. The protection 
requirements will be satisfied because of the constraints. Also, it is 
highly desirable if the As ( 1 ) values are linear functions of the orbital 
variables so that a simple optimization technique, e.g., linear 
programming, can be used. This is apparently not true from Table 4.2 
and Figure 4.3. Still, the As ( 1 > values can be approximated by 
piecewise linear functions, or the maximum value of As ( 1 ) within a given 
orbital range, denoted as AS, can be used as a constant parameter in the 
optimization formulation. The objective function is, in principle, 
entirely arbitrary, but we shall restrict it to a linear function of the 
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orbital locations since the whole purpose of this approach was to avoid 
the computational complexity of nonlinear optimization. 

The objective function used in the numerical examples below is the 
sum of the absolute deviations of the assigned satellite positions from 
an arbitrary prescribed set of such positions, an "ideal" set. Such an 
"ideal" set might arise from requests of the user administrations. 
Alternatively, it may be specified as a tool for the synthesis, e.g., if 
the westernmost end of the available orbital arc is selected as the 
"ideal" location for all satellites, the resulting scenario is likely to 
allow the insertion of additional satellites at the eastern end at a 
later time with a minimum of readjustment. 

Other objective functions which have been proposed for minimization 
are the length of the occupied orbital arc and the constant zero. The 
latter simply seeks to find a solution which satisfies the as(1) 
constraints. 

2. Mixed-Integer and Restricted-Basis Linear Programming 
Formulations 

% 

With the single-entry C/I e protection requirement enforced by the 
As (1 ) constraints and a linear objective function, the problem can be 
formulated as a mixed integer program (MIP) [88]. Either the piecewise 
linear as functions or the constant AS parameters might be used to 
formulate this program; in this chapter the constant AS parameters will 
be used. The set of satellite locations which satisfy the AS 
constraints and the feasible orbital range constraints constitutes the 
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feasible region. When the problem is solved via the branch-and-bound 
algorithm, the globally optimal solution is guaranteed [88]. However, 
the computational effort required to find this solution can be 
prohibitively long when the problem size, i.e., the number of 
satellites, is large. 

To ease the burden on computational effort, the same problem can be 
formulated as a linear program (LP) with a set of nonlinear side 
constraints; only the constant AS parameters can be used to formulate 
this program. A linear program is much more readily solvable than a 
nonlinear program or an integer program, and is most often solved by the 
simplex method [81]. However, in this problem, the nonlinear side 
constraints prevent one from using the simplex method in its most common 
form. The method can be modified to handle these nonlinear side 
constraints through the use of restricted basis entry. In doing so, one 
is certain to find a local, but not necessarily a global, optimum. The 
computational effort required by the LP technique should be acceptable 
when the problem size is large. 

The MIP and the restricted-basis entry LP (RBLP) formulations are 
given in Appendix E as formulated by Professor Charles H. Reilly of the 
Department of Industrial and Systems Engineering for both of the 
non-trivial objective functions discussed above. 

G. NUMERICAL EXAMPLES 

1. Definition of the Problem 

In this example, a model of six administrations in South America, 

as shown in Figure 4.5, is used; they are Argentina (ARG), Bolivia 
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(BOL), Chile (CHL ) , Paraguay (PRG), Peru (PRU), and Uruguay (URG). The 
feasible orbital ranges for all the administrations are taken to be from 
80 to 110 degrees west. It is assumed that all satellites have the same 
frequency assignments. The antenna reference patterns are from Figures 
2.4 and 2.5. The optimization requirement is to minimize the sum of the 
deviations of the actually assigned positions from the satellite 
preferred locations. Three sets of satellite preferred locations are 
considered, as listed in Table 4.4: in case 1, all the preferred 

satellite locations are at the center of the feasible range; in case 2, 
all the preferred satellite locations are at the western boundary; in 
case 3, every administration has its preferred satellite location at a 
longitude for which the azimuth angle of the satellite from the 
administration center is close to zero. 

First, all the As(l) values must be calculated. They were 
£.j"i at gd using a 3Q dB singls-sntrv protection recjui rement • the 
results are listed in Table 4.2. The AS values of the 
six-administration problem, obtained from Table 4.2 by choosing the 
maximum As(l) values over the feasible range 80 to 110 degrees, are 
listed in Table 4.5. They are denoted as AS^j for satellite i and j in 
Appendix E. Using the AS parameters instead of the approximate 
piecewise linear as functions results in a conservative design (C/I e 
will tend to be larger), at the expense of possibly not using the orbit 
resource with maximum efficiency; but the alternative, i.e., using the 
As functions instead of the AS parameters as constraints, would 
complicate the MIP formulations; note that the piecewise linear 
constraints can not be used in the RBLP formulation. 
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Table 4.4 


Satellite preferred locations of six 
administrations 


preferred 

ARG 

BOL 

CHL 

PRG 

PRU 

URG 

case 1 

95 

95 

95 

95 

95 

95 

case 2 

110 

110 

110 

110 

110 

no 

case 3 

87.5 

92.5 

97.5 

87.5 

102.5 

82.5 



Table 4.5 

AS parameters of six admini 

st rations 

AS 

ARG BOL CHL PRG 

PRU 

URG 

ARG 

* 4.17 4.19 4.32 

1.41 

4.14 

BOL 

* 4.57 4.04 

4.26 

0.94 

CHL 

* 2.00 

3.94 

1.59 

PRG 

* 

1.10 

2.46 

PRU 


* 

0.37 

URG 



* 
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2. NIP and RBLP Results 


The solutions of the MIP and RBLP fomulations for this objective 
function, i.e., Formulations III and IV in Appendix E, are listed in 
Table 4.6 together with the values of the total deviations, the occupied 
orbital arcs and the computer run times in seconds. These data were 
provided by Professor C. Reilly and Mr. D. Gonsalvez of the Department 
of Industrial and Systems Engineering. 

Because of the AS constraints, the solution is guaranteed to 
satisfy the single-entry C/I e protection requirements; still the 
aggregate C/I e values need to be checked. By assuming that the term 
l/(fn 2, [2[F(f n ,fm)/f m 2 ]} equals one in Equation (4.3), this becomes a 
co-channel interference calculation. Thp annrpnatp r/T„ rpcuitc nf tho 

_ ’ — — J J m ~ J / - C * — - — ~ * vnu 

MIP solution of case 1 (which happens to be the most densely packed 

solution), calculated from the streamlined SOUP code in Appendix A, are 
listed in Appendix F. The results show that of the total 54 test 
points, only three places in Chile, three places in Paraguay and one 
place in Peru have aggregate C/I e values between 27 to 30 dB while all 
the rest of the C/I e values are above 30 dB. Note that the AS values 
(Table 4.5) were calculated from single-entry C/I e protection 
requirement of 30 dB, with the objective of guaranteeing that the 

aggregate C/I e values be no less than 25 dB (see the discussion in 

Section IV. D). Clearly this objective has been achieved; in fact, the 
aggregate C/I e results are better than expected. 

A likely reason for the aggregate C/I e values to be better than 
expected is the following. First note that in this example every 
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Table 4.6 


Mixed integer and linear program results 



case 

1 

case 

! 2 

case 

3 

satellite 

MIP 

LP 

MIP 

LP 

MIP 

LP 

AR6 

88.68 

105.74 

101.35 

110.00 

88.76 

101.26 

BOL 

99.57 

101.57 

97.18 

104.33 

92.93 

92.50 

CHL 

95.00 

97.00 

105.54 

99.76 

97.50 

97.07 

PRG 

93.00 

95.00 

107.54 

97.76 

84.44 

87.50 

PRU 

91.06 

93.06 

109.63 

108.59 

102.50 

102.67 

URG 

96.59 

92.54 

110.00 

105.86 

81.98 

82.50 


deviation 

18.42 

23.71 

28.76 

33.69 

5.27 

14.36 

arc 

10.89 

13.20 

12.82 

12.24 

20.52 

20.17 

cpu(sec)* 

25.23 

1.31 

13.39 

1.30 

2.86 

1.25 


* IBM-3081 computer 
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administration requires only one satellite. For a satellite and its 
corresponding service area, usually only the first adjacent satellites 
on both sides are close enough to produce significant interference, and 
these satellites do not usually influence the same test point; 
interference from far-away satellites are in general negligible because 
the actual satellite spacings are much larger than the required minimal 
spacings. This can be seen from the C/I e results in Appendix F: among 

the seven test points with C/I e less than 30 dB, only two receive two 
strong and almost equal signal, i.e., the two in Chile with C/I e of 
27.52 dB and 28.37 dB, and none receive more than two. Therefore the 5 
dB margin between single-entry and total acceptable protection 
requirements is adequate for these sample problems; this confirms the 
argument made in Section IV. 0. 

3. Comparison Between MIP and RBLP Techniques 

Referring to Table 4.6, note that the total deviations in the MIP 
solutions are smaller than those of the RBLP solutions. This is not 
unexpected since the MIP solution guarantees a global optimum with 
respect to the objective function while the RBLP formulation does not. 
On the other hand, the computer run times for the MIP formulation are 
significantly longer, and they are known to increase more rapidly with 
problem size than is the case with the RBLP formulation. As stated in 
Section IV. F. 2, the computer run time could be prohibitively long for 
the MIP formulation when the problem size, i.e., the number of 
satellites, is large. Therefore, the RBLP formulation becomes more 
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attractive as a practical method of solving a real-world problem. The 
MIP formulation Is useful primarily for evaluating the performance of 
the RBLP formulation on small problems. 

4. Suggestions for Further Improvement 

The C/I e ratios in Appendix F suggest that a uniform 5 dB margin 
between the single-entry and total acceptable C/I e protection 
requirements may result in over-protection. This raises the concern 
that, for a limited orbit resource and a large number of requests, a 
satisfactory scenario might never be found if the orbit resource barely 
allows every service area to have only the threshold aggregate C/I e 
value. 

When a feasible solution of the MIP or RBLP formulation does not 
exist for a given set of AS constraints, the reason could be either that 
the AS constraints are too high, or that there exists no solution that 
could satisfy the total acceptable C/I e protection requirement. To find 
out which, a progressive testing process could be used by gradually 
decreasing the C/I e requirement level in the as ( 1 ) calculation, and 
using these values in the MIP or RBLP calculation. Note that as the 
single-entry C/I e requirement is gradually decreased, the first point at 
which a feasible solution (feasible in terms of the AS requirements) 
exists may be such that the total acceptable C/Ie requirement is not 
satisfied at some test points. However, this does not mean that a 
feasible solution (feasible in terms of total acceptable C/I e 
requirement) does not exist; this is because the margins between the 
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single-entry and total acceptable C/I e protection requirements may not 
have to be the same for all satellite pairs. When this is the case it 
is proposed that the gradient search method be used to fine-tune the 
solution and improve those unacceptable aggregate C/I e values; this 
approach is discussed below. If no solution is obtainable when the 
As (1 ) values are calculated using the total acceptable C/I e protection 
requirement and the AS parameters are the minimum (over 1) of the As(l) 
values, a feasible scenario definitely does not exist. 

To demonstrate how the gradient search method can fix the worst 
C/I e terms, the resulting satellite locations of the MIP solution of 
case 1 are laid out as shown in Figure 4.6; here B is for satellite of 
Bolivia, 1) of Uruguay, C of Chile, P of Paraguay, E of Peru, and A of 
Argentina. Indicated in the figure are the actual spacings (above the 
arrows) and the AS values (below the arrows) between a satellite and its 
first and second adjacent satellites. The star (*) sign means that the 

dCtus! SpdCi ng is lApngp th^ri thp rnrrp^r»nn(ji nn AS vsIhp Pi ret not a 

that there is no star sign for the Chile satellite, while other 
satellites have at least one star. This means that the spacings between 
the Chile satellite and its first and second adjacent satellites are all 
at their minimal required values. This may explain why the C/I e results 
for Chile are worse than the others. Then note that this scenario may 
be modified by moving the Peru satellite to the right (eastward) and 
moving the Uruguay satellite to the left (westward) by a small amount. 
Although this modification increases the sum of the absolute deviations 
from the "ideal" locations, it also improves the C/I e results for Chile 
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Figure 4.6. Satellite locations from mixed integer program result. 



by increasing the satellite spacings for the Chile satellite; note that 
the Peru and Uruguay satellites are its east and west adjacent 
satellites. This demonstrates that the gradient search technique can be 
used on the MIP or RBLP solution to improve the aggregate C/I e results. 

5. Possible Extensions of the Method 

The AS concept has many flexibilities. It was demonstrated here 
for elliptical satellite antenna beams, but it can also be applied to 
shaped beams if the shaped beam reference pattern is given. At present, 
such a pattern is not available. Also, to use the as approach the 
antenna reference patterns and the C/I e protection requirement need not 
be the same for all administrations. Such non -uniformity merely changes 
the interference calculations, margin calculations, and the resulting 
as(1) values, but the same optimization procedures are still applicable. 

The formulations may be useful when new satellites need to be added 
into an existing scenario in which the locations of existing satellites 
can not be changed. The computational burden of either technique 
depends greatly on the number of decision variables. For a task of 
adding more satellites, all the Information about the existing scenario 
constitutes fixed parameters, and the only decision variables correspond 
to the new satellites. Therefore, the problem size is small, and the 
computational burden is reasonable. 

This chapter has dealt only with the orbital assignment, not the 
frequency assignment. This is useful for the case that every satellite 
uses the full complement of the available spectral band. Still, the 
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same procedure should work if there is an a priori frequency assignment 
scenario; a Considerable body of literature exists on such frequency 
assignments [19,20,21,22]. The only modification is to use the proper 
protection requirement value (by referring to the frequency assignments 
and protection ratio) in the as(1) calculation; the planning process 
will be exactly the same. 

H. CONCLUSION AND RECOMMENDATION 

In this chapter, two techniques are presented to solve the orbital 
assignment problem. The MIP technique guarantees to find the globally 
optimal solution, but may require prohibitively long computing time when 
the problem size is large. Still it is very useful for testing other 
methods on small problems and may be applicable when a few satellites 
are added into an existing scenario. The RBLP technique guarantees to 
find a locally, but not necessarily a globally optimal solution, but is 
more practical in terms of the computational effort. 

Two approximations are used in these methods. The first one 
appears in the as(1) calculation; however, it is shown to be acceptable. 
The second one is to use maximum As(l) values, AS, in the MIP and RBLP 
formulations; this decreases the efficiency of orbit use, but probably 
not seriously. 

It is recommended that these methods be tried on larger scenarios, 
and that the extensions given in Section IV. 6. 5 be investigated. 
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CHAPTER V 


A SERVICE-AREA SPECIFICATION PROCEDURE 


A. INTRODUCTION 


In the previous chapters the optimization of orbital and frequency 
assignments was considered for communications satellites serving a given 
set of service areas defined by political boundaries. Now, procedures 
for choosing such service areas will be discussed in the fixed-satellite 
service (FSS) context, with the objective of increasing the 
communications capacity for all users. The need to study ways in which 
the concept of service areas should be defined was brought up in CCIR 
report 453-3 [89], and the study in this chapter is intended to give 
insight to that concept. By means of an illustrative example it will be 


i ^ U ^ ^ V W A A A A ^ ^ A 1 AJrt A 4* lA A ^ f 4 A .4 /t m A 1A 4 4 A A A 4 4- w A 1A J A 4 4 ^ 
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allowable spacing (AS, see Chapter IV) can serve as a basis for 
designing service areas for this objective. 


B. HISTORICAL PERSPECTIVE 

1. Service Area Assignment 

In the broadcasting-satellite service (BSS) planning, a satellite 
service area is usually specified by the territories of an 
administration, a subdivision of an administration, or (in some cases) a 
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grouping of administrations [9,13], Thus there appeared to be no need 
for a procedure for defining service areas for the WARC-77 and RARC-83 
deliberations. 

The situation is quite different with respect to FSS. The FSS 
satellites were initially used primarily for inter-administration 
communications. International consortia operated such satellites, with 
service areas chosen on the basis of demand (or market), technological 
considerations, and the composition of the consortium. The Intelsat 
satellites are examples of this type of operation. More recently, 
satellites have also been employed for intra-administration (Domsat) 
communications, e.g., satellite communications in U.S., Canada and the 
U.S.S.R. (if the latter is considered as a single entity; technically 
the International Telecommunication Union (ITU) regards each separate 
Soviet Republic as an administration) [90]. The simplest technology 
leads to large service areas; e.g., the whole contiguous continental 
U.S. (CONUS). While an experimental FSS satellite using regional and 
switched spot beams is at the heart of the NASA Advanced Communications 
Technology Satellite (ACTS) program [91], most civilian U.S. operational 
and planned systems are based on CONUS coverage. In current approaches 
to the orbital and frequency assignments for the FSS, the idea of 
specifying service areas by the territories of the administrations seems 
to be still prevalent. 
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2. Insufficiency of Communl cations Supply from Limited Spectrum 
and Orbit Resources 

As FSS planning proceeds, it has become clear that the 
communications demand is very large [92,93]. As an example, the U.S. 
fixed-service communications supply appears to be in great mismatch with 
the demand [80,94]. The projected U.S. domestic FSS communications 
demand in the year 2000, as estimated in a COMSAT Laboratory report, is 
shown in Table 5.1, and is also shown as a "pie chart" in Figure 5.1 
[94]; the total communications demand is equivalent to 2,474 
transponders. 

By international agreement the total bandwidth of the 6/4 GHz 
spectral band assigned to the FSS is 480 MHz [94,95] (note that at 
present this is the only band in the planning stage). This band is 
divided into 12 slots, each consisting of a frequency channel 36 MHz 
wide followed by a guard band of 4 MHz. Therefore, including the 
cross-polarization channels, there are 24 frequency channels in this 
band [94,95]; this is shown in Figure 5.2. A typical U.S. FSS satellite 
has 24 transponders transmitting the signals in these 24 channels. 
Therefore the total estimated demand of 2,474 transponders implies a 
need for 103 satellites to fulfill the demand with the 6/4 GHz band. 

Because of interference between satellite systems, satellites must 
be separated from each other to have proper protection [16,96,97]. For 
the antenna technology available up to the year 1982, satellites with 
CONUS beams for systems proposed at that time needed to be separated by 
more than three degrees [98]. However, the amount of communications 
capacity so provided in the 6/4 GHz band is far from enough to satisfy 
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U.S. Domestic FSS Communications Demand in the Year 2000 
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Customer Premises Service 


TRANSPONDER ALLOCATIONS 


(US Domestic Traffic) 


CPS (510) 



(142) 


TRUNKING (1822) 


Figure 5.1. Projected U.S. Domestic FSS coommuni cations demand for the 
year 2000. CPS stands for Customer Premises Service. The 
number in parentheses indicates number of transponders. 
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the FSS demand [99]. In order to increase the communications supply, 
the Federal Communication Commission (FCC) decided that the orbit 
resource should be more efficiently utilized and announced in the year 
1984 that future antenna technology should allow two-degree spacing for 
U.S. FSS satellites, and satellite orbital planning should be based on 
two-degree spacing after the year 1987 (This also applies to the 14/12 
GHz band) [80]. This will enable the U.S. to have about 25 
geostationary satellites in the 6/4 GHz band. Clearly this still falls 
far short of the projected demand of 103 satellites! 

Other spectral bands, e.g., the 14/12 GHz band, have been allocated 
to the FSS [31]. This should help increase the U.S. FSS communications 
supply. Still, it is apparent that the spectrum resource is limited, 
and that the orbit resource should be more efficiently utilized in order 
to maximally re-use the spectrum resource. 

3. Relation Between Frequency Re-Use and Service-Area 
Sped fi cation 

The advantage of using narrow antenna beams to achieve frequency 
re-use has been recognized for some time [89,100]. Thus it was noted in 
CCIR report 453-3 that it is necessary to study ways in which the 
concept of coverage area should be defined [89], Specifically, It was 
brought up in this report that under certain circumstances a satellite 
may transmit separate information on the same frequencies twice, or even 
a greater number of times, using antennas serving different parts of the 
world. This general idea has already found practical application with 
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the Intelsat series of satellites; while Intelsat IV satellites used 
global -coverage beams, Intelsat IV-A and V satellites employ 
hemispherically restricted beams to allow frequency re-use [101,102]. 

In view of the scarcity of the spectrum/orbit resources relative to 
demand, it becomes necessary to develop a method of specifying service 
areas such that the advantages* of the narrow-beam concept can be fully 
utilized for efficient spectrum/orbit management. 

C. SERVICE-AREA SPECIFICATION BY AS CONCEPT AND COMMUNICATIONS-DEMAND 
DENSITY 


1. Role of AS Concept In Service-Area Specification 

In order to provide more traffic supply, the AS concept discussed 
in Chapter IV suggests that the AS values between satellite pairs should 
be as small as possible so that more satellites can be allocated in the 
orbit. As discussed in detail in Section IV. C, the margin between the 
carrier power, C, and the single-entry interference, I, comes mainly 
from the transmitting and receiving antenna discrimination factors in 
the interference power; the former comes from the separation between the 
two service areas, the latter from the separation between the two 
satellites. Also, as discussed in Section IV. C, for a given 
single-entry C/I protection requirement, the larger the transmitting 
discrimination loss, the less the receiving discrimination loss and thus 
the less satellite separation is needed. For the special case when the 
transmitting discrimination loss is enough for the protection 
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requirement to be met, the satellites can be collocated, or equivalently 
a satellite can have two beams serving two service areas simultaneously; 
this is exactly what CCIR report 453-3 means by using the same frequency 
band twice at one satellite [89], 

Referring to Figure 2.4, the transmitting discrimination value 
depends on two factors: the half -power beam width (HPBW) and the 

off-axis angle. The value of HPBW depends on the size of the service 
area: the smaller the service area, the smaller the value of HPBW, and 
the narrower the beam; this is the narrow-beam idea stated in references 
[89,100]. The off-axis angle depends on the service-area separation. 

The combination of these two factors should be the key to the subject of 
service-area definition. If the system planner can control these two 
factors to reduce the AS values between satellite pairs, more satellites 
can be assigned in the orbit and a larger communications supply; i.e., 
more circuits, can be provided. To achieve this purpose, basically the 
service areas should be specified as small as possible, and their 
separations should be as large as possible. 

2. Role of Communications -Demand in Service-Area Specification 

The simple demand/supply concept implies that the supply should be 
used where the demand is; hence the communications -demand density must 
also be consulted in specifying service areas. There is precedent for 
using communications traffic-demand density quantitatively in the 
technical design of satellite systems. In time-division multiple access 
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(TDMA) system design, a simple rule is: the areas with heavy 

communications demand should be given more access to the communications 
system [103]. Similarly, service areas should specify regions with 
enough traffic demand to justify having satellite access as an area. 
Hence information about the communications traffic-demand density is an 
essential prerequisite to a reasonable choice of service areas. 

3. General Consideration 

Hence, the AS concept and the communications traffic-demand density 
should both be considered when specifying service areas. It is apparent 
that small and separated service areas allow the geostationary orbit to 
be used more efficiently; on the other hand, potential demand increases 
with the service-area size: a satellite which serves one metropolitan 

area (e.g., Boston) has less potential demand than one which serves a 
corresponding region (e.g., the Eastern U.S.) or the entire 
administration or a grouping of administrations. There has to be a 
compromise between these two factors. The overall rule is that the 
selected service areas must have enough communications demand to justify 
their own satellite beams. 

There are, of course, other factors which need to be taken into 
account in a practical plan. For example, multiple-beam satellites, 
which may also be beam-switching, require more advanced technology and 
are likely to be more expensive than single-beam ones. Such economic 
and perhaps other operational matters will not be considered here; only 
the communications capacity will be addressed. 
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D. APPLICATION OF SERVICE-AREA SPECIFICATION CONCEPT TO A LARGE 


ADMINISTRATION OR GROUP OF ADMINISTRATIONS 

1. General Description 

The idea of service-area specification by the AS concept and 
traffic-demand density may greatly benefit an administration, or a group 
of administrations, that has a large territory, when most of the traffic 
demand is between several small regions. When this idea is applied to 
such a situation, it is proposed to serve such administration(s) with a 
mixture of administration-coverage beams and intra-administration 
regional beams. The administration -coverage beams are intended to serve 
primarily areas of relatively low demand, while the majority of traffic 
demand from high-demand regions would be carried on the regional beams. 
The necessary separations between regional -beam satellites are usually 
smaller than those of the admini strati on -coverage satellites. Thus if, 
instead of assigning the entire available orbital arc to 
administration-coverage satellites, a part is used to accommodate 
regional -beam satellites, then more satellites can be allocated in this 
arc and more communications supply can be provided. 

Recall the scenario which has 25 CONUS beam satellites in the 6/4 
GHz band with every satellite using the full spectral band (Section 
V.B.2). At first look it seems that this scheme uses the spectrum and 
orbit resources to their full extent for Earth station antenna 
technology which requires two-degree satellite spacing. However, if 
some of the service areas are changed from CONUS to smaller, separated 
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regions, the AS values for their satellites will be smaller than two 
degrees, hence some orbit resource becomes available to allocate more 
satellites in the orbit. 

2. Procedure of Service-Area Specification and Satellite 
Assl gnment 

When regional service areas are to be specified, first the 
traffic-demand density should be considered. Because of the population 
distribution or the commercial activities, usually there will exist some 
regions of relatively small size and large traffic demand. Often such 
regions Include several large cities relatively close to each other. 

Such regions are good candidates for regional service areas within the 
admlni strati on (s). In selecting these service areas, three things 
should be considered simultaneously: the areas should be of small size, 

they should have large traffic demand, and the separations between these 
service areas should be large. 

After the regional service areas have been selected tentatively, 
the AS values between these service area pairs should be calculated by 
the method described in Chapter IV. Note that service areas for which 
the AS values are zero can be served by a single satellite; this is the 
multiple-beam design described in CCIR report 453-3. Different from the 
multiple-beam TDMA design, in this case the two beams can be active 
simultaneously [104]. Also note that multiple-beam satellites can carry 
both Inter- and intra-regional communications. 

A traffic-demand matrix can be formulated for these tentative 
service areas. The elements in this matrix are traffic demand between 
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the regions that can be carried by the regional -beam satellites: the 

diagonal elements are the intra-regional taffic demand, the off-diagonal 
elements are the inter-regional traffic demand. The latter are non-zero 
only when the corresponding service areas can be served by multiple-beam 
satellites. A complete demand matrix consists of the above-described 
matrix plus one additional element: the traffic demand which cannot be 

carried by regional satellites and must be carried by 
admi ni strati on -coverage satel 1 ites . 

Examination of the AS and demand matrices may suggest a revision of 
the regional service areas. For example, if many of the AS values for a 
particular region are large (e.g., approaching the AS value for 
administration-coverage satellites), this may indicate the region is too 
large, or too close to other regions, or both. On the other hand if 
many of the demand matrix elements corresponding to a region are very 
small, this may indicate the region is too small, or it is not a good 
candidate to be a regional service area. Thus a good choice of regional 
service areas becomes a compromise between achieving satisfactory AS 
values and a satisfactory demand matrix, as will be evident from the 
example below. 

A scenario including both regional and admi ni strati on -coverage 
beams may, in principle, be constructed by any suitable method, e.g., by 
extensions of the methods of Chapters III and IV. Here a scenario will 
be generated by deassigning some of a series of equally spaced 
administration-coverage satellites and then allocating regional -beam 
satellites in this vacated orbital arc, consistent with the AS 
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requirement. The number of administration -coverage satellites to be 

x. 

deassigned and the number of regional beams to be assigned are 
determined by consulting the traffic-demand matrix. This procedure has 
the advantage of being compatible with a network of uniformly spaced 
administration-coverage satellites; the new scenario might even evolve 
out of such a system. It also has the advantage of demonstrating the 
important concepts without the computational complexity of 
optimization. 

3. Traffic Distribution Between Narrow- and Hide-Beam Systems 

The traffic-demand problem is completely solved when the amount of 
traffic-supply is enough, or more than enough, to meet the demand. 
However, as stated earlier, the demand usually far exceeds the supply, 
and usually the adoption of the regional -service area idea can only 
improve, but not completely solve the demand/supply problem. Therefore, 
an important task of the procedure is to make sure all the demands have 
their proper share of the supply. 

It will be assumed that the satellites serving the regional service 
areas are dedicated satellites which are not designed to provide 
administration-wide service. Therefore some satellites need to be 
preserved to provide administration-wide communications service, even 
though the traffic demand which can not be served on a regional basis 
may be small. Many criteria could be used to decide the distribution of 
satellite beams. It is known that, compared to a scenario with only 
administration -coverage beams, there will be a larger communications 
supply when the narrow-beam idea is implemented. Therefore, one 
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criterion might be that all the demands have an equal percentage of 
supply-increase, another criterion might be that the total 
communications supply has the maximum amount of increase and no one 
suffers any supply decrease. The latter is used to develop the scenario 
in Section V.E. 

E. NUMERICAL EXAMPLE 

1. Description of Parameters 

Since the only available detailed traffic-demand density data 
available to the authors is the projected long-distance telephone 
traffic among major cities in the U.S., it will be used in an example to 
numerically explore the advantage of the service area specification 
concept. Because only the 6/4 GHz band is in the planning stage, this 
example considers the communi cations capacity only for this band. Also, 
only the down-link communications traffic regulation problem will be 
considered; the up-link problem can be implemented similarly in the 
up-link spectral band. 

Of the total projected FSS communications demand shown in Table 
5.1, more than 63% consists of long-distance telephone voice traffic 
[94]. The telephone voice demand will be taken as indicative of total 
demand. This assumption is justified in part because only telephone 
voice data is available to us, in part because it seems likely that the 
geographical distribution of other communications will be similar to 
that of telephone voice traffic, and in part because the objective here 
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is only to demonstrate a method, not to design a system. More detailed 
examination shows that a large portion of the voice long-distance 
telephone demand comes from major cities (or metropolitan areas), and 
there are 28 cities that share more than 75% of the voice long distance 
telephone demand. The list of these 28 cities, and the projected long 
distance voice telephone demand between them are shown in Table 5.2 
[94]. The U.S. map in Figure 5.3 shows the locations of these 28 cities 
with their communications rankings: the cities indicated by triple 

circles are in the top 5 rankings, double circles are for rankings from 
6 to 10, single circles for rankings from 11 to 15, and solid dots for 
rankings from 15 and up. 

In order to arrive at AS values and C/I ratios consistent with the 
antenna patterns shown in Figures 2.4 and 2.5, the C/I results of two 
CONUS beam satellites with 2.5-degree spacing were calculated and are 
listed in Table 5.3, and the worst single-entry C/I value is 25 dB. 

This value will be used as the single-entry C/I protection requirement 
for the U.S.; the total acceptable C/I protection requirement will then 
be 20 dB if 5 dB extra protection is needed to compensate for multiple 
interference. It is not implied here that 2.5-degree spacing or a 20 dB 
aggregate protection ratio is recommended. The purpose here is only to 
arrive at a consistent set of parameters for demonstrating a 
regional -coverage assignment procedure. 

The calculations in this section are on a co-channel basis. Since 
adjacent channels in the scheme of Figure 5.2 are cross-polarized, it is 
assumed that there is enough polarization and frequency discrimination 
that this is a satisfactory approximation. 
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Long Distance Voice Telephone Demand Between Top 28 Cities in U.S. in the Year 2000 
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Figure 5.3. Locations of top 28 cities in U.S. with largest 
long-distance telephone demand in the year 2000. 


Table 5.3 


C/I Results of Two CONUS-beam satellites with 
2.5-degree spacing 


COUNTRY 

SATELLITE 

(LON. ) 

FREQUENCY ( 

USA 

-100, 

,00 


4000.00 

USA 

-102. 

,50 


4000.00 

TEST COUNTRY s USA 


SATELLITE 

: -100.00 

TEST 

LON. 

POINT 

LAT. 


INT. SAT. 

C/I (dB) 

-69.20 

47.40 


USA 

-108.35 

-66.90 

44.60 


USA 

-108.29 

-69.90 

41.50 


USA 

=107.71 

-81.80 

24.40 


USA 

-108.44 

-97.20 

26.00 


USA 

-107.61 

-117.10 

32.30 


USA 

-108.47 

-124.20 

40.40 


USA 

-108.50 

-122.80 

49.00 


USA 

-108 , 39 

-95.10 

49.40 


USA 

-107.22 

TEST COUNTRY s USA 


SATELLITE 

: -102.50 

TEST 

LON. 

POINT 

LAT. 


INT. SAT. 

C/I (dB) 

-69.20 

47.40 


USA 

-108.35 

-66.90 

44.80 


USA 

-108.30 

-69.90 

41.50 


USA 

-107.72 

-81.80 

24.40 


USA 

-108.41 

-97.20 

26.00 


USA 

-107.60 

-117.10 

32.30 


USA 

-108.51 

-124.20 

40.40 


USA 

-108.54 

-122.80 

49.00 


USA 

-108.41 

-95.10 

49.40 


USA 

-107.19 
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2. Specifying Four Service Areas According to AS Consideration and 
Traffic-Demand Density 

Referring to Figure 5.3 about the communications traffic-demand 
density, four service areas that cover most of the top 28 cities may be 
specified as shown in Figure 5.4. The first service area, called the 
East region, covers the cities New York (1), Boston (5), Washington D.C. 
(7), Philadelphia (9), and Raleigh (18). The second service area, 
called North Central, covers Chicago (3), Detroit (6), Cincinnati (8) 
and Cleveland (10). The third service area, called West, covers Los 
Angeles (2), San Francisco (4) and Anaheim (12). The fourth service 
area, called South Central, covers Dallas (11), Houston (14), San 
Antonio (25) and New Orleans (27). These regional service areas were 
obtained by trial and error as reasonable compromises between the AS 
matrix and demand matrix requirements. For example, the traffic demand 
for the North Central region could be increased by including Syracuse 
(15), St. Louis (17), Minneapolis (20), Kansas City (22), and Milwaukee 
(24); however, this would increase its size and reduce the separation 
between regions, thus it would increase aSjt^nc and aS NC,SC sufficiently 
to prohibit sufficient frequency re-use. Even as it is, aS^nc turns 
out too large to use a multiple-beam design. Its size could be reduced 
by eliminating Chicago (3) and Detroit (6) from this region, but then 
the traffic-demand might be too low to justify having its own satellite 
beams. This illustrates the compromise between AS and traffic demand 
which must be the basis in specifying regional service areas. 
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ea se 


The AS matrix of the four service areas is shown in Table 5.4: the 

AS values between any two of the four service areas, except between East 
and North Central, are zero. To demonstrate that a multiple-beam scheme 
is feasible, the C/I ratio of a three-beam satellite serving East, West 
and South Central is shown in Table 5.5, and the C/I ratio of a 
three-beam satellite serving North Central, West and South Central is 
shown in Table 5.6; in both cases, the worst aggregate C/I values are 
larger than 25 dB. The AS value between East and North Central is 1.24 
degrees, as indicated in Table 5.7; however, for numerical convenience 
1.25-degree spacing will be used for the East and North Central 
satellites. 

With this AS matrix, there can be two classes of regional -beam 
designs for satellites serving these four regions. In the first one the 
service areas East, West and South Central, or any combination of these, 
can use one multiple-beam satellite. In the second, the service areas 
North Central, West, and South Central, or any combination of them, can 
use one multiple-beam satellite. There are two overall restrictions: 
satellites serving the East and North Central regions respectively need 
to be separated by at least 1.25 degrees, and any two satellites that 
serve the same service area must be separated by no less than 2.5 
degrees. Aside from these restrictions, a scenario can have any 
combination of the two regional -beam designs. For example, two adjacent 
satellites that serve the regions East and North Central, respectively, 
can also serve the service areas West and/or South Central. 
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Table 5.4 


AS matrix of four regions 


AS 

East 

N. Central 

West 

S. Central 

East 

2.50 

1.25 

0 

0 

N. Central 


2.50 

0 

0 

West 



2.50 

0 


S. Central 


2.50 


Table 5.5 


'-C/I Results of collocating satellites serving 


East, West and South Central regions 


REGION 

SATELLITE 

(LON. > 

FREQUENCY (MHz) 


EAS 

-100.00 


4000.00 


WES 

-100.00 


4000.00 


SCR 

-100.00 


4000.00 


TEST REGION 

: EAS 

SATELLITE 

: -100.00 


TEST 

POINT 

I NT. SAT 

. C/I (dB) 

MARGIN 

LON. 

LAT. 




-74.02 

40.72 

WES 

44 . 30 


-74.02 

40.72 

SCR 

34 . 35 


-74.02 

40.72 

TOTAL 

33.93 

8.93 

-71 .07 

42.35 

WES 

42.97 


-71 .07 

42.35 

SCR 

34 . 43 


-71 .07 

42.35 

TOTAL 

33.86 

8.86 

-77.02 

38.90 

WES 

44.29 


-77.02 

38.90 

SCR 

32 .49 


-77 .02 

38.90 

TOTAL 

32.21 

7.21 

-75.12 

39 . 95 

WES 

44.54 


-75.12 

39.95 

SCR 

33 . 90 


-75.12 

39 .95 

TOTAL 

33.54 

8 . 54 

-73.65 

35.78 

WES 

43.28 


-78.65 

35.78 

SCR 

28.91 


-78.65 

35.78 

TOTAL 

28.76 

3.76 

TEST REGION 

: WES 

SATELLITE 

s -100.00 


TEST POINT 

INT. SAT 

. C/I (dB) 

MARGIN 

-118.25 

34.05 

EAS 

44.13 


-1:3.25 

34.05 

SCR 

32 . 33 


-118.25 

34.05 

TOTAL 

32.05 

7.05 

-122.40 

37.80 

EAS 

43.56 


-122.40 

37.80 

SCR 

33 . 60 


-122.40 

37.80 

TOTAL 

33 . 37 

8.37 

-117.95 

33.85 

EAS 

43.97 


-117.95 

33.85 

SCR 

32.02 


-117.95 

33.85 

TOTAL 

31 .75 

6.75 

TEST REGION 

: SCR 

SATELLITE 

: -100.00 


TEST 

POINT 

INT. SAT 

. C/I (dB ) 

MARGIN 

-76. S0 

32.78 

EAS 

34 . 90 


-95.80 

32 . 78 

WES 

37.02 


-9:.S0 

32.73 

TOTAL 

o 9 p 

7 . 82 

-95.52 

29.47 

EAS 

35.71 


-53.52 

29.47 

WES 

35.22 


-93.52 

29.47 

TOTAL 

32.45 

7.45 

-90.12 

29.97 

EAS 

28.14 


-90.12 

29.97 

WES 

39.74 


-90.12 

29.97 

TOTAL 

27.85 

2.85 


15(1 


<dB) 


(dB> 


(dB ) 



Table 5.6 


C/I results of 

collocati ng 

satellites 

serving 

North Central , 

West and South Central 

regions 

REGION 

SATELLITE (LON. ) 

FREQUENCY 

(MHz ) 

WES 

- 100.00 

4000.00 


NCR 

- 100.00 

4000.00 


SCR 

- 100.00 

4060.00 


TEST REGION 

: WES 

SATELLITE 

: -106.00 


TEST POINT 

INT. SAT. 

C/I (dB) 

MARGIN 

LON. 

- 1 I S . 2 5 

LAT. 

34.05 

NCR 

43.45 

-113.25 

24.65 

SCR 

32.33 


-1 IS. 25 

3 4 . 05 

TOTAL 

32.00 

7 .00 

-122. 40 

37 . 80 

NCR 

43.47 


-122.40 

37 . 80 

SCR 

33 . 30 


-122.40 

37.80 

TOTAL 

33 . 36 

8.36 

-117.35 

33.85 

NCR 

43.23 


-117.95 

33.85 

SCR 

32 . 6 r 2 


-117.95 

33.85 

TOTAL 

31.70 

6.70 

TEST REGION 

: NCR 

SATELLITE 

: - 100.00 


TEST 

POINT 

INT. SAT. 

C/I (dB) 

MARGIN 

-87.63 

41.88 

WES 

42.72 

-87.83 

41.88 

SCR 

30.48 


-87.63 

41.88 

TOTAL 

30.23 

5.23 

-81.68 

41.50 

WES 

43.22 


-81.68 

41.50 

SCR 

31.60 


-81.68 

4 1.50 

TOTAL 

3 1.50 

6 . 50 

-83.05 

42.33 

WES 

43.97 


-83.05 

42.33 

SCR 

32 . 84 


-83.05 

42 . 33 

TOTAL 

32.52 

7.52 

-84 .52 

39 . 10 

WES 

43.73 


-84 . 52 

39 . 10 

SCR 

29 . 55 


-84.52 

39 . 10 

TOTAL 

29 . 39 

4.39 

TEST REGION 

: SCR 

SATELLITE 

: - 100.00 


TEST 

POINT 

INT. SAT. 

C/I (dB) 

MARGIN 

-36.80 

32.78 

WES 

37.02 

-95 . 80 

32.78 

NCR 

31.08 


-36.80 

32.78 

TOTAL 

30.09 

5.09 

-35.27 

29,77 

WES 

36 . 67 


-95.37 

29.77 

NCR 

34.02 


-RE. 37 

23.77 

TOTAL 

3 2 . t c 

7 . 76 

-98.52 

29.47 

WES 

35.22 


-38.52 

29.47 

NCR 

34 . 37 


-58.52 

29.47 

TOTAL 

31.77 

6.77 

-30. 12 

29 . 97 

WES 

39.74 


-50. 12 

29.97 

NCR 

29.50 


-90. 12 

29. 97 

TOTAL 

29.11 

4.11 


<d B) 


(dB ) 


(dB) 
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Table 5.7 


As calculation for East and North Central regions 


REGION 

SATELLITE 

(LON. ) FREQUENCY (MHz) 


NCR 

-70.00 


4000 . 00 


EST 

-70.00 


4000.00 


TEST REGION 

: NCR 

SATELLITE : 

-70.00 


TEST 

LON. 

POINT 

LAT. 

I NT. SAT. 

C/I (dB) 

MARGIN 

-87.63 

41.88 

EST 

28.66 

3.66 

-83.05 

42.33 

EST 

14.71 

-10.29 

-84.52 

39. 10 

EST 

17.10 

-7.90 

-81.68 

41.50 

EST 

8.59 

-16.41 

TEST REGION 

? EST 

SATELLITE : 

-70.00 


TEST 

POINT 

I NT. SAT. 

C/I (dB) 

MARGIN 

-71.07 

42.35 

NCR 

29.00 

4.00 

-75.12 

39 .95 

NCR 

22.11 

-2.89 

-78.65 

35.78 

NCR 

12.64 

-12.36 

-74.02 

40.72 

NCR 

29.50 

4.50 

-74.02 

40.72 

NCR 

29.50 

4.50 

WORST MARGIN 

IS -16.41 

dB AT NCR 

( -81.68, 41 

.50) 


REQUIRED SATELLITE SPACING : 1.24 AT - 73.00 FOR C/I 25.0 dB 


152 


The projected long-distance telephone demand matrix between the 
major cities in these four areas is given in Table 5.8(a). When traffic 
demand for a region or between regions is computed below, we shall 
simply add matrix elements from Table 5.2. For instance, internal 
demand for the West region would be the demand from L.A. to San 
Francisco (2)-(4), San Francisco to L.A. (4)-(2), L.A. to Anaheim 

(2) -(12), Anaheim to L.A. (12)-(2), San Francisco to Anaheim (4)-(12), 
Anaheim to San Francisco (12)-(4). For demand from the West region to 
the North Central region, one would add the matrix elements 
corresponding to transmission from each of (2), (4), (12) to each of 

(3) , (6), (8), (10). To make these numbers more applicable to the 
satellite planning task, they are converted to a corresponding number of 
satellite beams in the following way. From Figure 5.2, there are 24 
frequency channels in the 6/4 GHz spectral band, and each channel is 36 
MHz wide. Assuming that a high quality telephone channel occupies 8 KHz 
bandwidth, one MHz bandwidth can carry 120 telephone channels, and hence 
a full spectral band can carry 103,680 telephone channels. The number 
of satellite beams is therefore obtained from the number of telephone 
channels by dividing by 103,680. Table 5.8(a) is re-stated as satellite 
beam requirements in Table 5.8(b). 

A complete requirement matrix, R, for the U.S. is shown in beam 
units in Table 5.9. In this matrix, the Inter- and intra-regional 
communications demand of the four regions that can use the two 
regional -beam designs are listed individually, the rest of the demand, 
i.e., 73.15 beams, that must go through the CONUS-beam satellites is 
listed in one category. Note that the inter-regional demand between 
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Table 5.8 

Telephone communications demand between major 
cities in four regions 
(a) In channel units 


channels 

East 

N. Central 

West 

S. Central 

East 

558400 

514000 

400200 

198700 

N. Central 

514000 

240400 

253300 

124300 

West 

400200 

253300 

124800 

97800 

S. Central 

198700 

124300 

97800 

31200 


(b) In beam units 


beams 

East 

N. Central 

West 

S. Central 

East 

5.38 

4.95 

3.85 

1.91 

Central 

4.95 

2.31 

2.44 

1.19 

West 

3.85 

2.44 

1.20 

0.94 

Central 

1.91 

1.19 

0.94 

0.30 
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Table 5.9 

Complete requirement matrix 


beams 

East 

N. Central 

West 

S. Central 

CONUS 

East 

5.38 

* 

3.85 

1.91 


Central 

* 

2.31 

2.44 

1.19 

73.15 

West 

3.85 

2.44 

1.20 

0.94 


Central 

1.91 

1.19 

0.94 

0.30 
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East and North Central is included in the CONUS beam category because it 
can not use any of the two regional -beam designs. 

3. Improvement of traffic-supply matrix with service-area 
specification 

A supply matrix is used to express the traffic-supply arrangement; 
each element in this matrix is the amount of traffic supply from one 
region to another offered from a scenario, up to the maximum amount of 
demand from the corresponding regions. 

If the communications demand is supported by 25 CONUS-beam 
satellites, the average percentage of satisfaction, s, is 
s = 25/103 
= 24.27% . 

Assuming that each user in the continental U.S. has equal probability of 
accessing all satellites, then each user will have the same percentage 
of satisfaction and the corresponding supply matrix in beam numbers is 
given by 

s = s • R , (5.1) 

which is shown in Table 5.10. This represents the total supply expected 
from the al 1 -CONUS system. While the origin and destination of traffic 
by regions East, N. Central, S. Central and West are not particularly 
meaningful with the all CONUS-beam design, the same partitioning as in 
Table 5.9 has been retained in Table 5.10; this will prove convenient 
when this scheme is compared to the multiple-beam designs to be 
discussed. 
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Table 5.10 


Supply matrix without regional -beam 
satellites 


beams 

East 

N. Central 

West 

S. Central 

CONUS 

East 

1.305 

* 

0.934 

0.463 


N. Central 

* 

0.560 

0.592 

0.289 

17.753 

West 

0.934 

0.592 

0.291 

0.228 


S. Central 

0.463 

0.289 

0.228 

0.073 
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When Implementing the two regional -beam designs, some of the 25 
CONUS beam satellites must be preserved; the criterion here Is to have 
the maximum total increase of traffic supply and no loss for any 
individual. The projected total demand is 103 beams. The requirement 
matrix that can use the two regional -beam designs, which is simply the 
regionally served part of Table 5.9, is shown in Table 5.11; the 
regional sub-total requirement is found to be 29.85 beams. Suppose 
there are m CONUS satellites being replaced by the regional -beam 
designs. If the taffic which can be carried only by the CONUS beams 
does not suffer any loss after implementing the regional -beam designs, 
the number m must satisfy the requirement 
( 25— m)/ ( 103-29.85) > (25/103). 

The maximum value of m, denoted by M, is 7. Replacing M consecutive 
CONUS satellites with regional -beam satellites will lead to the maximum 
total increase of traffic supply. For a general case, the value of M is 
calculated from 

M - IFIX[n-(n/t).(t-r)], ( 5 . 2 ) 

where the symbols are 

IF IX : a function whose value is the largest integer equal 
to or less than the argument, 

n : number of satellites in a scenario before implementing 
regional -beam designs, 
t : total communications demand, 

r : communications demand that is to be satisfied by the 
regional -beam designs. 


158 



Table 5.11 


Requirement matrix of four regions 


beams 

East 

N. Central 

West 

S. Central 

sub-total 

East 

5.38 

★ 

3.85 

1.91 

11.14 

Central 

* 

2.31 

2.44 

1.19 

5.94 

West 

3.85 

2.44 

1.20 

0.94 

8.43 

Central 

1.91 

1.19 

0.94 

0.30 

4.34 


total 29.85 
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Reference to the AS matrix in Table 5.4 shows the orbital arc that 
accommodates seven CONUS-beam satellites can support thirteen satellites 
with 1.25-degree spacing. Of these thirteen satellites, there could be 
maximally seven beams illuminating the region East since satellites 
serving the same service area must be separated by no less than 2.5 
degrees. Then, because of the AS restriction which does not allow 
collocation of East and North Central beams, there can be at most six 
beams illuminating North Central. (The reason that East has been 
allotted more beams than North Central is that it has more demand.) The 
maximum numbers of beams for West and South Central are both six. They 
can not be seven because the beam arrangement should allow 
inter-regional communications from East or North Central to West and 
South Central; the choice of seven West beams would result in 
beam-collocation of the East and West beams and no communications 
between West and North Central. (If these statements do not seem clear, 
the reader is encouraged to refer to Figure 5.5, which shows the 
sequence of the thirteen regional -beam satellites, and to experiment 
with other arrangements, keeping in mind both restrictions and 
regional -beam designs discussed in Section V.E.2.) With these beam 
assignments, the supply of regional beams to these four regions is 
listed in Table 5.12. 

The corresponding percentages of satisfaction for the four regions 
are listed in Table 5.13; the corresponding supply matrix, obtained by 
multiplying elements of each row with the respective percentage of 
satisfaction, is shown in Table 5.14(a). It is assumed that the 
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Table 5.12 


Supply of regional beams to four regions 

regi on beams 
East 7 

N. Central 6 

West 6 

S. Central 6 

total 25 
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Table 5.13 


Percentage of satisfaction of demand by 
means of regional beams 


region 
East 
N. Central 
West 
S. Central 


p. o. s. 

7/11.14 = 0.6283 
1 

6/8.43 = 0.7117 
1 
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satellite traffic is likely to be symmetrical, hence this matrix needs 
to be symmetrized. The adjustment should be dominated by the first row 
because the demand from/to East is the largest and the rate of 
successful connection depends therefore most strongly on the 
availability of the channels to the East. The first step is to 
symmetrize the first column and row, and the matrix should become as 
shown in Table 5.14(b). Next, the adjustment should be dominated by the 
third row because the communications from/to West is the second largest. 
Note that since the value of its first element has been decreased in 
step one, the available communications supply to the other three 
elements will increase. The total demand represented by the last three 
elements in row three can be seen from Table 5.9 to be 4.58 beam units. 

The beam supply corresponding to these elements are the six beams 

allotted to the West region (see Table 5.12) minus the 2.419 beam units 

assigned to the first element of row three in Table 5.14(b), or 3.581 

beams. The percentage of satisfaction is therefore 2.419/3.581, or 
0.7818, which means the matrix is adjusted as shown in Table 5.14(c). 
This adjustment procedure continues, but is less complicated for the 
remaining two rows because the amount of beam supply for North Central 
and South Central is larger than the demand. The adjustment is 
completed by satisfying the demand, i.e., inserting the corresponding 
elements from Table 5.11, and the final result is shown in Table 
5.14(d). At this point it should be recalled that the design began by 
converting seven of 25 satellites to regional -beam use, see discussion 
above Equation (5.2); thus 18 remain for CONUS coverage. For a clear 
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Table 5.14 

Supply matrix adjustment 
(a) Step one 


beams 

East N. 

Central 

West 

S. Central 

East 

3.380 

★ 

2.419 

1.200 

N. Central 

★ 

2.310 

2.440 

1.190 

West 

2.740 

1.736 

0.854 

0.669 

S. Central 

1.910 

1.190 

0.940 

0.300 




(b) Step two 


beams 

East 

N. Central West 

S. Central 

East 

3.380 

* 2.419 

1.200 

N. Central 

* 



West 

2.419 



S. Central 

1.200 
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Table 5.14 (Continued) 
(c) Step three 


beams 

East 

N. Central 

West 

S. Central 

East 

3.380 

* 

2.419 

1.200 

N. Central 

* 


1.908 


West 

2.419 

1.908 

0.938 

0.735 

S. Central 

1.200 


0.735 





(d) Fi 

nal step 



beams 

East 

N. Central 

West 

S. Central 

sub-total 

East 

3.380 

★ 

2.419 

1.2QQ 

7 

N. Central 

* 

2.310 

1.908 

1.190 

5.408 

West 

2.419 

1.908 

0.938 

0.735 

6 

S. Central 

1.200 

1.190 

0.735 

0.300 

3.425 
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comparison of how the regional -beam designs increase the communications 
supply, the supply matrix without these regional -beam designs (as shown 
in Table 5.10) is re-listed in Table 5.15, and the supply matrix with 
regional -beam designs is given in Table 5.16. It is obvious that the 
design has accomplished its objective: increasing the communications 

supply (i.e., satellite availability) to the four regions while 
decreasing it for no user. 

One possible beam arrangement that could provide the number of 
beams listed in Table 5.12 is shown in Figure 5.5. One possible 
communications traffic distribution that could provide the supply matrix 
as shown in Table 5.14(d) is shown in Figure 5.6, and is derived as 
follows. First, one should saturate the capacity of the one-beam 
satellites with intra-regional traffic supply, as in satellites S 2 , S 7 , 
and Sj 2 * Then one should try to saturate the capacity of the two-beam 
satellites with inter-regional traffic supply, as in S 3 , S 4 , S 5 , Sg, S 3 , 
Sg, Sjo, and S^; because of the limited supply given in Table 5.14(d), 
S 5 , S 3 , and Sg are not fully used by inter-regional traffic, and the 
remaining capacity in each is assigned to intra-regional traffic. The 
traffic of the three-beam satellites are then distributed to fulfill the 
remaining traffic supply quota in Table 5.14(d). Also note that as 
shown in Figure 5.6, additional intra-regional beams are available for 
North Central and South Central. Comparison of the supply matrix in 
Table 5.16, which is implemented by Figure 5.6, with the requirements in 
Table 5.9 show that this situation can be viewed as more than 100 % 
satisfaction of the demand for these regions. 
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Table 5.15 



Supply matrix without regional 

beams 


beams 

East 

N. Central 

West S. 

Central 

CONUS 

East 

1.305 

* 

0.934 

0.463 


N. Central 

* 

0.560 

0.592 

0.289 

17.753 

West 

0.934 

0.592 

0.291 

0.228 


S. Central 

0.463 

0.289 

0.228 

0.073 





Table 

5.16 




Supply matrix with regional 

beams 


beams 

East 

N. Central 

West S. 

Central 

rest 

East 

3.380 

* 

2.419 

1.200 


N. Central 

* 

2.310 

1.908 

1.190 

18 

West 

2.419 

1.908 

0.938 

0.735 


S. Central 

1.200 

1.190 

0.735 

0.300 
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Figure 5.6. Possible communications distribution for four regions. 


4. Discussion of Nunerlcal Example 

The above example was formulated by observing the communi cations - 
demand density on a map and selecting service areas intuitively 
according to the following criteria: the communications demand within 

and between the service areas should be high, the service areas should 
be small, and the distance between these service areas should be large. 
To some extent these criteria conflict, and a compromise was necessary. 
Although the service areas were specified intuitively, there should be 
ways to automate this process. The ultimate criterion for the choice of 
service area specification Is that it should result in the maximum 
communications supply. Intermediate criteria in terms of the AS and 
demand matrix elements would be useful if they can be developed. 

It was fortunate that the required spacing between the East and 
North Central satellites turned out slightly less than half the value 
for CONUS-beam satellites; this made the satellite “conversion" from 
CONUS to regional-beam particularly simple and straightforward. For 
other choices of service areas, the AS matrix will be different, so will 
the beam arrangement and the resulting communications supply. Still the 
principle is the same: smaller AS values means more satellites. It is 

hoped that this example will provide insight for generating more general 
and automated procedures. 

The concept of replacing administration-coverage satellites with 
regional -beam ones may have great practical merit. At present, 
satellite planning work may take place many years before the actual 
satellite implantation. It is very difficult to make a plan that can 
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cope with the technologies available twenty years later; therefore, a 
good plan should be flexible so that the administrations will not be 
tied to old technologies. For example, if CONUS-beam designs are used 
in the U.S. planning process with the option of "converting" some to 
multiple-beam designs later, higher information capacity can be obtained 
within the framework of an existing scenario, i.e., without changes in 
the remaining satellites. Therefore, from the viewpoint of flexibility, 
the multiple-beam replacement option should be attractive as a component 
of the U.S. planning. 

F. DISCUSSION AND CONCLUSION 

It was suggested that service areas of the FSS system should be 
specified according to the communications traffic-demand density in 
conjunction with the AS concept, because this could enable the system 
planner to specify more satellites and provide more communi cations 
supply. When applying this concept to specify satellites for an 
administration which has heavy traffic between several small and 
separated regions, it was shown that a mixture of administration 
coverage and intra-administration regional coverage can increase the 
communications capacity compared with only administration coverage. A 
numerical example was used to illustrate the design procedure for 
replacing several of a series of uniformly spaced administration- 
coverage satellites with regional -beam ones. It was shown that a 
substantial communications capacity increase could be obtained for many 
users without decreasing the capacity for any user. The procedure was 


171 



intuitive In part, and it Is recommended that techniques for formalizing 
and automating it be investigated. 

The regional service-area concept in this chapter Is closely 
related to the AS concept In Chapter IV. In Chapter IV, the 
specification of service areas Is given and fixed and the focus of the 
study is allocating satellite locations in the orbit. The focus in 
Chapter V is the specification of service areas. It is found that 
consideration of AS and traffic demand leads to useful techniques for 
this specification. 

An important advantage in specifying smaller (i.e., regional as 
opposed to administration-coverage) service areas in the case of 
geographically large administrations is the reduction of interference to 
other communications systems. A smaller service area means a smaller 
beam, and hence a faster drop of the field strength away from the 
service area. Thus, the use of regional service areas can give better 
interference protection to other service areas, not only within but also 
outside the same administration. 

In this chapter ,the inter-satellite service (ISS) was not 
considered. The adoption of the ISS might affect this study 
significantly. 
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CHAPTER VI 


CONCLUSION 

The purpose of this study has been to explore orbit and frequency 
assignment methods with which the spectrum/orbit resources can be 
efficiently utilized for satellite communications. 

In Chapter III, the mathematical nature of the orbital /frequency 
assignment problem is investigated by analyzing an objective function 
used by the extended-gradient and cyclic-coordinate search techniques. 

It is shown that the permutation of the orbital /frequency assignments is 
an important part of the problem. This indicates that a necessary 
condition for a technique to be able to find the globally optimal 
scenario is that it should to be able to deal with both the continuous 
aspects of satisfying the signal quality requirement and with the 
permutation problem. It is also shown that, at least when the frequency 
variables are fixed, for a given orbital permutation this objective 
function is likely to have only one local minimum. This suggests that a 
sufficient condition to obtain the globally optimal solution by a 
extended gradient search method is that it should terminate as the 
ordinary gradient search procedure with the optimal permutation in 
orbital and frequency assignments. 

In Chapter IV, a technique for obtaining the optimal orbital 
assignments is presented. The idea is to convert the signal quality 
requirement to minimum satellite spacing requirements and use them as 
constraints on the relative satellite locations. With a set of these 
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constraints, the assignment problem can be formulated as a mixed integer 
program (integer programming is often used to solve a permutation 
problem) and solved by the branch-and-bound method. The globally 
optimal solution is guaranteed; however, the required computational time 
may be prohibitively long when the problem size (number of satellites to 
be assigned) is large. To overcome this difficulty, the same problem is 
formulated as a linear program and solved by a version of the simplex 
method with restricted basis entry; only a locally optimal solution is 
guaranteed, but the method is more practical in terms of computational 
effort. The solutions of both programs are guaranteed to have 
satisfactory single-entry C/I e (carrier-to-effective interference ratio) 
results. 

In Chapter V, a basis of specifying service areas for the FSS 
system is proposed. It is suggested that while some satellites should 
cover large territories, some satellites should cover smaller regions 
where the communications demands are high. Smaller and separated 
service areas require smaller necessary satellite spacings, thus more 
satellites can be allocated in the orbit and higher communications 
supply can be obtained. The method involves simultaneous consideration 
of a requirement matrix and the AS matrix for the proposed regions. A 
numerical example, in which projected voice telephone demand is used for 
the requirement matrix, demonstrates the validity of this approach. 
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APPENDIX A 


STREAMLINED SOUP CODE 


C*** » MAIN PROGRAM « 

C THIS IS A MINI SOUP PROGRAM 

(^**«* ******************************************************** ******* 

c 

IMPLICIT INTEGERS(I-N) , REALS <A-H, O-Z) 

C 

CHARACTERS NAMESA 


C 


C 

C 


C 


C 


C 


C 


COMMON /CONSTS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10 , COIMIN, OOISING 

COMMON /PAR AMS/ NOMSAR, NAMESAUO ,2) ,NTPSA(10> ,CPHI0 (10) 


COMMON /VECTOR/ 

1 

2 


DSLON(IO) ,RSLON (10) ,XO<10> ,YO(10) , 
ROAIKJ t 

XOAC(IO) ,YOAC(10) , ZOAC (10) ,ROAC(10) 


COMMON /VARBLS/ 

1 


FREQ ( 10) r I POLAR ( 10) ,GAINR(10) ,GAINT(10) , 
EIRr(IO) , IFTNST(IG) r I PTN ER (10) 


COMMON /MIN ELL/ 

1 

1 


BCLAT (10) f BCLON (10) ,DBCLAT (10) ,DBCLON(10) , 
REFLAT (10) ,REFLON (10) ,AXR<10) , 

ORIENT (10) , AXMAJ ( 10) 


COMMON /TPOINT/ RELON ( 10 ,20) ,RELAT (10 , 20) , DELON ( 10 , 20) , 

1 DEL AT (10/20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 

COMMON /ANGLES/ YPH IT, YPHIR, PHITK, YPHIO 


C 

COMMON /REAL/ PIKJ f PKKJ , WFQM, WDRCK, YFWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

C 

q********* *****************************:. >**************************** 

C 

OPEN (UNIT- 20, FILE*' INPUTO.DAT' , TYPE-' OLD 1 ) 

OPEN (UNIT-6, FILE-'NASAP. DAT* ,TYPE- 1 NEW ' ) 

C 

CALL ICON ST 
CALL INPUTO 
C 

C ASSUME POWER DENSITIES AT ALL BEAM AIM POINTS ARE CONSTANT 
C 

DO 2 K - 1 ,NUMSAR 
CALL GAINER (K) 

EIRP(K) - PFD+10 • *DLOG10 ( 4 • *PI*ROAC (R) *ROAC (K) ) -ERDB 
2 CONTINUE 
C 

CALL ZFUNCT 


C 


C 

C 

C 

C 

C 

C 

C 


*** 

*** 


STOP 

END 

SUBROUTINE ZFUNCT 

>> THIS ROUTINE IS THE OVERALL CONTROL ROUTINE AT EACH 
>> STEP IN THE LINE SEARCH PROCEDURE 


IMPLICIT INTEGERS (I-N) , REALS (A-H, O-Z) 


« 

<< 
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nnn 


c 

CHARACTERS NAMES A 
C 

COMMON /CON STS/ 

1 
C 

COMMON /PAR AMS/ 

C 

COMMON /VECTOR/ 

1 
2 
C 

COMMON /VARBLS/ 

1 
C 

COMMON /MINELL/ 

1 
1 
C 

COMMON /TPOINT/ 

1 
C 

COMMON /ANGLES/ 

C 

COMMON /REAL/ 

1 
C 

C**** ******************************************************* ********* 

*** » INITIALIZE PARAMETER 

DO 10 1-1 ,NUMSAR 
CALL REFCAL ( I ) 

10 CONTINUE 
C 

C*** » OUTER SUMMATION (OVER K) FOR ALL SERVICE AREAS « 

C 

DO 1000 K «* 1 , NUMSAR 
JNTP - NTPSA(K) 

WRITES, 705) NAMESA(K,1) ,DSLON(K) 

C 

C*** >> MIDDLE SUMMATION (OVER J) FOR ALL TEST POINTS IN AREA K « 

C 

DO 900 J « 1 , JNTP 
CALL KPHI (K, J) 

CALL XPHI0 (K, K, J) 

CALL KPWDRC (K,J) 

CALL XPWFQ (FREQ (K) , FREQ (K) ) 

CALL XPOWER (K,K,J) 

C 

C*** » CALCULATE INTERFERENCE POWER 

C 

SUMP*0 . 

DO 500 I - 1 ,NUMSAR 

IF (I.EQ.K) GO TO 500 
CALL ZPHI (I , K, J) 

CALL XPHI0 (I , K, J) 

IF (IPOLAR(I) .EQ. IPOLAR(K) ) THEN 
CALL ZWDRC (I,K,J) 

ELSE 


E, PI, RADIAN, DEGREE, GCR, ER, ERE®, EAP, 

PFD, ALOG E, ALN1 0 , COIMIN , CO IS ING 

NUMSAR, NAMES A (10,2) ,NTPSA(10) ,CPHI0(10) 

DSL ON (10) ,RSLON (10) ,XO(10) ,YO(10) , 

ROAIKJ, 

XOAC(IO) ,YOAC (10) , ZOAC(IO) ,ROAC(10) 

FREQ (10) ,IPOLAR(10) ,GAINR(10) ,GAINT(10) , 
EIRP(10) , IPTNST (10) ,IPTNER(10) 

BCLAT(IO) ,BCLON (10) ,I»CLAT(10) ,DBCLON(10) , 
RE FLAT (10) ,REFLON(10) ,AXR(10) , 

ORIENT (10) ,AXMAJ (10) 

RELGN (10 ,20) ,RELAT(10 ,20) , DELON (10,20) , 
DELAT (10 ,20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 

YPH IT, YPH IR, PH ITK, YPH 10 

PIKJ , PKKJ , PWFQM, PWDRCK, YPWDRC, YPWDRX, 
XOAKKJ , YOAKKJ , ZOAKKJ , ROARKJ 
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CALL ZPWDRX <I,K,J) 

END IF 

CALL XFWFQ (FREQ (K) ,FREQ C I ) ) 

CALL XPOWER (I,K,J) 

COI * PKKJ-PIKJ 
DMGN * COI - COISING 
C 

WRITE (6 ,706) DELON (K,J) , DEL AT (K r J) , NAMES A (1 ,1) ,PIKJ , COI, DMGN 
C 

IF (PIKJ .LT. -700.) PIKJ « -700. 

C 

TEN1 - 10. ** (PIKJ/10.) 

SUMP • SUMP + TEN1 

500 CONTINUE 

C 

IF (NUMSAR. NE.2) THEN 
SUMPDB « 10 • *DLOG10 (SUMP) 

CO IT « PKKJ- SUMPDB 
TDMGN « CO IT- CO I MIN 

WRITE (6 ,708) DELON <K,J) , DEL AT (K, J) , SUMPDB, CO IT , TDMGN 
END IF 
C 

900 CONTINUE 
1000 CONTINUE 
C 

705 FORMAT (//,12X, 'TEST COUNTRY s ! ,A6 ,5X, 8 SATELLITE : 8 ,F8 .2 ,// ,15X, 
1*TEST POINT' ,10X,*INT. SAT. ■ ,2X, * INT. PWR* ,3X,«C/I (dB) 1 , 

25x •MARGIN*) 

706 FORMAT(15X,F7.2,2X,F7.2,5X,A6,4X, F8.2,4X, F6.2,6X, F6.2) 

708 FORMAT (13X, 2 ( 2X, F7 . 2) ,6X, •TOTAL* ,4X, F8 .2 ,4X, F6 .2 ,6X, F6 .2 ,/) 

RETURN 

END 

C 

SUBROUTINE ICON ST 
C 

C*** >> THIS ROUTINE INPUTS CONSTANTS THAT ARE USED IN THE PROGRAM <<; 

C 

c ************************** ***************************************** 

c 

IMPLICIT INTEGERM(I-N) ,REAL*8 (A-H, O-Z) 

C 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10 ,COIMIN, COISING 

C 

Q *+ ************************* ***************************************** 

C 

E * 2.7182818285 
PI = 3 .1415926536 
RADIAN * PI / 180.0 
DEGREE « 180.0 / PI 
GCR - 6.6134 
ER - 6.371E+06 
ALOGE - 0.4342944819 
ALN10 - 2.3025850930 
ERDB = -20.0 * DLOGIO(ER) 

PFD = -90. 

EAP * 0.6 
C 

RETURN 

END 
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c 

c 


SUBROUTINE INPUTO 


c*** 

» 

THIS ROUTINE INPUTS 

<< 


» 

- SERVICE AREA DATA 

« 

c*** 

» 

- INITIAL SATELLITE LOCATIONS 

<< 

c*** 

>> 

- IN ITIAL FREQUENCY ASSIGNMENTS 

<< 

c*** 

» 

- INITIAL POLARISATIONS 

<< 

c*** 

>> 

- ANTENNA PATTERN TYPE 



C******************************************************************* 

C 

IMPLICIT INTEGER*4 (I-N) ,REAL*8(A-H,0-Z) 

C 

C 


C 

C 

C 

C 


CHARACTER *6 NAMES A 

COMMON /CONSTS/ E , PI , RADIAN, DEGREE, GCR, ER, ERE®, EAP, 

1 PFD, ALOGE, ALN10 ,COIMIN, GOISING 

COMMON /PAR AMS/ NUMSAR, NAMESA(10 ,2) ,NTPSA(10) ,CPHI0(10) 

COMMON /VECTOR/ DSLONdO) ,RSL0N(10) ,XO(10) ,YO(10) , 

1 ROAIKJ, 

2 XOAC(IO) ,YOAC(10) ,ZOAC(10) ,ROAC(IO) 

COMMON /VARBLS/ FREQdO) ,IPOLAR(10) ,GAINR(10) ,GAINT(10) , 

1 EIRP(IO) , IPTNSTdO) , IPTNER(IO) 

COMMON /MINELL/ BCLAT(IO) ,BCLON(10) ,DBCLAT(10) ,DBCLON(10) , 
1 REFLAT (10) ,REFLGN(10) ,AXR(10) , 

1 ORIENT (10) ,AXMAJ( 10) 

COMMON /TPOINT/ RELCN (10 ,20) ,RELAT(10,20) ,DEL(W(10,20) , 

1 DEL AT (10 ,20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 


******************************************************************** 


*** » INPUT DATA FOR SATELLITE, BEAM, EARTH STATION CHARACTERISTICS 

READ (20,*) COISING 

CO I MIN - COISING - 5. 

READ (20,*) NUMSAR 
DO 10 N “ 1, NUMSAR 

READ(20 ,801) ( NAMES A (N, I) ,1b 1,2) 

READ (20 ,*)DSLON (N) ,FREQ(N) ,IPOLAR(N) ,IPTNST(N) , IPTNER(N) 

READ (20 ,*) DBCL(W(N) , DBCLAT(N) , ORIENT (N) ,AXMAJ(N) ,AXMIN 
READ (20,*) NTP 
DO 12 N1 =1 , NTP 

READ(20,*)DELCW(N,N1) ,DELAT (N, Nl) 

RELON (N, N1 ) - DELON (N,N1) ‘RADIAN 
RELAT (N, Nl) - DEL AT (N, Nl) ‘RADIAN 
12 CONTINUE 
C 

NTPSA(N) ■ NTP 

FREQ(N) » FREQ (N) * 1000. 

RSLON(N) - DSL CM (N) * RADIAN 
AXR(N) - AXMAJ ( N) /AXMIN 
ORIENT (N) - OR IENT(N) ‘RADIAN 
AXMAJ (N) - AXMAJ (N)*RADIAN 
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noon on non no non 


BCLON(N) - DBCLON(N) ‘RADIAN 
BCLAT(N) - DBCLAT (N) ‘RADIAN 
C 

10 CONTINUE 
C 

C*“ >> REFLECT INPUT DATA « 

C 

WRITE (6, 901) 

DO 20 N ■= 1 , NUMSAR 

WRITE (6 ,90 2) NAMESA(N,1) ,DSL0N(N) , 

$ FRBQ(N) , I POLAR (N) ,IPTNST(N) , IPTNER(N) 

20 CONTINUE 
C 

801 FORMAT ( 2 A6) 

901 FORMAT (//12X, 'COUNTRY' ,5X, ' SATELLITE' ,5X, ' FREQUENCY' ,3X, 

$ 'POLAR. ' ,3X, 'PTN(ST) ' ,3X, • FTN (ER) ',/) 

902 FORMAT (13X, A6 ,5X,F9 .2 ,5X, F9 .2 ,5X, 12 ,7X, 12 ,8X, 12) 

*“ >> SET UP THE VECTORS FOR TOE TEST POINTS « 

DO 3 0 K - 1, NUMSAR 

XOI - DCOS (RSLON (K) ) 

YOI - DSIN (RSLON (K) ) 

XO(K) ■ XOI 
YO (K) - YOI 

BCLT *= BCLAT(K) 

BCLN = BCLON(K) 

COSBLT = DCOS (BCLT) 

CXE = COSBLT*DCOS ( BCLN) 

CYE - COSBLT*DSIN (BCLN) 

CZE - DSIN ( BCLT) 

XOACI - CXE— GCR*XOI 
YOACI - CYE-GCR*YOI 
ZOAC(K) = CZE 

ROACI = DSQRT (XOACI *XOACI + YOACI ‘YOACI + CZE*CZE) 

XOAC(K) = XOACI 
YOAC(K) « YOACI 
ROAC(K) * ROACI 

DO 30 J = 1 » NTPSA (K) 

COSLAT = DCOS(RELAT(K,J)) 

XE(K,J) «= COSLAT * DCOS ( RELON (K, J) ) 

YE(K,J) - COSLAT * DSIN ( RELON (K, J) ) 

ZE (K, J) - DS IN ( REL AT (K , J) ) 

30 CONTINUE 

RETURN 
END 

SUBROUTINE GAINER(K) 

THIS IS TO CALCULATE EARTH RECEIVER GAIN FROM ANTENNA 
DIAMETER, AND TRANSMITTER GAIN FROM HPBW. 

IMPLICIT INTEGER*4 (I— N) ,REAL*8 (A-H, O-Z) 

C 
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CHARACTER *6 NAMESA 


COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10,COIMIN, COISING 

COMMON /VARBLS/ FREQ (10) ,IPOLAR(10) ,GAINR( 10) ,GAINT(10) , 

1 EIRP(10) ,IPTNST(10) ,IPTNER(10) 

COMMON /PAR AMS/ NUMSAR, NAMESA (10 ,2) ,NTPSA(10) , CPHIO(IO) 

COMMON /MINELL/ BCLAT(IO) ,BCLON(10) , DBCLAT(IO) ,DBCLON(10) , 
1 REFLAT(IO) ,REFLON(10) ,AXR(10) , 

1 ORIENT (10) ,AXMAJ( 10) 


XO - 223 ./180. 

GOTO (10,20,30,40 ,50) ,IPTNER(K) 

FSS, DIAMETER 3 METERS 
D « 3. 

WAVEL ■ 300. / FREQ (K) 

XI ■ D/ WAVEL 

GAINR(K) - 10 . *DLOG10 (PI *PI *X1*X1*EAP) 
CPHIO(K) - X0/X1 
GO TO 60 

FSS, DIAMETER 3 METERS 
D ■ 3. 

WAVEL « 300. / FREQ(K) 

XI - D/WAVEL 

GAINR(K) ■ 10.*DLOG10(PI*PI*X1*X1*EAP) 
CPHio(K) * X0/X1 
GO TO 60 

FSS, DIAMETER 4.5 METERS 
D •= 4.5 

WAVEL - 300. / FREQ(K) 

XI - D/WAVEL 

GAINR(K) - 10.*DLOG10(PI*PI*X1*X1*EAP) 
CPHIO(K) « X0/X1 
GO TO 60 

FSS, DIAMETER 4.5 METERS 
D *= 4.5 

WAVEL - 300. / FREQ (K) 

XI - D/WAVEL 

GAINR(K) * 10 . *DLOG10 (PI*PI*X1*X1*EAP) 
CPHIO(K) - X0/X1 
GO TO 60 

BSS, DIAMETER 1 METER 
D ■ 1. 

WAVEL - 300. / FREQ(K) 

XI « D/WAVEL 

GAINR(K) - 10.*DLOG10(PI*PI*X1*X1*EAP) 
CPHIO(K) - X0/X1 



<2° n n n non non on nonn nnnnnn n n nnn 


CALCULATE TRANSMIT GAIN FROM HPBW 

60 GAINT(K) - 10 . ‘DLOGIO (EAP*AXR (K) * (PI ‘XO/AXMAJ (K) ) “2) 

RETURN 

END 

SUBROUTINE KPHI (K,J) 

“* » THIS ROUTINE COMPUTES ME VECTOR COMPONENTS OF THE LINK << 

“* » K-K-J , AND THE CORRESPONDING TRANSMITTING PHI ANGLE << 

******************************************************************* 


IMPLICIT INTEGERS (I-N) , REAL *8 (A-H, O-Z) 

CHARACTER *6 NAMES A 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFDf ALOGE, ALN1 0 , CO I MIN, COISING 

COMMON /PAR AMS/ NUMSAR, NAMESA (10 ,2) ,NTPSA(10) ,CPH 10 ( 10) 

COMMON /VECTOR/ DSLCN (10) r RSLON (10) ,XO (10) ,YO (10) , 

1 ROAIKJ, 

2 XOAC(IO) ,YOAC(10) ,ZOAC(10) ,ROAC(10) 

COMMON /VARBLS/ FREQ (10) ,IPOLAR(10) ,GAINR(10) ,GAINT(10) , 

1 EIRP(10) rlPTNST(lO) , IPTNER(IO) 

COMMON /MIN ELL/ BCXAT(IO) ,BCLON(10) , DBCLAT(IO) ,DBCLON(10) , 
1 REFLAT (10) ,REFLON(10) f AXR(10) , 

1 ORIENT (10) ,AXMAJ (10) 


COMMON /T POINT/ BllCNi 


COMMON /ANGLES/ YPHIT,YPHIR, PHITK, YPHI0 


COMMON /REAL/ PI KJ , PKKJ , PWFQM, IWDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 


****** ******************************************************* ******* 


XOAKKJ - XE(K,J) - GCR * XO(K) 

YOAKKJ « YE(K,J) - GCR * YO(K) 

ZOAKKJ - ZE(K,J) 

ROAKKJ - DSQRT (XOAKKJ ‘XOAKKJ + YOAKKJ ‘YOAKKJ + ZOAKKJ ‘ZOAKKJ) 

COSPHI « (XOAC(K) * XOAKKJ + YOAC(K) * YOAKKJ + ZOAC(K) * ZOAKKJ) 
$ / (ROAC(K) * ROAKKJ) 

PHITK - DAR COS (COS PH I) 

RETURN 

END 

SUBROUTINE ZPHI (I,K,J) 

* >> THIS ROUTINE COMPUTES THE PHI ANGLES « 
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C*** » (FOR I INTERFERING WITH R TEST POINT J) « 

C 

C ************* ****************** ************************************ 

c 

IMPLICIT INTEGERM (I-N) ,REAL*8(A-H,0-Z) 

C 

CHARACTER*6 NAMESA 
C 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10 ,OOIMIN, OOISING 

C 

COMMON /PAR AMS/ NUMSAR, NAMESAdO ,2) ,NTPSA(10) ,CPHI0 (10) 

C 

COMMON /VECTOR/ DSLON(IO) ,RSL0N(10) ,XO(10) ,YO(10) , 

1 ROAIKJ, 

2 XOAC(IO) ,YOAC (10) ,ZOAC(10) ,ROAC (10) 

C 

COMMON /VARBLS/ FREQ(IO) , IPOLAR(IO) , GAINR(IO) , GAINT(IO) , 

1 EIRP(IO) ,IPTNST(10) ,IPTNER(10) 

C 

COMMON /MIN ELL/ BCLAT(IO) , BCLON(IO) ,DBCLAT(10) , DBCLON(IO) , 

1 REFLAT(IO) ,REFLON(10) ,AXR(10) , 

1 ORIENT ( 10) ,AXMAJ( 10) 

C 

COMMON /TPOINT/ RELON ( 10, 20 ) ,RELAT( 10,20) , DELON (10,20 ) , 

1 DELAT (10 ,20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 

C 

COMMON /ANGLES/ YPHIT, YPHIR, PHITK, YPH 10 
C 

COMMON /REAL/ PI KJ , PKKJ , WFQM, IWDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

C 

£*•**«** ***************************. ********************************* 

C 

C*** » CALCULATE OFF AXIS VECTOR COMPONENTS ( IKJ) << 

C 

XOAIKJ - XE (K, J) - GCR * XO(I) 

YOAIKJ - YE (K, J) - GCR * YO(I) 

ZOAIKJ = ZE(K,J) 

ROAIKJ «= DSQRT (XOAIKJ *XOAIKJ + YOAIKJ *YOAIRJ + ZOAIKJ *ZOAIKJ) 

C 

C*** » COMPUTE DISCRIMINATION ANGLES « 

C*** » FOR THE TRANSMITTING ANTENNA « 

C 

TNUMER = XOAC(I) * XOAIKJ + YOAC(I) * YOAIKJ + ZOAC(I) * ZOAIKJ 
TDENOM = ROAC(I) * ROAIKJ 
TEMPU - TNUMER / TDENOM 
C 

YPHIT = 0 .0 

IF (DABS(TEMPU) .LT. 1.0) YPHIT » DARCOS (TEMPU) 

C 

C*** » COMPUTE DISCRIMINATION ANGLES « 

c *** >> poR THE RECEIVING ANTENNA « 

c 

TNUMER - XOAKKJ * XOAIKJ + YOAKKJ * YOAIKJ + ZOAKKJ * ZOAIKJ 
TDENOM = ROAKKJ * ROAIKJ 
TEMPU * TNUMER / TDENOM 
C 

YPHIR *0.0 

IF (DABS (TEMPU) .LT. 1.0) YPHIR « DARCOS (TEMPU) 
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c 

RETURN 

END 

C 

SUBROUTINE XPHIO (I f K,J) 

C 

C*** » THIS ROUTINE COMPUTES TOE ELLIPTICAL BEAM HALF POWER 
C*** >> BEAM WIDTH USING THE METHOD GIVEN IN TOE SOUP-3 MANUAL 
C 

£**** ******************************************************* ******** 


C 

c 

c 


c 

c 


c 


c 


c 


c 

c 


c 


IMPLICIT INTEGER *4 ( I-N) , REAL *8 (A-H,0-Z) 
CHARACTERS NAMES A 


COMMON /CONSTS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10 , COIMIN, OOISING 

COMMON /PAR AMS/ NUMSAR, NAMESA<10,2) ,NTPSA(10) ,CPHI0 (10) 


COMMON /VECTOR/ 

1 

2 


DSL ON (10) ,RSLON(10) ,XO(10) ,YO(10) , 
ROAIKJ, 

XOAC(IO) ,YOAC(10) , ZOAC(IO) ,ROAC(10) 


COMMON /VARBLS/ FREQ(IO) , IPOLAR(IO) ,GAINR(10) , GAINT(IO) , 
1 EIRP(IO) ,IPTNST(10) ,IPTNER(10) 


COMMON /MIN ELL/ 

1 

1 


B CL AT (10) ,BCLON(10) ,DBCLAT(10) ,DBCLON(10) , 
REFLAT (10) ,REFLON(10) ,AXR(10) , 

ORIENT(IO) r AXMAJ (10) 


COMMON /TPOINT/ RELON (10 ,20 ) ,RELAT (10 ,20) , DELON ( 10 ,20) , 

1 DELAT (10,20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 

COMMON /ANGLES/ YPH IT, YPH IR, PH ITK, YPH 10 


COMMON /REAL/ PI KJ , PKKJ , PWFQM, WDRCK, YPWDRC, YPWDRX, 
1 XOAKRJ , YOAKKJ , ZOAKKJ , ROAKKJ 


£************************ ************************************* ******* 


C 

PHIS - RSLON(I) 

PHIC - BCLON(I) 

THETAC - BCLAT(I) 

C 

SINPHS - DSIN (PHIS) 

COSPHS - DCOS(PHIS) 

COSTC - DCOS (THETAC) 

C 

COSW *= COSTC * DCOS (PH IS- PHIC) 

SINK « DSIN (DAR COS (COSW) ) 

C 

RSLAT - 0. 

IF ( BCLON(I) .EQ.RSLON(I) ) GO TO 1 
IF ( BCLAT(I) .EQ.RSLAT ) GO TO 2 

ARG = OOSTC * DSIN (DABS ( BCL ON ( I ) -RSLON (I) ) ) / SINW 
A * DAS IN (ARG) 

IF (BCLON(I) .GT.RSLON(I) .AND. BCLAT (I ) .GT. RSLAT) A«2.*PI-A 
IF (BCLON(I) .LT.RSLON(I) .AND. BCLAT (I ) .LT. RSLAT) A«PI-A 
IF (BCLON(I) .GT.RSLON(I) .AND. BCLAT (I ) .LT. RSLAT) A-PI+A 
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COSA - DCOS (A) 

SINA - DSIN(A) 

GO TO 3 

1 CONTINUE 
COSA * -1 . 

SINA - 0 . 

A - PI 

IF (BCLAT(I) . LT. RSLAT) GO TO 3 
COSA - 1. 

A * 0 . 

GO TO 3 

2 CONTINUE 
COSA - 0. 

SINA - 1. 

A - PI/2. 

IF (BCLON(I) .LT.RSLCN(I)) GO TO 3 
SINA « -1 . 

A - 1 ,5*PI 

3 CONTINUE 
C 

TANT « SINW / (GCR - COSW) 

TRAD - DATAN (TANT) 

SINT - DS IN (TRAD) 

COST « DCOS (TRAD) 

C 

A21 - - COSA * SINFHS 

A22 « COSA * COS PH S 

A23 - SINA 

C 

A31 - SINT * COSPHS + COST * SINA * SINPHS 

A32 - SINT * SINPHS - COST * SINA * COSPHS 

A3 3 - COST * COSA 

C 

FLAT » REFLAT (I) 

FLON - REFLON (I) 

CSFLAT - DCOS (FLAT) 

VRI - CSFLAT * DCOS (FLON) 

VR2 - CSFLAT * DS IN (FLON) 

VR3 - DS IN (FLAT) 

C 

BLAT « BCLAT(I) 

BLON - BCLON(I) 

CSBLAT - DCOS (BLAT) 

VC1 « CSBLAT * DCOS ( BLON ) 

VC2 - CSBLAT * DS IN ( BLON ) 

VC3 - DSIN(BLAT) 

C 

PH IE « RELON (K, J) 

THETAE - RELAT (K» J) 

COSTE « DCOS (THETAE) 

VE1 - COSTE * DCOS (PH IE) 

VE2 - OOSTE * DS IN (PH IE) 

VE3 - DSIN (THETAE) 

C 

VS1 “ GCR * COSPHS 
VS2 - GCR * SINPHS 
VS3 - 0.0 
C 

VRMVC1 * VRI - VC1 
VRMVC2 - VR2 - VC2 
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VRMVC3 « VR3 - VC3 
C 

VEMVS1 - VEl - VS1 
VEMVS2 « VE2 - VS2 
VEMVS3 - VE3 
C 

S1NUMR - A3 1*VRMVC1 + A32*VRMVC2 + A33*VRMVC3 
SI DEN R - A21*VRMVC1 + A22*VRMVC2 + A23*VRMVC3 
C 

S2NUMR » A31*VEMVS1 + A32*VEMVS2 + A33*VEMVS3 
S2DENR « A21*VEMVS1 + A22*VEMVS2 + A23*VEMVS3 
C 

IF (S1DENR. NE .0.0) GO TO 10 

51 - PI / 2.0 
GO TO 15 

10 SI - DAT AN ( SI NO HR/ SI DEN R) 

C 

15 IF (S2DENR. NE .0.0) GO TO 20 

52 « PI / 2.0 
GO TO 25 

20 S2 « DATAN ( S2NUMR/ S2DENR) 

C 

25 SIGMA « S2 - SI 
CS - DCOS( SIGMA) 

SS - DS IN (SIGMA) 

AR - AXR(I) 

P0 » AXMAJ(I) / DSQRT (CS *CS + AR*AR * SS*SS) 

C 

ypHio « po 

RETURN 

END 

C 

SUBROUTINE REFCAL(N) 

*** >> THIS ROUTINE CALCULATES THE REFERENCE POINT LAT. & LON. << 
*** >> BASED ON THE ALGORITHM IN SOUP MANUAL 3.4, MAY 1983 << 

******************************************************************* 


c 

c 


c 

c 


c 


c 


c 


IMPLICIT INTEGERS (I-N) f REAL *8 (A-H, O-Z) 
CHARACTER *6 NAMES A 


COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, AL0GE, ALN1 0 , COIMIN, OOISING 

COMMON /PAR AMS/ NUMSAR, NAMESAdO ,2) ,NTPSA(10) ,CPH 10 ( 10) 


COMMON /VECTOR/ 

1 

2 


DSL ON (10) ,RSLON(10) ,XO(10) ,YO(10) , 
ROAIKJ, 

XOAC(IO) ,YOAC (10) , ZOAC (10) ,ROAC(10) 


COMMON /VARBLS/ FREQ(IO) ,IPOLAR(10) ,GAINR(10) ,GAINT(10) , 
1 EIRP(IO) ,IPTNST(10) rIPTNERdO) 


COMMON /MINELL/ 

1 

1 


BCLAT(IO) , BCLON (10) ,DBCLAT(10) ,DBCLON(10) , 
REFLAT (10) ,REFL(N(10) ,AXR(10) , 

ORIENT (10) , AXMAJ (10) 


COMMON /TPOINT/ RELCN (10 ,20) ,RELAT (10 ,20) , DELON (10 ,20) , 
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1 DELAT (10,20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 

C 

COMMON /ANGLES/ YPHIT, YPHIR, PHITK, YPHIO 
C 

COMMON /REAL/ PIKJ, PKKJ , FWFQM, IWDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

C 

C*********** ********************************************************* 

c 

PGMPS *= BCLON(N) - RSLON(N) 

COSTG - DCOS ( BCL AT ( N) ) 

COSPP - DCOS (PGMPS) 

Q1 « OOSTG * DSIN (PGMPS) 

Q2 - GCR - COSTG * COSPP 
Q3 - DSQRT (Ql**2+Q2**2) 

C 

SX > DATAN2 (Q1 ,Q2) 

SY - DATAN2(DSIN(BCLAT(N) ) ,Q3) 

C 

AMAJ2 - AXMAJ(N) * 0.5 

SX2 - AMAJ2 * DCOS (ORIENT (N) ) 

SY2 - AMAJ2 * DS IN (ORIENT (N) ) 

C 

XI - SX + SX2 
Yl - SY + SY2 
C 

Q4 - DAR COS (DCOS (XI) * DCOS (Yl)) 

Q5 - DS IN ( Q4 ) 

T « GCR * Q5 

IP (T .LE. 1.0) GO TO 10 
C 

XI ■ SX - SX2 
Yl * SY - SY2 

Q4 - DAR COS (DCOS (XI) * DCOS(Yl)) 


Q5 

- DS IN ( Q4 ) 


T » 

> GCR * Q5 


IF 

(T .LE. 1.0 

.AND. T . GE. -1.0) GO TO 10 

IF 

(T .GT. 0.) 

T « 1.0 

IF 

(T .LT. 0.) 

T - -1.0 


C 

WRITE (6, 901) T 

901 FORMAT (/10X, ******* POSSIBLE ERROR IN ELLIPSE SELECTION 
$ F5.1) 

C 

10 PX = DARSIN (DSIN (Yl) / Q5) 

IF (XI .LT. 0.0) PX « PI - PX 
D - DARSIN (T) 

BLAM ■ D - Q4 

REFLAT (N) ■ DARSIN (DSIN ( BLAM) * DSIN(PX)) 

AL - DAR COS (DCOS ( BLAM) /DCOS ( REFLAT (N) ) ) 

IF ( DABS(PX) .GT. PI/2.) AL - -AL 
REFLON(N) - RSLON(N) + AL 
C 

RETURN 

END 

C 

C 

FUNCTION DARSIN (X) 

C 

IMPLICIT REAL*8(A-H,0-Z) 
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DARSIN - DAS IN (X) 

RETURN 

END 

C 

FUNCTION DAROOS(X) 

C 

IMPLICIT REAL *8 (A-H f O-Z) 

DARCOS « DA COS (X) 

RETURN 

END 

C 

SUBROUTINE KPWDRC (K,J) 

C 

C*** » THIS ROUTINE COMPUTES 1HE CARRIER SIGNAL DISCRIMINATING << 
C*** >> FACTOR (DESIRED POWER - CO-POLARISED ONLY - PHIR«0 ) « 

C 

q******************************************************************* 

C 

IMPLICIT INTEGERM ( I-N) , REAL *8 (A- H, O-Z) 

C 

CHARACTER *6 NAMES A 
C 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10 , COIMIN, OOISING 

C 

COMMON /PAR AMS/ NUMSAR,NAMESA(10,2> ,NTPSA(10) ,CPHI0 (10) 

C 

COMMON /VECTOR/ DSLON(IO) ,RSLON(10) ,XO(10) ,YO(10) , 

1 ROAIKJ, 

2 XOAC(IO) ,YOAC (10) ,ZOAC(10) ,ROAC(10) 

C 

COMMON /VARBLS/ FREQ(IO) ,IPOLAR(10) ,GAINR(10) ,GAINT(10) , 

1 EIRP(IO) ,IPTNST(10) ,IPTNER(10) 

C 

COMMON /MIN ELL/ BCLAT QO) ,BCLON ( 10) ,DBCLAT (10) ,DBCLON ( 10) , 
i REFLAT(IG) ,REFLCN (10) ,AXR(iC) , 

1 ORIENT (10) ,AXMAJ( 10) 

C 

COMMON /TPOINT/ RELON (10 ,20) ,RELAT (10 ,20) , DELON ( 10 ,20) , 

1 DELAT (10 ,20) ,XE (10,20) ,YE(10,20) ,ZE(10,20) 

C 

COMMON /ANGLES/ YPHIT, YPHIR, PHITK, YPHIO 
C 

COMMON /REAL/ PIKJ , PKKJ , WFQM, PWDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

C 

£** ********************************************* ********************* 

c 

PT - PHITK 
PO *= YPHIO 
C 

GO T0(10,20,30,40,50) , IPTNST (K) 

C 

10 CALL PTNST1 (PT, PO ,GAINT (K) ,DISC) 

GOTO 60 
C 

20 CALL FTNST2 (PT, PO ,GAINT(K) ,DISC) 

GOTO 60 
C 

30 CALL PTNST3 (PT, PO ,GAINT (K) ,DISC) 
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GOTO 60 
C 

40 CALL PTNST4 ( PT , PO ,GAINT (K) , DISC) 

GOTO 60 
C 

50 CALL PTNST5 (PT, P0 ,GAINT (K) ,DISC) 

C 

60 YDSCT « DISC 
C 

C PREVENTING INPERFECT MINIMUM ELLIPSE NOT COVERING ALL POINTS 
C 

C IF (YDSCT.LT.-3.) YDSCT— 3 . 

C 

C*** » COMPONENT OF DESIRED POWER DEPENDENT ON ORBIT LOCATION « 

C 

PWDRCK « EIRP(K) + YDSCT + GAINR(K) - 20 .0 * DLOG10 (ROAKKJ) 
RETURN 
C 

END 

C 

SUBROUTINE ZIWDRC (I,K,J) 

C 


C*** 

>> 

THIS ROUTINE 

COMPUTES THE 

INTERFERENCE SIGNAL 

<< 

C*** 

>> 

TRANSMITTING 

AND 

RECEIVING 

DISCRIMINATIONS FOR 

<< 

C*** 

>> 

THE CO- POLAR IS ED 

CASE 


<< 


c 

£**** ******************************************************* ******** 

c 

IMPLICIT INTEGER*4 (I-N) ,REAL*8 (A-H, O-Z) 

C 

CHARACTER *6 NAMESA 


C 

C 

c 

c 


c 

c 


COMMON /CONSTS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD# ALOGE, ALN10 , COIMIN, COISING 

COMMON /PAR AMS/ NUMSAR, NAMESA (10 , 2) ,NTPSA(10) ,CPHI0 ( 10) 

COMMON /VECTOR/ DSLON ( 10) ,RSLON (10) , XO (10) r YO (10) f 

1 ROAIKJ, 

2 XOAC(IO) r YOAC(10) ,ZOAC(10) ,ROAC(10) 

COMMON /VARBLS/ FREQ (10) r IK)LAR(10) ,GAINR(10) ,GAINT(10) , 

1 EIRP(10) f IPTNST (10) ,IPTNER(10) 

COMMON /MINELL/ BCLAT(IO) #BCLON(10) f DBCLAT(10) , DBCLON(IO) , 
1 REFLAT (10) ,REFLON(10) ,AXR(10) , 

1 ORIENT (10) ,AXMAJ( 10) 

COMMON /TPOINT/ RELON (10 ,20) ,RELAT (10 ,20) ,DELON ( 10 , 20) f 
1 DELAT (10 f 20) ,XE(10,20) ,YE(10,20) r ZE(10 r 20) 

COMMON /ANGLES/ YPH IT, YPHIR, PHITK, YPHI0 


COMMON /REAL/ 


PIKJ , PKKJ , PWFQM, PWDRCKr YPWDRC, YPWDRX, 
XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 


C** ******* *********************************************************** 

O 

C*** » TRANSMITTING DISCRIMINATION << 

C 
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PT * YPHIT 
PO « YPHIO 

GO T0(10,20,30,40,50) ,IFTNST(I) 

C 

10 CALL PTNST1 (PT f PO ,GAINT (I ) , DISC) 

GOTO 60 
C 

20 CALL PTNST2(PT,P0,GAINT(I) ,DISC) 

GOTO 60 
C 

30 CALL FTNST3(PT,P0,GAINT(I) ,DISC) 

GO TO 60 
C 

40 CALL FTNST4 (PT, P0 ,GAINT (I ) ,DISC) 

GO TO 60 
C 

50 CALL FTNST5 (PT, P0 ,GAINT (I ) ,DISC) 

C 

60 YDSCT - DISC 
C 

C*** » RECEIVING DISCRIMINATION << 

C 

PR * YPHIR 
P0 - CPHIO(K) 

C 

GO TO (110 ,120 ,130,140 ,150) , IPTNER(K) 

C 

110 CALL PTNER1 (PR, P0 , FREQ ( I ) ,GAINR(K) ,DISC) 

GOTO 160 
C 

120 CALL PTNER2 (PR, P0 ,FREQ (I ) ,GAINR(K) ,DISC) 

GOTO 160 
C 

130 CALL PTNER3 (PR, P0 ,FREQ (I ) , GAINR(K) ,DISC) 

GO TO 160 

/-* 

140 CALL FTNER4(PR,P0,FREQ(I) ,GAINR(K) ,DISC) 

GO TO 160 
C 

150 CALL PTNER5(PR,P0,FREQ(I) ,GAINR(K) ,DISC) 

C 

160 YDSCR - DISC 
C 

C*** » COMPONENT OF POWER DEPENDENT ON TOE ORBIT LOCATION << 

*** » CO-POLARISED CASE « 

YPWDRC - EIRP(I) + YDSCT + YDSCR 
$ - 20. * DLOG1 0 ( ROAIKJ ) 


RETURN 

END 

SUBROUTINE ZIWDRX (I,K,J) 

*** » THIS ROUTINE COMPUTES THE INTERFERENCE POWER << 

*** >> TRANSMITTING AND RECEIVING DISCRIMINATIONS FOR « 

*** >> THE CROSS- POLARISED CASE « 

******************************************************************* 
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C 


IMPLICIT INTEGER *4 (I-N) ,REAL*8(A-H,0-Z) 
CHARACTER*6 NAMES A 


C 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD,ALOGE,ALN10,COIMIN,COISING 

C 

COMMON /PAR AMS/ NUMSAR, NAMESA (10 ,2) ,NTPSA (10) ,CPH 10 ( 10) 

C 

COMMON /VECTOR/ DSLCN(IO) ,RSLON(10) ,XO(10) ,YO(10) , 

1 ROAIKJ , 

2 XOAC(IO) ,YOAC (10) ,ZOAC(10) ,ROAC(10) 

C 

COMMON /VARBLS/ PRBQ(IO) ,IPOLAR(10) ,GAINR(10) ,GAINT(10) , 

1 EIRP(10) ,IPTNST(10) , IPTNER(IO) 

C 

COMMON /MIN ELL/ BCLAT (10) ,BCLON (10) ,DBCLAT (10) ,DBCLON ( 10) , 

1 REFLAT(IO) ,REFLON(10) ,AXR(10) , 

1 ORIENT (10) ,AXMAJ (10) 

C 

COMMON /TFOINT/ RELON(10,20) ,RELAT (10,20) ,DELON(10,20) , 

1 DELAT (10 ,20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 

C 

COMMON /ANGLES/ YPHIT, YPHIR, PH ITK, YPHIO 
C 

COMMON /REAL/ PIKJ, PKKJ, WPQM, WDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

******************************************************************** 


*** » TRANSMITTING DISCRIMINATION << 

PT = YPHIT 
PO ■ YPHIO 

GOTO (10,12,14,16,18) ,IPTNST(I) 

10 CALL PTNST1(PT,P0,GAINT(I) , DISC) 

GO TO 20 

12 CALL PTNST2(PT,P0,GAINT(I) , DISC) 

GO TO 20 

14 CALL FTNST3(PT,P0,GAINT(I) , DISC) 

GO TO 20 

16 CALL PTNST4(PT,P0,GAINT(I) , DISC) 

GO TO 20 

18 CALL PTNST5(PT,P0,GAINT(I) , DISC) 

20 YDSCT - DISC 

'** » RECEIVING DIRECTIVITY « 

PR - YPHIR 
PO - CPHIO(K) 

C 

CALL XPTNER1 (PR, PO ,GAINR(K) ,DISC) 
YDSXR «= DISC 
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c 

C*** >> COMPONENT OF POWER DEPENDENT ON TOE ORBIT LOCATION << 

C*** » CROSS- POLARISED CASE << 

C 

100 YPWDRX - EIRP(I) + YDSCT + YDSXR 

$ - 20,0 * DLOGIO (ROAIKJ) 

C 

RETURN 

END 

C 

SUBROUTINE XPWFQ <FQD,FQI) 

C 

C*** » THIS ROUTINE COMPUTES TOE FREQUENCY DEPENDENT PORTION « 
C*** » IN TOE POWER EQUATION « 

C 

C ********************** ********************************************* 

C 

IMPLICIT INTEGERS (I-N) , REAL*8 (A-H f O-Z) 

C 

CHARACTER *6 NAMESA 
C 

OOMMON /CONSTS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10,COIMIN,OOISING 

C 

COMMON /FAR AMS/ NCMSAR, NAMESA (10 ,2) ,NTPSA (10) ,CPH 10 ( 10) 

C 

COMMON /VECTOR/ DSLON ( 10) ,RSLON ( 10) , XO (10) ,YO (10) , 

1 ROAIKJ, 

2 XOAC(IO) ,YOAC(10) ,ZOAC(10) ,ROAC<10) 

C 

COMMON /VARBLS/ FREQ ( 10) , IFOLAR ( 10) ,GAINR( 10) ,GAINT (10) , 

1 EIRP(IO) ,IFTNST(10) ,IPTNER(10) 

C 

COMMON /MIN ELL/ BCLAT (10) ,BCLON ( 10) ,DBCLAT (10) ,DBCLON ( 10) , 

1 REFLAT(IO) ,REFLON(10) ,AXR(10) , 

1 ORIENT (10) ,AXMAJ(10) 

C 

COMMON /TPOINT/ RELON (10 ,20) ,RELAT (10 ,20) , DELON ( 10 ,20) , 

1 DELAT (10 ,20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 

C 

COMMON /ANGLES/ YPH IT, YPH IR, PH ITK, YPHI0 
C 

COMMON /REAL/ PIKJ , PKKJ , PWFQM, FWDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

C 

C* ******************************************************************* 

C 

C*** » POWER BEING CALCULATED IS INTERFERING POWER << 

C 

X = (FQI - FQD) 

ABSX = DABS(X) 

IF (ABSX. LE. 1 5 •) THEN 
FF * 0. 

ELSE 

FF = -(ABSX-15.) *1.6 
END IF 
C 

PWFQM *= FF - 20.0 * DLOGIO (FQI) 

C 
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RETURN 

END 

C 

SUBROUTINE XPOWER (I,K,J) 

*** » THIS ROUTINE CALCULATES ALL TOE POWERS AFTER THE << 

*** » DIRECTIVITY AND FREQUENCY PORTIONS ARE COMPUTED « 

******************************************************************* 

IMPLICIT INTEGER*4(I-N) ,REAL*8 (A-H, O-Z) 

CHARACTER*6 NAMESA 

COMMON /CONSTS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10 ,COIMIN, OOISING 

COMMON /PAR AMS/ NUMSAR,NAMESA(10,2) ,NTPSA(10) ,CPHI0(10) 

COMMON /VECTOR/ DSLON(IO) ,RSLON(10) ,XO(10) ,YO(10) , 

1 ROAIKJ, 

2 XOAC(IO) r YOAC (10) , ZOAC(IO) ,ROAC(10) 

COMMON /VARBLS/ FREQ(IO) ,IPOLAR(10) ,GAINR(10) ,6AINT(10) , 

1 EIRP(IO) , IPTNST(IO) , IFTNER(IO) 

COMMON /MINELL/ BCLAT(IO) ,BaON(10) ,I»CLAT(10) ,I»CLON(10) , 

1 REFLAT(IO) ,REFLON(10) , AXR(IO ) , 

1 ORIENT (10) ,AXMAJ( 10) 

COMMON /TPOINT/ RELCN ( 10 ,20) ,RELAT (10 ,20) , DELON (10 ,20) , 

1 DELAT (10,20) ,XE(10,20) ,YE(10,20) ,ZE(10,20) 

COMMON /ANGLES/ YPH IT, YPHIR, PHITK, YPHIO 

COMMON /REAL/ PIKJ, PKKJ , FWFQM, FWDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

******************************************************************** 


IF (K. BQ .1) GO TO 100 

<** » POWER BEING CALCULATED IS AN INTERFERING POWER « 

IF (IPOLAR(I). EQ . IPOLAR(K) ) GO TO 50 

** » CROSS- POLARIZED « 

PIKJ ■= YPWDRX + PWFQM + ERE© + 27.56 
RETURN 

** >> CO-POLARIZED « 

50 PIKJ - YPWDRC + PWFQM + ERDB + 27.56 
RETURN 

** » POWER BEING CALCULATED IS A DESIRED POWER « 

100 PKKJ « FWDRCK + IWFQM + ERE© + 27.56 
RETURN 
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c 

END 

C 

SUBROUTINE PTNST1 (PT, PO ,G, DISC) 

FSS SATELLITE TX PATTERN FROM CCIR REPORT 558-2 

IMPLICIT INTEGER *4 ( I-N) ,REAL*8(A-H,0-Z) 

COMMON /CONSTS/ E , PI , RADIAN , DEGREE, GCR, ER, ERDB, EAP, 
1 PFD, ALOGE, ALN10 , COIMIN, COISING 

X*PT/P0 

IF (X .LE. 1.3) GO TO 10 
IF (X .LE. 3.15) GO TO 20 

DISC * -(7.5 + 25.0 * DLOGIO (X) ) 

IF (DISC. LE. ( -G-10 . ) ) DISC* -G-10. 

GO TO 40 

10 DISC * -12.0 * X * X 
GO TO 40 

20 DISC « -20.0 

40 RETURN 
END 

SUBROUTINE FTNST2 (PT, P0 ,G, DISC) 

FSS SATELLITE TX PATTERN FROM RARC 83 $ 5.1.10.1 

IMPLICIT INTEGERMU-N) , REAL *8 (A-H,0-Z) 

C 

COMMON /CONSTS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 
1 PFD, ALCGE, ALN1G ,OGIM!N, COISING 

C 

DP0 =P0 ‘DEGREE 
XI =PT/P0 
X2 =DP0/ 0.8 
X3 *.5* (l.-l ./X2) 

C 

PI * .4/DP0+X3 
P2 =1 .155/DPO +X3 
P3 *1 .60/DP0+X3 
P4 = 4.0/DP0+X3 
P5 = 6.968/DPO +X3 


P6 

* 10 

.** ( 

(G-ll 

.5) /25, 

.)/X2 

IP 

(XI 

• LE. 

0.5) 

GO TO 10 

IF 

(XI 

• LE. 

P2) 

GO 

TO 

12 

IF 

(XI 

.LE. 

P3) 

GO 

TO 

14 

IF 

(XI 

• LE. 

P4) 

GO 

TO 

16 

IF 

(XI 

. LE. 

P5) 

GO 

TO 

18 

IF 

(XI 

.LE. 

P6) 

GO 

TO 

20 


C 

DISC = -G 
GO TO 40 
C 

10 DISC = -12.0 * XI * XI 
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GO TO 40 


C 

12 

C 

14 

C 

16 

C 

18 

C 

20 

C 

40 

C 


c 


c 


c 


c 


c 

10 

c 

12 

C 

14 

C 

16 

C 

40 

C 


DISC - -18.75*DP0*DP0*(X1-X3)*(X1-X3) 

GO TO 40 

DISC - -25 
GO TO 40 

DISC » -17.5 - 25 • *DLOG10 ( (X1-X3) *X2) 

GO TO 40 

DISC * -35. 

GO TO 40 

DISC- -11.5 - 25.* DLOG10( (X1-X3) *X2) 

RETURN 

END 

SUBROUTINE FTNST3 (PT, P0 ,G, DISC) 

PSS SATELLITE TX PATTERN PROM RARC 83 $ 5.1.10.1 
WITH MODIFICATION 

IMPLICIT INTEGERS (I-N) ,REAL*8 (A-H f O-Z) 

COMMON /CONST S/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 
1 PFD, ALOGE, ALN1 0 , COIMIN, COISING 

DP0 « P0*DEGREE 
XI - PT/P0 
X2 « DP0/0.8 
X3 - .5* (1 .-1 «/X2) 

IF (XI .LE. 0.5) GO TO 10 
P2 - 1.265/DPO +X3 
IF (XI .LE. P2) GO TO 12 
P3 - 10 • **( (30.-24 .) /30 .) 

IF (XI .LE. P3) GO TO 14 
P4 - 10 . *M (G-24 . ) /30.) 

IF (XI .LE. P4) GO TO 16 

DISC « -G 
GO TO 40 

DISC « -12.0 * XI * XI 
GO TO 40 

DISC - -18 .75*DP0*DP0* (X1-X3) * (X1-X3) 

GO TO 40 

DISC = -30. 

GO TO 40 

DISC- -24.-30. *DLOG10(X1) 

RETURN 

END 

SUBROUTINE ETNST4 (PT, P0 ,G, DISC) 


194 



noon 


c 

C PSS SATELLITE TX PATTERN FROM RARC 83 $5.10.1 
C WITH MODIFICATION 
C 

IMPLICIT INTEGER *4 ( I-N) ,REAL*8(A-H,0-Z) 

C 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 
1 PFD, ALOGE, ALN10 , COIMIN, OOISING 

C 

DPO - PO *DEG REE 
XI - PT/PO 
X2 - DPO/ 0 .8 
X3 « .5M1.-1./X2) 

C 

IF (XI .LE. 0.5) GO TO 10 
P2 - 1 .26 5/ DPO +X3 
IF (XI .LE. P2) GO TO 12 
P3 ■ 10 . ** ( (30.-24 .) /30.) 

IF (XI .LE. P3) GO TO 14 
P4 «= 10.**((G-24.)/30.) 

IF (XI .LE. P4) GO TO 16 
C 

DISC - -G 
GO TO 40 
C 

10 DISC = -12.0 * XI * XI 
GO TO 40 
C 

12 DISC - -18 .7 5*DP0*DP0* (X1-X3) * (X1-X3) 

GO TO 40 
C 

14 DISC = -30. 

GO TO 40 
C 

16 DISC = -24.-30. *DLOG10 (XI) 

C 

40 RETURN 
END 
C 

SUBROUTINE PTNST5 (PT, PO , G, DISC) 

SATELLITE TX PATTERN FROM RARC 83 P.111, 

BSS PATTERN 

IMPLICIT INTEGER*4 (I-N) ,REAL*8 (A-H, O-Z) 

C 

COMMON /CONSTS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 pfd,aloge,alnio,coimin,coising 

C 

XI - PT/PO 
C 

IF (XI .LE. 1.58) GO TO 10 
IF (XI .LE. 3.16) GO TO 12 
IF (XI .LE. 10.) GO TO 14 
C 

DISC = -42.5 
GO TO 40 
C 

10 DISC = -12.0 * XI * XI 
GO TO 40 
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c 

12 DISC - -30. 

GO TO 40 
C 

14 DISC « -17 .5-25.*DLOG10(Xl) 

GO TO 40 
C 

40 RETURN 
END 
C 

SUBROUTINE FTNER1 (PR, P0 ,F, G, DISC) 

FSS EARTH RECEIVER PATTERN FROM CCIR REPORT 391-4 
ANTENNA DIAMETER 3 METERS, MAIN LOBE NOT GAUSSIAN 

IMPLICIT INTEGER*4(I-N) ,REAL*8(A-H,0-Z) 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 
1 PFD,ALOGE,ALN10,OOIMIN, GOISING 

DPR » PR ‘DEGREE 
DPO = P0 * DEGREE 
D - 3. 

WAVEL - 300. /F 
XI - D/WAVEL 

G1 ■ 2. + 15.* DLOGIO(XI) 

PM - 20 ./XI* DSQRT(G-Gl) 

PS - 15.85 / X1**0 .6 

IF (DPR .LE. PM) GO TO 50 
IF (DPR .LE. PS) GO TO 60 
IF (DPR .LE. 48.) GO TO 70 

DISC = -10. 

GO TO 80 

50 DISC « G- 2.5E-3 * X1*X1*DPR*DPR 
GO TO 80 

60 DISC = G1 
GO TO 80 

70 DISC - 3 2 .-25 . *DLOG10 (DPR) 

80 RETURN 
END 

SUBROUTINE PTNER2 (PR, P0 ,F, G, DISC) 

FSS EARTH REVEIVER PATTERN FROM CCIR REPORT 391-4 
MAIN LOBE GAUSSIAN, ANTENNA DIAMETER 3 METERS, 
MODIFIED FOR NON US COUNTRIES 

IMPLICIT INTEGER*4(I-N) ,REAL*8 (A-H, O-Z) 

C 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER,ERDB, EAP, 
1 PFD, ALOGE, ALN10 ,COIMIN, COISING 

C 

DPR - PR ‘DEGREE 
X - PR/P0 
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P2 - 10.“((32.+10.)/25.) 

C 

IF (X.LE.0.5) THEN 
DISC - G - 12 . *X*X 
GO TO 80 

ELSE IF (DPR .GE. P2) THEN 
DISC - -10. 

GO TO 80 
END IF 
C 

DISC - G - 12.*X*X 

D1 « 32. - 25 . ‘DLOGIO (DPR) 

IF (D1 .GE. DISC) DISC * D1 
C 

80 RETURN 
END 
C 

SUBROUTINE PTNER3 (PR, P0 ,F, G, DISC) 

FSS EARTH REVEIVER PATTERN FROM CQR REPORT 391-4 
MAIN LOBE GAUSSIAN, ANTENNA DIAMETER 4.5 METERS, 
MODIFIED FOR US ONLY 

IMPLICIT INTEGER*4 (I-N) ,REAL*8 (A-H, O-Z) 

COMMON / CON STS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 
1 PFD,AL0GE,ALN1 0 , CO I MIN, CO IS ING 

DPR = PR ‘DEGREE 
X * PR/PO 

P2 - 10.“((29.+10.)/25.) 

IF (X.LE.l.) THEN 
DISC = G - 12 . *X*X 
GO TO 80 

ELSE IF (DPR .GE. P2) THEN 
DISC ■ -10. 

GO TO 80 
END IF 

DISC - G - 12 . *X*X 
D1 « 29. - 25 . ‘DLOG1 0 (DPR) 

IF (D1 .GE. DISC) DISC » D1 

80 RETURN 
END 

SUBROUTINE PTNER4 (PR, P0 ,F, G, DISC) 

FSS EARTH REVEIVER PATTERN FROM CCIR REPORT 391-4 
MAIN LOBE GAUSSIAN, ANTENNA DIAMETER 4.5 METERS, 
MODIFIED FOR US CNLY 

IMPLICIT INTEGER*4 (I-N) ,REAL‘8(A-H,0-Z) 

C 

COMMON /CONSTS/ E, PI , RADIAN, DEGREE, GCR, ER, ERDB, EAP, 
1 PFD, ALOGE, ALN1 0 , COIMIN, OOIS ING 

C 

DPR - PR ‘DEGREE 
X « PR/PO 
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P2 - 10 .** ( (29.+10 .) /25.) 

C 

IF (X.LE.l.) THEN 
DISC - G - 12.*X*X 
GO TO 80 

ELSE IF (DPB . GE. P2) THEN 
DISC - -10. 

GO TO 80 
END IF 
C 

DISC - G - 12.*X*X 

D1 « 29. - 25 »*DLOG10 (DPR) 

IF (D1 .GE. DISC) DISC - D1 
C 

80 RETURN 
END 
C 

SUBROUTINE PTNER5(PR,P0,F,G,DISC) 

EARTH RECEIVER PATTERN FROM RARC-83 P.115 CURVE B 
BSS PATTERN 

IMPLICIT INTEGERM(I-N) ,REAL*8(A-H, O-Z) 

COMMON /CONSTS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 
1 PFD, ALOGE, ALN10 ,COIMIN, GOISING 

X « PR/PO 

IF (X .LE. 0.25) GOTO 10 
IF (X .LE. 0.94) GOTO 12 
IF (X .LE. 18.88) GOTO 14 

DISC - G -43.2 
RETURN 

10 DISC « G 
RETURN 

12 DISC « G -12. * X * X 
RETURN 

14 DISC - G -11.3 -25. * DLOGIO(X) 

RETURN 

END 

SUBROUTINE XPTNER1 (PR, P0 ,G, DISC) 

CROSS POLARIZATION RECEIVER PATTERN, BSS 

IMPLICIT REAL *8 (A-H, O-Z) 

C 

X - PR/PO 

IF (X. LE. 0 .25) GO TO 40 

IF (X. LE. 0 .44) GO TO 50 

IF (X.LE.l. 28) GO TO 60 

IF (X.LE.3 .22) GO TO 70 

IF (X.LE.5.60) GO TO 80 

IF (X. LE. 18 .88) GO TO 90 
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c 

DISC - G-43.2 
GO TO 100 
C 

40 DISC « G-25. 

GO TO 100 
C 

50 DISC « G-(30.+40.*DLOG10(1.0-X)) 
GO TO 100 
C 

60 DISC - G-20. 

GO TO 100 
C 

70 DISC « G-<17 .3+25. *DLOG10(X)) 

GO TO 100 
C 

80 DISC - G-30 . 

GO TO 100 
C 

90 DISC - G-(11.3+25.*DLOG10 (X) ) 

C 

100 RETURN 
END 
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APPENDIX B 


CONCAVE, QUASI -CONCAVE AND PSEUDO-CONCAVE FUNCTIONS 

The content in this appendix is from the book "Nonlinear 
Programming" written by 01 vi L. Mangasarian [1053. All the page numbers 
appearing in the content are referred to this book. Here only the 
concave functions are discussed; the same discussions apply to convex 
functions with obvious substitution. 

Definition of convex set (p. 39) 

A set IcR n is a convex set if, for x.yel, aeR, 0<a<l, one has 
(l-a)x+ayel , (B.l) 

here R n is the set of n-dimensional vector space, R is the set of real 
numbers. 

Definition of concave function (p. 56) 

A numerical function f defined on a set IcR n is said to be 
concave at xel if, for yel, 0<ta<l, (l-a)x+ayel, one has 

(l-a)f (x)+af (y) < f[(l-a)x+ay] ; (B.2) 

f is said to be concave on I if it is concave at each xel. 
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Definition of strictly concave function (p. 57) 

A numerical function f defined on a set IcR n is said to be strictly 
concave at xel if, for yel, xfy, 0<a<l, (l-a)x+ayel, one has 

(l-a)f (x)+af (y) < f[(l-a)x+ay] ; (8.3) 

f is said to be strictly concave on I if it is strictly concave at each 
xel. 

Definition of quasi -concave function (p. 132) 

A numerical function f defined on a set IcR n is said to be quasi -concave 
at xel if, for yel, f(x)<f(y), Oca<l, (l-a)x+ayel, one has 

f(x) < f[(l-a)x+ay] ; (B.4) 

f is said to be quasi -concave on I if it is quasi -concave at each xel. 

Definition of strictly quasi -concave function (p. 137) 

n 

A numerical function f defined on a set IcR is said to be strictly 
quasi -concave at xel if, for yel, f(x)<f(y), 0<a<l, (l-a)x+ayel, one has 

f(x) < f[(l-a)x+ay] ; (B.5) 

f is said to be strictly quasi -concave on I if it is strictly quasi - 


concave at each xel 



Theorem 1 (p. 139) 


Let f be a numerical function defined on the convex set I in R n , 
and let xel be a local maximum. If f is strictly quasi -concave at x, 
then f(x) is a global maximum of f on I. 

Definition of pseudo-concave function (p. 141) 

Let f be a numerical function defined on some open set in R n 
containing the set I. f is said to be pseudo-concave at xel if it is 
differentiable at x, and for yel, vf(x)(y-x)<0, one has 

f(y) < f(x) ; (B.6) 

f is said to be pseudo-concave on I if it is pseudo-concave at each xel. 

Theorem 2 (p. 143) 

Let I be a convex set in R n , and let f be a numerical function 
defined on some open set containing I. If f is pseudo-concave on I, 
then f is strictly quasi -concave on I and hence also quasi -concave on I. 
The converse is not true. 

Theorem 3 (p. 144) 

Let f be a numerical function defined on some open set I in R n , let 
xel, and let f be differentiable at x. If f is concave at x, then f 
is pseudo-concave at x, but not conversely. 
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Theorem 4 (p. 145) 


Let I be a convex set in R n , and let f be a numerical function 
defined on some open set containing I. If f is pseudo-concave on I, 
then each local maximum of f on I is also a global maximum of f on I. 

A final note: a function that is strictly concave is also concave, 

a function that is concave is also pseudo-concave, a function that is 
pseudo-concave is also strictly quasi -concave, a function that is 
strictly quasi -concave is also quasi -concave; the converse is not true. 
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PRU 


APPENDIX C 


CONTOUR PLOTS OF OBJECTIVE-FUNCTION SURFACE 

ARG=62 . BGL=68. 
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PRG 


ARG=62 . 


PRU = 68 . 



BOL 
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PRG = 62 . 


B0L=68 . 
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RRG=68. 





63. 
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APPENDIX D 


AS CALCULATION CODE 


C*** >> MAIN PROGRAM << 

THIS IS TO CALCULATE THE NECESSARY SATELLITE SPACING 
FOR WO SERVICE AREAS. 

INITIAL SATELLITE LOCATIONS SHOULD BE THE SAME. 

SUBROUTINES SAfffi AS MINI-SOUP PROGRAM, EXCEPT DIMENSION OF 
SERVICE AREA IS WO 

»*•***••»*******•••*•••**•*•***•**•*••**••***•********••******•***, 
IMPLICIT REAL *8 (A-H, 0- Z) 

CHARACTER *6 NAMES A 

COMMON /CON STS/ E , PI , RADI AN . DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN1 0 .COIMIN , OOIS ING 

OOMMON /PAR AMS/ NOMSAR . NAMESA(2 ,2) ,NTPSA<2) ,CPHIO ( 2) 

COMMON /VECTOR/ DSL0NC2) ,RSL0N(2) ,XD<2) ,*0(2) , 

1 ROAIKJ, 

2 X0AC(2) ,Y0AC<2) ,Z0AC(2> ,R0AC<2> 

COMMON /VARBLS/ FREQ < 2) ,IFOLAR < 2) ,GAINR< 2) ,GAINT<2) , 

1 EIRP<2) ,IPTNST<2) ,IPTNER<2) 

COMMON /MIN ELL/ BCLAT(2) ,BCL0N(2) ,DBCLAT(2) , DBCL0N(2) , 

1 REFLAT12) ,REFL0N<2) ,AXR(2) , 

I ORIENT (2) ,AXMAJ<2> 

COMMON /TPOINT/ RELON ( 2 .20) ,RELAT (2 .20) ,DELON<2 ,20) , 

1 DELAT (2.20) ,XE(2.20) ,YE<2,20) r ZE(2.20) 

COMMON /ANGLES/ YPHIT, YPHIR. PHITK, YPHIO 

COMMON /REAL/ PIKJ, PKKJ , WFQM, WDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

OOMMON /SR03/ DMGN( 2 >20) 

t*l**t*B*******************B*******«***t****t***i*****«*Bt*******«** 

OPEN (UNIT-20 .FILE* * rNPUTO.DAT* ,TYPE- * OLD* ) 

OPEN (UNIT-6 , FILE- * DSOUT. DAT' , TYPE-' NEW* ) 

CALL I CON ST 
CALL INPUTO 

DO 2 K - 1,NUMSAR 
CALL GAINER (K) 

EIRP(K) - PFD+10 . *DLOGl 0(4. *PI *ROAC (K) *ROAC (K) ) -ERDB 
2 CONTINUE 

CALL ZFUNCT 
CALL SEPAR 

STOP 

END 

SUBROUTINE ZFUNCT 

*** » THIS ROUTINE IS THE OVERALL CONTROL ROUTINE FOR << 
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t<**t*i**«t**t*******i*t**«*«o*****m*mo********ti*t***««»* t tt 
IMPLICIT REAL*8(A-H,0-Z) 

CHARACTER*6 NAMESA 

COMMON /CON STS/ E, PI, RADI AN. DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN1 0 , OOIMIN, GOISING 

OOMMON /PAR AMS/ NUMSAR, NAMESA(2. 2) ,NTPSA(2) , CPHI0(2) 

COMMON /VECTOR/ DSL0N(2) , RSLON(2) , X0(2) ,YO(2) , 

1 ROAIKJ, 

2 XOAC(2) , YOAC(2) , ZOAC(2) ,ROAC(2) 

COMMON /VARBLS/ FREQ (2) r IF0LAR(2) , GAINR(2) ,GAINT(2) , 

1 EIRP(2) ,IPTNST(2) ,IPTNER(2) 

COMMON /MIN ELL/ BCLAT(2) r BCLON(2) ,DBCLAT(2) , DBCLON(2) , 

1 REFLAT (2) ,REFLON(2) ,AXR(2) , 

1 ORIENT (2) ,AXMAJ(2) 

COMMON /TPOINT/ REL0N<2.20) ,RELAT(2.20) ,DELOtJ(2,20) , 

1 DEL AT (2,20) r XE(2,20) ,YE<2,20) ,ZE<2,20) 

COMMON /ANGLES/ YPHIT,YPHIR. PHITK, YPHIO 

COMMON /REAL/ PIKJ, PKKJ, FWFQM, FWDRCK, YPWDRC, YPWDRX, 

1 XO AKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 

OOMMON / SRCH/ DMGN(2.20) 

k ** >> INITIALIZE PARAMETER 

DO 10 1-1 .NUMSAR 
CALL REFCAL ( I ) 

10 CONTINUE 

k ** >> OUTER SUMMATION (OVER K) FOR ALL SERVICE AREAS << 

DO 1000 K « 1 , NUMSAR 
JNTP - NTPSA(K) 

WRITE(6,705) NAMESA(K,1) ,DSLON(K) 

*** >> MIDDLE SUMMATION (OVER J) FOR ALL TEST POINTS IN AREA K << 

DO 900 J - 1 , JNTP 
CALL KPHI (K,J) 

CALL XPHI0 (K, K, J) 

CALL KPWDRC (K,J) 

CALL XPWFQ (FREQ(K) , FREQ(K) ) 

CALL XPOWER (K,K, J) 

'** » CALCULATE INTERFERENCE POWER 

SUMP-0. 

DO 500 I - 1, NUMSAR 

IF (I.EQ.K) GO TO 500 
CALL ZPHI (I,K, J) 
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CALL XPHIO (I,K,J) 

IF (IPOLAR(I) .BQ.IFOLAR(K)) THEN 
CALL ZFWDRC (I,K,J) 

ELSE 

CALL ZIWDRX (I,K,J) 

END IF 

CALL XPWFQ (FREQ(K) , FREQ ( I ) ) 

CALL XPOWER (I,K,J) 

CO I - PKKJ-PI KJ 
DMGN(K,J) - OOI - COISING 
C 

WRITE (6 ,706) DELON (K,J) ,DELAT(K,J) ,NAMESA(I,1> , PIKJ , OOI, DMGNCK, J) 
500 CONTINUE 

C 

900 CONTINUE 
1000 CONTINUE 
C 

705 FORMAT (//,12X, 'TEST COUNTRY ,A6 ,5X, ' SATELLITE : » ,F8 .2 ,// ,15X, 

1 'TEST POINT' ,1 OX, ' INT. SAT. ' ,2X, ' INT. WR',3X, 'C/I (dB) ' , 

25X,' MARGIN') 

706 FORMAT (/ ,15X, F7.2 ,2X, F7.2 ,5X, A6 ,4X, F8 .2 ,4X, F6.2 ,6X, F6.2) 

RETURN 

END 


SUBROUTINE SEPAR 

THIS IS TO FIND THE MINIMAL REQUIRED SATELLITE SPACING FOR TWO 
SERVICE AREAS 

************************************** 


c 

c 


c 

c 


c 


c 


c 


c 

c 


c 


IMPLICIT REAL*8(A-H,0-Z) 
CHARACTER *6 NAffiSA 


COMMON /CON STS/ E, PI, RADIAN. DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN1 0 » 00 1 MIN . COISING 

COMMON /PAR AMS/ NUMSAR, NAMESA (2 .2) ,NTPSA<2) ,CPHI0 ( 2) 


COMMON /VECTOR/ DSLON ( 2) ,RSLON ( 2) ,XO (2) ,YO (2) , 

1 ROAIKJ, 

2 XOAC(2) ,YOAC(2) ,ZOAC(2) ,ROAC(2) 


COMMON /VARBLS/ FREQ ( 2) , IFOLAR ( 2) ,GAINR( 2) ,GAINT (2) , 
1 EIRP(2) ,IPTNST(2) ,IPTNER(2) 


COMMON /MINELL/ 

1 

1 


BCLAT(2) , BCLON ( 2) ,DBCLAT(2) ,I»CLON(2) , 
REFLAT (2) ,REFLON ( 2) ,AXR(2> , 

ORIENT (2) ,AXMAJ ( 2) 


COMMON /TPOINT/ RELON(2,20) ,RELAT(2,20) ,DELON<2,20) , 

1 DELAT(2,20) ,XE(2,20) ,YE(2.20) ,ZE(2,20) 


COMMON /ANGLES/ YPHIT, YPHIR, PHITK, YPHIO 

COMMON /REAL/ PIKJ , PKKJ , IWFQM, WDRCK, YPWDRC, YPWDRX, 

1 XOAKKJ , YOAKKJ , ZOAKKJ , ROAKKJ 
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COMMON / SRCH/ DMGN(2,20) 

******************************** 

NTP1 - NTPSA(l) 

NTP2 - NTPSA(2) 

SELECT WORST MARGIN AMONG ALL C/I VALUES AT ALL TEST POINTS 

RMIN1 - EMGN(l.l) 

DO J-2,NTP1 

RMIN1 - DMINKRMINl ,DMGN(1,J) ) 

END DO 

RMIN2 - DMGN( 2,1) 

DO J=2,NTP2 

RMIN2 = DM INI ( RMIN2 ,DMGN( 2, J) ) 

END DO 

NO SATELLITE SPACING REQQUIRED 

IF ( RMIN1 . GE. 0 . . AND. RMIN2 . GE. 0 . ) THEN 
WRITE (6,1) 

1 FORMAT (// ,11X, ' **** SATELLITE SEPARATION NOT NEEDED ****') 
RETURN 

END IF 

CALCULATE TOPOCENTRIC AND GEOCENTRIC ANGLE (SATELLITE SPACING) 
FROM WORST TEST POINT 

PRECAUSION MADE IF TWO SYSTEMS USE DIFFERQIT RECEIVING PATTERNS 

IF(IPTNERd) .EQ. IFTNER(2) )THEN 
IF (RMIN1.LE. RMIN2) THEN 
CALL MINMGN( 1 . J1 , RMIN1) 

WRITE (6,2) RMIN1 ,NAMESA(1,2) ,DELON(l,Jl) ,DELAT(1,J1) 

CALL FTNANGd,RMINl ,TOPOANG) 

CALL GEOANGt 2, 1 , J1 ,TOPOANG, DELTAS) 

ELSE 

CALL MINMGN( 2 , J2 ,RMIN2) 

WRITE (6,2) RMIN2 ,NAMESA(2,2) ,DELON(2,J2) ,DELAT(2,J2) 

CALL PTNANG(2,RHIN2,TOPOANG) 

CALL GEOANGd ,2 , J2 ,TOPOANG, DELTAS) 

END IF 
ELSE 

CALL MINMGNd , J1 ,RMIN1) 

WRITE (6,2) RMIN1 ,NAMESA(1,2) , DELON (1,J1) , DELAT(1,J1) 

CALL PTNANGd , RMIN1 , TOPOANG1) 

CALL GEOANG ( 2 ,1 , J1 ,TOPOANGl , DELTAS 1) 

CALL MINMGN(2,J2,RMIN2) 

WRITE (6,2) RMIN2,NAMESA(2,2) , DELON(2,J2) ,DELAT(2,J2) 

CALL PTNANG ( 2 , RMIN2 , TOPOANG2 ) 

CALL GEOANGd, 2 , J2 ,TOPOANG2 , DELTAS2) 

DELTAS « DMAX1 (DELTAS 1 , DELTAS2) 

END IF 
C 

WRITE(6, 100) DELTAS. DSLONd) ,00ISING 

2 FORMAT (//,1 IX, 'WORST MARGIN IS ',F6.2,' KP',A6,'(', 

1 F7 .2 , ' , ' ,F7 .2 , ' ) 1 ) 

100 FORMAT (/,11X, 'SATELLITE SEPARATION s ',F5.2,' AT ', 

1 F7.2 , ' FOR C/I ' ,F4.1 , ' dB') 

C 
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RETURN 

END 

C 

SUBROUTINE MINMGNUCl. J1 ,RMIN) 

THIS IS TO FIND THE MINIMUM MARGIN AMONG THE TEST POINTS 
IMPLICIT REAL*8(A-H,0-Z) 

CHARACTER *6 NAMES A 

COMMON /CON STS/ E, PI, RADIAN. DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10, COIMIN, OOISING 

COMMON /PAR AMS/ NUMSAR. NAMESA(2 .2) ,NTPSA (2) ,CPHI0 ( 2) 

COMMON / SRCH/ DMGN(2,20) 

DO J=1 , NTPSA(Kl) 

IF(DMGN(K1, J) . LE.RMIN) THEN 
J1 = J 
RETURN 
END IF 
END DO 
RETURN 
END 

SUBROUTINE PTNANG (K1 ,RMIN . TOPOANG) 


THIS IS TO CALCULATE THE TO PO CENTRIC ANGLE (FROM GROUND 
RECEIVING REFERENCE PATTERN) NECESSARY TO PROVIDE 
THE DISCRIMINATION 


C 

C 


C 

C 


C 

C 


C 


IMPLICIT REAL*8(A-H,0-Z) 
CHARACTER *6 NAMES A 


OOMMON /CON STS/ 

1 

OOMMON /PARAMS/ 
COMMON /VARBLS/ 

1 

COMMON /SRCH/ 


E, PI , RADIAN . DEGREE, GCR, ER, ERDB, EAP, 
PFD, ALOGE, ALN1 0 , COIMIN, OOISING 

NUMSAR, NAMES A (2, 2) ,NTPSA(2) ,CPHI0(2) 

FREQ (2) ,IPOLAR(2) ,GAINR(2) ,GAINT(2) , 
EIRP (2) ,IPTNST(2) ,IPTNER(2) 

DMGN(2,20) 


IF (IFOLAR(l) ,NE. IPOLAR(2) ) THEN 
CALL XRPTNER1(PR,CPHI0(K1) ,GAINR(K1) ,RMIN) 
GO TO 60 
END IF 


GOTO (51 ,53 ,55,57,59) ,IFTNER(K1) 

51 CALL RPTNER1 (PR, CPHIO (Kl) , FREQ(Kl) ,GAINR(K1) ,RMIN) 
GO TO 60 

53 CALL RPTNER2 (PR. CPHIO (Kl) , FREQ (Kl) ,GAINR(K1) ,RMIN) 
GO TO 60 

55 CALL RPTNER3 (PR. CPHIO (Kl) ,FREQ(K1) ,GAINR(K1) ,RMIN) 
GO TO 60 
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57 CALL RPTNER4 (PR. CPHIO (Kl) , FREQ(Kl) .GAINR(Kl) ,RMIN) 

GO TO 60 

59 CALL RPTNER5 (PR. CPHIO (Kl) , FREQ(Kl) ,GAINR(K1) ,RMIN) 

C 

60 TOPOANG - PR * DEGREE 
RETURN 

END 

C 

SUBROUTINE GEOANG ( I, K,J, TOPOANG, DELTAS) 

THIS IS TO ITERATE TO CALCULATE THE NECESSARY GEOCENTRIC 
ANGLE FOR THE REQUIRED DISCRIMINATION 


C 

C 

C 

C 


C 


C 

C 

C 

C 


C 

20 


IMPLICIT REAL*8(A-H,0-Z) 


CHARACTER *6 NAMES A 


COMMON /CON STS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALNIO.OOIMIN. OOISING 

COMMON /PAR AMS/ NUMSAR,NAMESA(2,2) ,NTPSA(2) ,CPHI0 (2) 


COMMON /VECTOR/ 

1 

2 


DSLON ( 2) , RSLON ( 2) ,XO (2) ,YO(2) , 
ROAIKJ , 

XOAC (2) (YOAC (2) ,ZOAC(2) ,ROAC(2) 


COMMON /MINELL/ 

1 

1 


BCLAT(2) , BCLON ( 2) ,DBCLAT(2) ,DBCLON(2) , 
REFLAT (2) ,REFLON(2) ,AXR(2) , 

ORIENT (2) ,AXMAJ ( 2) 


COMMON /TPOINT/ RELON ( 2 .20) ,RELAT (2 ,20) ,DELON ( 2 ,20) , 

1 DELAT(2,20) ,XE(2,20) ,YE(2,20) ,ZE(2,20) 


COMMON /SRCH/ DMGN(2.20) 


DSLONG = DSLON (I) 


IF (DBCLON (K) . GE. DBCLON (I ) ) THEN 
DELTA * 0.01 
DSK - DSLONG+ TOPOANG/ 2. 

DSI ■ DSLONG-TOPOANG/2. 

ELSE 

DELTA - -0.01 

DSK »= DSLONG- TOPOANG/ 2 . 

DSI ■ DSL ONG+ TOPOANG/ 2. 

END IF 


DSK - DSK-DELTA 

DSI - DSI+ DELTA 

RSK « DSK “RADIAN 

RSI - DSI “RADIAN 

XOAK - XE(K,J) -GCR*DCOS(RSK) 

YOAK « YE (K, J) -GCR*DSIN (RSK) 

ZOAK - ZE(K,J) 

ROAK - DSQRT (XOAK*XOAK+YOAK*YOAK+ZOAK*ZOAK) 

XOAI - XE(K,J)-GCR*DCOS(RSI) 

YOAI - YE (K, J) -GCR*DSIN(RSI) 

ZOAI - ZE (K, J) 

ROAI - DSQRT (XOAI*XOAI+YOAI*YOAI+ZOAI*ZOAI) 

ARG - (XOAK*XOAI+YOAK*YOAI+ZOAK*ZOAI) / (ROAK*ROAI) 
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TOPO - DA COS (ARG) ‘DEGREE 
IF(TOPO. GT. TOPOANG) GOTO 20 
DELTAS - DABS ( DSK-DS I ) +0 . 0 2 
RETURN 
END 

SUBROUTINE RPTNER1 (PR, P0 ,F, G, DISC) 

FSS EARTH RECEIVER REFERENCE PATTERN FROM CQR REPORT 391-4 
ANTENNA DIAMETER 3 METERS, MAIN LOBE NOT GAUSSIAN 

IMPLICIT REAL *8 (A-H, O-Z) 

COMMON /CON STS/ E, PI, RADIAN, DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10,OOIMIN> OOISING 

D ■ 3. 

WAVED *= 300. /F 
XI - D/WAVEL 

G1 = 2. + 15.* DLOGIO(XI) 

IF (DISC. GE. ( Gl-G) ) GO TO 50 
IF (DISC. GE. (-G-10 .) ) GO TO 60 

PR - PI 
GO TO 80 

50 PR = RADIAN * DSQRT(-DISC*400./ (X1*X1) ) 

GO TO 80 

60 PR « RADIAN *10. ** (-(DISC+G-32 .) /25. ) 

80 RETURN 
END 

SUBROUTINE RFTNER2 (PR, P0 , F, G, DISC) 

FSS EARTH RECEIVER PATTERN FROM CQR REPORT 391-4 
MAIN LOBE GAUSSIAN, ANTENNA DIAMETER 3 METERS, 

MODIFIED FOR NON US COUNTRIES 

IMPLICIT REAL *8 (A-H, O-Z) 

C 

COMMON /CON STS/ E, PI, RADIAN. DEGREE, GCR, ER, ERDB, EAP, 

1 PFD, ALOGE, ALN10,COIMIN. OOISING 

C 

IF (DISC. GE. -3 . ) THEN 
PR - P0 * DSQRT( -DISC/1 2.) 

GO TO 80 

ELSE IF (DISC. LE. - ( GflO .) ) THEN 
PR - PI 
GO TO 80 
END IF 
C 

PR1 - P0 * DSQRT( -DISC/1 2.) 

PR2 - RADIAN * 10. ** (-( Gf DISC-32 .) /25. ) 

PR = DMAX1 ( PR1 , PR2 ) 

C 

80 RETURN 
END 
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SUBROUTINE RETNER3 (PR. PO ,F, G, DISC) 

PSS EARTH RECEIVER PATTERN FROM CCXR REPORT 391-4 
MAIN LOBE GAUSSIAN. ANTENNA DIAMETER 4.5 METERS, 
MODIFIED FOR US CNLY 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON /CON STS/ E, PI, RADI AN. DEGREE, GCR, ER, ERDB, EAP, 
1 PFD, ALOGE,ALN10,COIMIN, OOISING 

IF (DISC. GE. -12.) THEN 
PR * PO * DSQRT ( -DIS C/1 2 . ) 

GO TO 80 

ELSE IF (DISC. LE. -(GflO.) ) THEN 
PR « PI 
GO TO 80 
END IF 

PR1 ■ PO * DSQRT ( -DIS C/1 2.) 

PR2 - RADIAN * 10. ** (-(GfDISC-29.) /25.) 

PR - DMAX1 ( PR1 , PR2 ) 

80 RETURN 
END 

SUBROUTINE RETNER4 (PR. PO ,F, G, DISC) 

FSS EARTH RECEIVER PATTERN FROM COR REPORT 391-4 
MAIN LOBE GAUSSIAN. ANTENNA DIAMETER 4.5 METERS. 
MODIFIED FOR US CNLY 

IMPLICIT REAL*8(A-H» O-Z) 

wnflUw /U/BBIO/ El, r± , kauiaH . LiKfeKSG, ViUC. £K, £KUU, bAP, 

1 PFD, ALOGE,ALN10.OOIMIN. OOISING 

IF (DISC. GE. -12.) THEN 
PR “ PO * DSQRT (-DISC/ 12.) 

GO TO 80 

ELSE IF (DISC. LE.-(Gfl0.) ) THEN 
PR - PI 
GO TO 80 
END IF 

PR1 - P0 * DSQRT ( -DIS C/1 2.) 

PR2 ■= RADIAN * 10. ** (-(GfDISC-29.) /25.) 

PR - DMAX1 ( PR1 , PR2 ) 

80 RETURN 
END 

SUBROUTINE RFTNER5 (PR, P0 ,F, G, DISC) 

BSS PATTERN, RARC-83 . CPM, P.115 CURVE-B, 

ANTENNA DIAMETER 1 METER 

IMPLICIT REAL*8(A-H. O-Z) 
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COMMON /CON STS/ E , PI , RADIAN . DEGREE, GCR, ER, ERDB, EAP 
1 PFD, ALOGE, ALN1 0 » OOI MIN# OOISING 

IF (DISC. GE. -10. 6032) THEN 
PR - PO * DSQRT( -DISC/1 2.) 

GO TO 80 

ELSE IF (DISC. LE. -43. 2) THEN 
PR - PI 
GO TO 80 
END IF 

PR - PO * 10. ** (-(DISC+11 .3) /25. ) 

30 RETURN 
END 

SUBROUTINE XRPTNER1 (PR, PO ,G, DISC) 

CROSS POLARIZATION RECEIVER PATTERN 
IMPLICIT REAL*8(A-H#0-Z) 

X » PR/PO 

IF (DISC. GE. -20.) GO TO 60 
IF (DISC. GE. -30.) GO TO 70 
IF (DISC. GE. -43. 2) GO TO 80 

PR - PI 
GO TO 100 

60 PR - 0. 

GO TO 100 

70 PR - 10.**(-(17.3+DISC)/25.) * PO 
GO TO 100 

BO PR - 10 . ** (-(11.3+DISC) /25. ) * PO 


100 RETURN 



APPENDIX E 


FORMULATIONS OF MIXED-INTEGER AND LINEAR PROGRAMS 
A. ALGORITHMS 

A mixed integer program (MIP) can be solved by a branch-and-bound 
algorithm [88], In such an algorithm, a linear objective function is to 
be optimized; the constraints of the problem are expressed as linear 
equalities or inequalities. Some, but not necessarily all, of the 
decision variables are integers. The set of feasible solutions that 
satisfy the constraints constitutes the feasible region. For our 
purpose, we can assume that the integer variables are bounded. Because 
the integer variables have a finite number of feasible values, the 
number of feasible solutions is finite; therefore an enumerative 
approach can be used to theoretically test all the feasible solutions in 
order to find the globally optimal solution(s). 

To perform the enumeration by the branch-and-bound algorithm, the 
set of feasible solutions can be structured as a tree, and every branch 
represents one possible value of a particular integer variable. For the 
branch-and-bound concept, "branching" means testing a path that leads to 
a subset of feasible solutions, "bounding" means calculating the upper 
and lower bounds of the objective function value associated with the 
tested path. In the enumerative process, the upper and lower bounds of 
the objective function value are updated whenever more favorable values 
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are found in the hranch-and-bound process. In the branch-and-bound 
process a path is tested against the updated upper and lower bounds: it 

finds a more favorable upper or lower bound of the objective function 
value and updates it, or it terminates testing that path and all its 
associated feasible solutions once it determines that this path can not 
yield a more favorable upper or lower bound of the objective function 
value. The process terminates when the updated upper and lower bounds 
are equal, or when it determines that an optimal solution does not 
exist. All the feasible solutions will have been considered implicitly, 
but hopefully very few will have been examined explicitly. A globally 
optimal solution is guaranteed, if one exists, by this process. 

The linear program (LP) is commonly solved by the simplex method 
[81]. In this program a linear objective function is to be optimized; 
the constraints of the problem are expressed as linear equalities or 
inequalities. All the decision variables are continuous variables, and 
they must have nonnegative values. The set of feasible solutions 
constitutes a feasible region that is a convex set; its boundaries are 
the hyperplanes representing the linear constraints and nonnegativity 
restrictions. Because the objective function and the constraints are 
all linear equations, a locally optimal solution is always at a vertex 
which is the intersection of the bounding hyperplanes. The simplex 
method examines a sequence of locally optimal feasible solutions of the 
linear program. Each solution examined shares at least one boundary 
with the previous one, and has an objective function value no less 
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favorable than that of the previous solution. The process terminates 
when it is determined that no improved solution can be found. 

The simplex method can be modified to handle nonlinear 
complementarity constraints through the use of restricted basis entry: 
when a complementary variable enters the basis at a non-zero value, its 
complement is forced to be non-basic and cannot enter the basis (except 
when the variable leaves). The effect of this modification is that the 
continuous feasible region is divided into many distinct subregions. 
Although the simplex method can still work, it only guarantees a local, 
but not necessarily a global, optimum, provided a feasible solution is 
found. 

As the number of satellites increases, the computational time 
needed to find a solution typically increases exponentially for the MIP 
technique, but only polynomially for the restricted basis LP (RBLP) 
technique. Therefore, the MIP technique may take a prohibitively lono 
time to solve a large problem; and the RBLP technique becomes an 
acceptable alternative, even though sub-optimal, rather than global ly- 
optimal, solutions may be found. 

B. PARAMETERS AND VARIABLES 

The coordinate system here is in the reverse longitude direction of 
the common global system, hence the longitude values increase as one 
moves in the westerly direction. Or simply speaking, the new coordinate 
system uses the magnitude of west longitude; in this way, there are only 
non-negative variables in both of the new formulations. 
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The parameters in the formulations are: 

W j( E j) : westernmost (easternmost) feasible location for satellite 
j in degrees west, 

d. : preferred location for satellite j relative to E., 

3 J 

AS ij : rec l u ^ re< l satellite separation between satellites i and j, 

1, if satellite i is west of j, 

x ii = { 

J 0, otherwise, 

m : number of satellites. 

E = min {E }, 
j 3 

W = max {W.}. 
j J 

The decision variables are: 

Xj : relative location of satellite j with respect to E^. , 
p i j ( n i j) : degrees west(east) of satellite i that satellite j is 
located, 

x j + ( x j~) : degrees west(east) of its preferred location that satellite 
j is located, 

Y : length of the occupied orbital arc. 
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Note that the nonnegative variable alone can not represent the 

separation between satellites i and j, and another nonnegative variable, 

n.j, is needed when i is east of j. Therefore when satellite i is west 

of j, the value of p^. is positive, and n^ should be set equal to zero; 

when satellite i is east of j, n. . is positive and p. . should be zero. 

' J 1 J 

Also note from the definition of x., the coordinate system for every 

J 

satellite is re-originated at the eastern-most bound, Ej. 


C. FORMULATION I 


If the optimal criterion is to minimize the total amount of orbital 
arc occupied by the satellites to be synthesized, the MIP formulation 
is 


Minimize f = Y 
subject to 


x. 


. - x. - p. . + n. . = 0, 

1 1 r ll 1 1 * 


P i j + n i j > AS 1j’ 

x. < W. - E., 

J J J 

P ij + n ij - Y < °* 
x. > 0, 

P ij’ n ij > °* 

Y > 0, 

(E - W)- Xlj + Pij < 0, 

(W - E)^. + n 1d < (W - E), 

j e {0,1}, 


for i < j 
for i < j 
for all j 
for i < j 
for all j 
for i < j 

for i < j 
for i < j 
for i < j 


(E.l) 


(E. 2) 

(E.3) 

(E.4) 

(E.5) 

(E.6) 

(E.7) 

(E.8) 

(E.9) 

(E. 10) 

(E. 11) 
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(E. 12) 


Equations (E.9) and (E.10) guarantee that 
< 0, or n^. < 0. 

Together with Equation (E.7), one has 

p.j = 0, or n.j = 0. (E. 13) 

If the optimal solution can be found, the optimal values of the 
variables x.'s specify the optimal orbital locations for the satellites. 
Otherwise, the code will tell the user that a feasible solution does not 
exi st. 

It is possible that the objective is only to have a feasible 
solution that satisfies the C/I protection requirements, as described in 
Chapter IV. Then one could reformulate the objective function of 
Equation (E.l) and ignore the variable Y, Equations (E.5) and (E.8). 

The calculation process will either stop when it finds a feasible 
solution, or determines that none exists. 

Also, the same formulation (Equations (E.l) through (E. 11) ) is 
still applicable if the system requirement is to have the maximum C/I 
results for all the service areas, and there is no interest in the 
conservation of the orbit resource. One could progressively adjust the 
C/I protection requirement and repeat the AS and MIP calculations until 
the resulting scenario uses up the whole feasible arc. This result is 
the scenario that offers the maximum C/I. 
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D. FORMULATION II 


For the same objective as stated above, the LP formulation with the 
nonlinear side constraints is 


Minimize f = Y 




(E. 14) 

subject to 





x i - x j - p u + "u - °- 

for 

i < 

j 

(E. 15) 

p fj + "ij * “ir 

for 

i < 

j 

(E.16) 

x j<V E r 

for 

all 

j 

(E. 17) 

p.j - Y < 0, 

for 

i < 

j 

(E. 18) 

x j » °* 

for 

all 

j 

(E. 19) 

p ij* n 1J > °- 

for 

1 < 

j 

(E. 20) 

V > 0, 




(E.21) 

■ 0. 

for 

i < 

j 

(E.22) 


The solution to this problem is an orbital assignment. 


In this formulation. Equation (E.22) is not a linear equation, 
thus the simplex method needs to be modified through the use of 
restricted basis entry: the variable p^ can not be a basic variable if 
n^. is a basic variable, and vice versa. Although the simplex method 
can still work, it only guarantees a locally optimal solution, provided 
a feasible solution is found. 
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E. FORMULATION III 


It is possible that every administration has a preferred satellite 
location and hopes the actual assigned location will be near the 
preferred location. Solutions become less and less attractive as the 
actual location becomes further and further removed from the preferred 
location. A suitable objective is to minimize the sum of the absolute 
deviations of the satellites' actual locations from their preferred 
locations, i.e., to minimize the total deviation. The MIP formulation 
is 


Minimize f = Z (x; + + x-j") 
i J J 



(E .23) 

subject to 




x j “ x j + + X J" = d j’ 

for 

all j 

(E.24) 

x i ' x j “ Pij + n ij = °» 

for 

i < j 

(E.25) 

p ij + n ij > AS ij» 

for 

i < j 

(E .26) 

x j < W j ~ E j * 

for 

all j 

(E.27) 

x j> x j + > x j" > °» 

for 

all j 

(E .28) 

Pij, nij > 0, 

for 

i < j 

(E .29) 

(E - W) • x 1 j + Pij < 0, 

for 

1 < J 

(E.30) 

(W - E) • Xi j + nij < (W - E), 

for 

i < j 

(E.31) 

x i j e {0,1), 

for 

i < j 

(E.32) 

The optimal values of the variables 

Xj prescribe optimal 

locations for 

the satellites which minimize the total 

deviation of the 

prescribed 


locations from the preferred locations. 
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(E.12) 


Equations (E.9) and (E.10) guarantee that 
p ij < °- or n ij 4 °- 


Together with Equation (E.7), one has 

p. . =0, or n. . =0. 
ij ij 


(E .13) 


If the optimal solution can be found, the optimal values of the 

variables x. specify the optimal orbital locations for the satellites. 

J 

Otherwise, the code will tell the user that a feasible solution does not 
exist. 

It is possible that the objective is only to have a feasible 
solution that satisfies the C/I protection requirements, as described in 
Chapter IV. Then one could reformulate the objective function of 
Equation (E.l) and ignore the variable Y, Equations (E.5) and (E.8). 

The calculation process will either stop when it finds a feasible 
solution, or determines that none exists. 

Also, the same formulation (Equations (E.l) through (E . 11 ) ) is 
still applicable if the system requirement is to have the maximum C/I 
results for all the service areas, and there is no interest in the 
conservation of the orbit resource. One could progressively adjust the 
C/I protection requirement and repeat the AS and MIP calculations until 
the resulting scenario uses up the whole feasible arc. This result is 
the scenario that offers the maximum C/I. 
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Table E.l 


Satellite preferred locations of six 
administrations 


country 

ARG 

BOL 

CHL 

PRG 

PRU 

URG 

case 1 

95 

95 

95 

95 

95 

95 

case 2 

no 

110 

110 

no 

no 

no 

case 3 

87.5 

92.5 

97.5 

87.5 

102.5 

82.5 


Table E.2 

AS parameters of six administrations 


AS 

ARG 

BOL 

CHL 

PRG 

PRU 

URG 

ARG 

* 

4.17 

4.19 

4.32 

1.41 

4.14 

BOL 


* 

4.57 

4.04 

4.26 

0.94 

CHL 



* 

2.00 

3.04 

1.59 

PRG 




* 

1.10 

2.46 

PRU 





* 

0.37 

URG 






* 
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The AS matrix for these six countries is in Table E.2. The 
resulting assignments, together with the total occupied arc Y, the total 
deviation Z of the prescribed locations from the preferred locations, 
and the C.P.U. times in seconds are listed in Table E.3. 

It is clear that the solutions obtained by solving the MIP 
formulation are better than those found by solving the LP formulation. 
The MIP solutions can never be worse than the corresponding LP 
solutions, but more computer time is required to solve the MIP 
formulation. 
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Table E.3 



Mixed-integer and 

1 linear 

program 

results 



case 

1 

case 2 

case 

3 

satellite 

MIP 

LP 

MIP 

LP 

MIP 

LP 

ARG 

88.68 

105.74 

101.35 

110.00 

88.76 

101.26 

BOL 

99.57 

101.57 

97.18 

104.33 

92.93 

92.50 

CHL 

95.00 

97.00 

105.54 

99.76 

97.50 

97.07 

PRG 

93.00 

95.00 

107.54 

97.76 

84.44 

87.50 

PRU 

91.06 

93.06 

109.63 

108.59 

102.50 

102.67 

URG 

96.59 

92.54 

110.00 

105.86 

81.98 

82.50 

deviation 

18.42 

23.71 

28.76 

33.69 

5.27 

14.36 

arc 

10.89 

13.20 

12.82 

12.24 

20.52 

20.17 

cpu(sec)* 

25.23 

1.31 

13.39 

1.30 

2.86 

1.25 


* IBM-3081 computer 
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APPENDIX F 


C/I CALCULATION OF NIXED INTEGER PROGRAM RESULT 


COUNTRY 

SATELLITE (LON. ) 

FREQUENCY 

ARG 

-88 . 68 

4000 . 00 

EOL 

-99.57 

4000 . 00 

CHL 

-95.00 

4000.00 

r rg 

-93.00 

4000.00 

f RU 

-91 .06 

4000.00 

URG 

-96.59 

4000.00 


TEST COUNTRY : ARG SATELLITE : -3S.68 


TEST POINT 


LON. 

LAT. 

-6 3. 

20 

-21 . 

80 

-ez . 

20 

-21. 

80 

-65. 

20 

-21 . 

80 

-6 5. 

20 

-21 , 

80 

-66. 

20 

-21 . 

80 

-65. 

20 

-21 . 

80 

-63 . 

80 

-22. 

.00 

-63. 

80 

-22. 

,00 

-62. 

30 

-22. 

.00 

- c 2 . 

80 

-22. 

,00 

-6 3. 

B0 

-22. 

,00 

-6 2 . 

80 

-22. 

.00 

-53. 

80 

-27 . 

, 20 

-53 . 

80 

-27 . 

. 20 

-53 . 

80 

-27. 

.20 

-53 . 

80 

-27. 

,20 

-53. 

80 

-27. 

.20 

-53. 

80 

-27. 

.20 

-56, 

70 

-36, 

.90 

-55. 

70 

-36, 

.90 

-56 . 

70 

-36, 

.90 

-55. 

70 

-36. 

.90 

-55. 

70 

-36. 

.90 

-55. 

70 

-36, 

.90 

-63. 

80 

-54. 

,70 

-S3. 

80 

-54, 

.70 

-62. 

80 

-54, 

,70 

-63. 

80 

-54. 

.70 

-63. 

80 

-54. 

.70 

-63. 

80 

-54. 

.70 

-68. 

30 

-54. 

.80 

-63. 

30 

-54. 

.80 

-63. 

30 

-54. 

.80 

-63 . 

30 

-54. 

.80 

-63. 

30 

-54, 

.80 

-63 . 

30 

-54. 

.80 

-73. 

20 

-50. 

,90 

-73, 

20 

-50, 

,90 

-73. 

20 

-50. 

.90 

-73, 

20 

-50 

.90 

-73. 

20 

-50. 

.90 


NT. SAT. 

C/I <dl 

BOL099 

40.83 

CHL095 

37.25 

PRG093 

39.39 

PRU091 

41.47 

URG096 

68.02 

TOTAL 

33.39 

BOL099 

41.37 

CHL095 

51.35 

PRG093 

32.03 

PRU091 

53.86 

URG096 

67.24 

TOTAL 

31.48 

BOL 099 

53.05 

CM! 03 5 

53 \ A 

PRG093 

31 .53 

PRU091 

52.74 

URG096 

49.37 

TOTAL 

31.39 

BOL099 

69.73 

CHL095 

55.57 

PRG093 

59.61 

PRU091 

53.45 

URG096 

41 . 77 

TOTAL 

41.25 

BOL099 

72.73 

CHL095 

42.01 

PRG093 

67.29 

PRU091 

52.85 

URG096 

74.65 

TOTAL 

41.65 

BOL099 

73.38 

CHL095 

39.26 

PRG093 

63.21 

PRU091 

53.16 

URG096 

75.74 

TOTAL 

39.08 

BOL099 

73.25 

CHL095 

35.44 

PRG093 

68.41 

PRU091 

52.91 

URG096 

75.82 


) MARGIN (dB) 

10.83 

7.25 

9.39 
11.47 
38*02 

8.39 

11.37 
21.35 

2.03 
23 .86 

37.24 
6.48 

23.05 
33 # i a 
1 ^53 
22.74 

19.37 

6.39 

39.73 
25.57 
29.61 
23.45 
11.77 

16.25 

42.73 
12.01 
37.29 
22.85 

44.65 

16.65 

43.38 

9.26 
38.21 
23.16 

45.74 
14.08 

43.25 
5.44 

38.41 

22.91 

45.82 


237 



original page is 
OF POOR QUALITY 


-73.20 

-50.90 

TOTAL 

35.36 

10.36 

-71.40 

-39 .00 

BOL039 

70.45 

40.45 

: .40 

-39.00 

CI0.095 

35.-15 

5.45 

-71.40 

-39.00 

PR3093 

6 5 . 6 

35.76 

-71.40 

-39.00 

PR 3091 

5 1.18 

24.13 

-71.40 

-39.00 

URG036 

71 . ?6 

41 .36 

- 71 .40 

-39.00 

TOTAL 

35 . 18 

10.38 

r . 50 

-31 .40 

BOL039 

70.35 

40.35 

-7./. 50 

-31 . 40 

CHL095 

35 . 10 

5.10 

-70.50 

-31.40 

PR 0093 

61.53 

31.63 

- 7. / . 50 

-31.40 

PRU091 

5-< .08 

24 ,:78 

-70.50 

-31.40 

URG096 

63.29 

38 . 2 9 

-70.50 

-31.40 

TOTAL 

35.03 

10.03 

-63.60 

-24.80 

BOL099 

49.00 

19.00 

-63.60 

-24.80 

CHL095 

35.23 

5.23 

-63.60 

-24.80 

PR '309 3 

59.37 

29.37 

-63.60 

-24.80 

PRU091 

53.21 

23.21 

-63.60 

-24.80 

URG096 

67.38 

37.38 

-63.60 

-24.80 

TOTAL 

34.96 

9.96 


TEST COUNTRY 

s BOL 

SATELLITE t 

-99.57 


TEST POINT 

I NT. SAT. 

C/I <dB> 

MARGIN 

-65.00 

-12.20 

ARG088 

70.32 

40.32 

-65.00 

-12.20 

CHL095 

61 .01 

31 .01 

-65.00 

-12.20 

P n 0093 

64.35 

34.85 

-65.00 

-12.20 

PRU091 

45.32 

15.32 

-65.00 

-12.20 

URG096 

66.08 

36.08 

-65.00 

-12.20 

TOTAL 

45.11 

20. 11 

-65.50 

-9.80 

ARG088 

69.00 

39.00 

-65.50 

-9.80 

CHL095 

59.68 

29.68 

-63.50 

-9.80 

PRG093 

63.52 

33.52 

-65.50 

-9.80 

PRU091 

42.75 

12.75 

-63.50 

-9.80 

URG096 

66.37 

36.37 

-63.50 

-9.80 

TOTAL 

42.60 

17.60 

-69.00 

-11.20 

ARG088 

69.81 

39.81 

-69.00 

-11.20 

CHL095 

60.49 

30.49 

-69.00 

-11.20 

PRG093 

64.34 

34.34 

-69.00 

-11.20 

PRU091 

38.97 

8.97 

-69.00 

-11.20 

URG096 

66.93 

36.93 

-69.00 

-11.20 

TOTAL 

38.92 

13.92 

-60.00 

-16.10 

ARG088 

47.93 

17.93 

-60.00 

-16.10 

CHL095 

60.57 

30.57 

-60.00 

-16.10 

PRG093 

43.69 

13.69 

-60.00 

-16.10 

PRU091 

67.54 

37.54 

-50.00 

-16.10 

URG096 

62.26 

32.26 

-60.00 

-16.10 

TOTAL 

42.18 

17.18 

-57.50 

-18.00 

ARG088 

41.71 

11.71 

-57.50 

-18.00 

CHL095 

59.57 

29.57 

-57.50 

-18.00 

PRG093 

38.56 

8.56 

-57.50 

-18.00 

PRU091 

66.53 

36.53 

-57.50 

-18.00 

URG096 

59.66 

29.66 

-57.50 

-18.00 

TOTAL 

36.80 

11.80 
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9«ejNAL PAGE IS 

w P00 '* quality 


-6 7 

♦ 50 

-22.70 

-6 '* 

.50 

-22.70 

-6 ' 

. 50 

-22.70 

-6 7 

.50 

-22.70 

-c 

. 50 

-22.7 J 

-e 

.50 

-22.70 


TEST COUNTRY : CHL 

TEST 

POINT 

-6 > . 50 

-17.50 

-62.50 

-17.50 

-63.50 

-17.50 

-63.50 

-17.50 

-6^.50 

-17.50 

-63.50 

-17.50 

-67. 10 

-23.00 

-67.10 

-23.00 

-67.10 

-23.00 

-67. 10 

-23 .00 

-67. 10 

-23.00 

-67.10 

-23.00 

-70.00 

-34.20 

-70.00 

-34.20 

-70.00 

-34.20 

-70.00 

-34.20 

-70.00 

-34.20 

-70.00 

-34.20 

-71.70 

-43.20 

-71.70 

-43.20 

-71.70 

-43.20 

-71 .70 

-43.20 

-71.70 

-43 . 20 

-71.70 

-43.20 

-72.80 

-51 .30 

-72.80 

-51 .30 

-72.80 

-51 .30 

-72.80 

-51.30 

-72.80 

-51 .30 

-72.80 

-51.30 

-75.70 

-46.80 

-75.70 

-46.80 

-75.70 

-46.80 

-75.70 

-46.80 

-75.70 

-46.80 

-75.70 

-46.80 

-74.00 

-28.90 

-74.00 

-28.90 

-74.00 

-28.90 

-74.00 

-28.90 

-74.00 

-28.90 

-74.00 

-28.90 

-70.40 

-18.30 

-70.40 

-18.30 


AP.G0S8 
CHL095 
PP.G09 3 
PRL091 
URG096 
TOTAL 


SATELLITE j 

INT. SAT. 
ARG083 
BOL099 
PR 009 3 
PRU09 1 
URC096 
TOTAL 

ARG088 

BOL099 

PRG093 

PRU091 

URG096 

TOTAL 

ARG088 

BOL099 

PRG093 

PRU091 

URG096 

TOTAL 

ARG088 
BOL099 
PRG093 
PRU09 1 
URG096 
TOTAL 

ARG088 

BOL099 

PRG093 

PRU091 

URG096 

TOTAL 


41.44 

11.44 

31.31 

1.91 

49 . 

:5 

19.35 

54.71 

24.71 

5 7 . JjJ 

29 .00 

3 1 . 

3 5 

6 . JS 


-35.00 
C/I (dB) 

MARGIN 

62.92 

32 . 92 

30. 38 

0 . 38 

50.44 

. 20.44 

30.53 

0.33 

55.17 

25.17 

27.52 

2.52 

35.33 

5.33 

32.18 

2. 18 

35.49 

5 .49 

48.86 

18.86 

49.75 

19.75 

29.19 

4. 19 

35.39 

5.39 



61.16 

31 ! 16 

54.61 

24.61 

59 . 79 

29 . 79 

51.30 

21.30 

35.20 

10. 20 

35.24 

5.24 

62 . 26 

32 . 26 

58.53 

28.53 

59.09 

29.09 

56.00 

26.00 

35. 15 

10. 15 

34.91 

4.91 

63.54 

33.54 

59.73 

29.73 

58.11 

28.11 

58.10 

28. 10 

34.85 

9.85 


ARG088 

BOL099 

PRG093 

PRU091 

URG096 

TOTAL 


36.11 

63.76 

60.28 

58.44 

58.50 

36.04 


6.11 

33.76 

30.28 

28.44 

28.50 

11.04 


ARG088 41 

BOL099 62 

PRG093 55 

PRU091 60 

URG096 55 

TOTAL 41 


43 

11.43 

06 

32.06 

97 

25.97 

69 

30.69 

35 

25.35 

03 

16.03 


ARG088 

BOL099 


63.43 
31 .30 


33.43 

1.80 
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-70.40 

-18.30 

PRG093 

50.94 

20.34 

-70.40 

-18.30 

PRU091 

31 .07 

1 .07 

-70.40 

-18.30 

URC096 

55.51 

25.51 

-70.40 

-18.30 

TOTAL 

28.37 

3.37 


COUNTRY 

: PRG 

SATELLITE : 

-93. .90 


TEST POINT 

INT. SAT. 

C/I 'dB> 

MARGIN 

-57.60 

-25.30 

ARG088 

32.10 

2. 10 

-57.60 

-25.30 

BOL099 

40.4 1 

10.41 

-57.60 

-25.30 

CHL095 

52.58 

22.58 

-57.60 

-25.30 

PRU091 

52.47 

22.47 

-57.60 

-25.30 

URG096 

46.66 

16.66 

-57.60 

-25.30 

TOTAL 

31 .31 

6.31 

-58.60 

-27.30 

ARG088 

30.07 

0.07 

-53.60 

-27.30 

BOL099 

43.33 

13.33 

-53.60 

-27.30 

CHL095 

50.95 

20.95 

-53.60 

-27.30 

PRU091 

50.34 

20.84 

-53.60 

-27.30 

URG096 

35.20 

5 . 20 

-53.60 

-27.30 

TOTAL 

28.70 

3.70 

-56.20 

-27.20 

ARG088 

31.10 

1 . 10 

-56.20 

-27.20 

BOL099 

45.13 

15. 13 

-55.20 

-27.20 

CHL095 

51 .31 

21.31 

-55.20 

-27.20 

PRU091 

51.20 

21.20 

-55.20 

-27.20 

URG096 

37.50 

7 . 50 

-55.20 

-27.20 

TOTAL 

30.00 

5.00 

-54.70 

-25.50 

ARG088 

31.36 

1 .86 

-54.70 

-25.50 

BOL099 

43.23 

13.23 

-54.70 

-25.50 

CHL095 

51.46 

21 .46 

-54.70 

-25.50 

PRU091 

51.34 

21.34 

-54.70 

-25.50 

URG096 

49.09 

19.09 

-54.70 

-25.50 

TOTAL 

31.39 

6.39 

-54.20 

-24. 10 

ARG088 

31.72 

1.72 

-54.20 

-24.10 

BOL099 

41.60 

11.60 

-54.20 

-24.10 

CHL095 

51.03 

21 .03 

-54.20 

-24.10 

PRU091 

50.91 

20.91 

-54.20 

-24.10 

URG096 

57.11 

27.11 

-54.20 

-24.10 

TOTAL 

31.19 

6.19 

-53.10 

-20.20 

ARG088 

32.21 

2.21 

-53.10 

-20.20 

BOL099 

36.69 

6.69 

-53.10 

-20.20 

CHL095 

51.87 

21.87 

-53.10 

-20.20 

PRU091 

51.75 

21.75 

-53.10 

-20.20 

URG096 

60.60 

30.60 

-58.10 

-20.20 

TOTAL 

30.31 

5.81 

-59.10 

-19.30 

ARG088 

31.99 

1.99 

-59.10 

-19.30 

BOL099 

35.81 

5.81 

-59.10 

-19.30 

CHL095 

51.46 

21.46 

-59.10 

-19.30 

PRU091 

51.35 

21.35 

-59.10 

-19.30 

URG096 

61.11 

31.11 

-59.10 

-19.30 

TOTAL 

30.41 

5.41 

-62.20 

-20.50 

ARG088 

31.31 

1.81 

-62.20 

-20.50 

BOL099 

35.22 

5.22 

-62.20 

-20.50 

CHL095 

45.90 

15.90 

-62.20 

-20.50 

PRU091 

51.45 

21 .45 
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-S:\20 

-62.20 

-V0.50 

-20.50 

URG096 

TOTAL 

60.07 

30.02 

30.07 

5.02 

-62.70 

-62.70 

-62.70 

-62.70 

-62.70 

-62.70 

-22.20 

-22.20 

-22.20 

-22.20 

-22.20 

-22.20 

ARG088 

BOL099 

CHL095 

PRU091 

URG096 

TOTAL 

31.03 
35 . 50 
38.46 
51.18 
57.93 
29.13 

1.03 
5.50 
8 . 46 
21.18 
27.93 
4.13 

-53.70 

-53.70 

-53.70 

-53.70 

-53.70 

-53.70 

-27.20 

-27.20 

-27.20 

-27.20 

-27.20 

-27.20 

ARG088 

BOL099 

CHL095 

PRU091 

URG096 

TOTAL 

30.12 
43.01 
51 .02 
50.90 
35.62 
28.32 

0. 12 
13.01 
21 .02 
20.90 
5.62 
3.82 


TEST COUNTRY s PRU 


TEST 

POINT 

-70.40 

-18.30 

-70.40 

-18.30 

-70.40 

-18.30 

-70.40 

-18.30 

-70.40 

-18.30 

-70.40 

-18.30 

-69.00 

-12.30 

-69.00 

-12.30 
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